Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,533,061 T→G 53.3% K124N (AAA→AAC ftsE ← predicted transporter subunit: ATP‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,533,0610TG53.3% 4.0 / 93.1 60K124N (AAA→AACftsEpredicted transporter subunit: ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (13/15);  new base G (13/19);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

CCGACAACGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTG  >  REL606/3532913‑3533195
                                                                                                                                                    |                                                                                                                                      
ccGACAACGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACtt                                                                                                                                       <  1:649102/150‑1 (MQ=255)
     aaCGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATcc                                                                                                                                  <  2:601920/150‑1 (MQ=255)
       cgcgTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAg                                                                                                                                <  2:310398/150‑1 (MQ=255)
          gtcgtcCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTgc                                                                                                                             <  2:362638/150‑1 (MQ=255)
           tcgtcCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTgcc                                                                                                                            >  2:685268/1‑150 (MQ=255)
             gtcCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTgccgc                                                                                                                          >  2:221942/1‑150 (MQ=255)
                    ttACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACAcc                                                                                                                   <  2:311297/150‑1 (MQ=255)
                              gTTCGTCCGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGGAGAGGGGAGTTCTTCGCTTTGTTCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGGCGaa                                                                                                         >  1:305660/1‑150 (MQ=255)
                                     cGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAACCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGATATCCAGTGCCGCCGACACCCGGCGACGAATAGCGTc                                                                                                  >  1:516114/1‑150 (MQ=255)
                                          gcagTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTATCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGc                                                                                             >  1:202056/1‑150 (MQ=255)
                                                     ggCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCg                                                                                  >  1:57850/1‑150 (MQ=255)
                                                          aTTCCCCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGCAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACCAATATCGTCACCGATAGCACCAGCGATAAt                                                                             >  1:69822/1‑150 (MQ=255)
                                                          aTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAAt                                                                             >  2:513865/1‑150 (MQ=255)
                                                          aTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAAt                                                                             >  1:290627/1‑150 (MQ=255)
                                                          aTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAAt                                                                             >  2:670938/1‑150 (MQ=255)
                                                          aTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCGGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATGCAGTTCCGCCGACTCCCGGCGACGACTATCGTCCCGGCTGGCACCAGCCATAAt                                                                             >  1:501116/1‑150 (MQ=255)
                                                           ttCAGCACCGCCCGGGCCATGCCAACACGCTGTTGTTCACCGCCCGAAAGCCGGATAGGGAAGTTCTTCGCTGTGTCCAGTAGCCCGTCGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATc                                                                            <  1:289910/150‑1 (MQ=255)
                                                           ttCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATc                                                                            <  2:58019/150‑1 (MQ=255)
                                                                  ccGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGa                                                                     >  2:274150/1‑150 (MQ=255)
                                                                         ggCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCAc                                                              >  2:477627/1‑150 (MQ=255)
                                                                              tGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTAt                                                         >  2:100391/1‑150 (MQ=255)
                                                                                acagcacGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCg                                                       <  1:222506/146‑1 (MQ=255)
                                                                                    cacGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCgctgt                                                   >  2:223143/1‑146 (MQ=255)
                                                                                     gacaggttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  1:427611/145‑1 (MQ=255)
                                                                                     gacaggttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  2:213675/145‑1 (MQ=255)
                                                                                     gacaggttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  1:349670/145‑1 (MQ=255)
                                                                                     ccGCTGTTGTTCACCGCCCGAAAGCTGAATATGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGCATATCGTCACCGCTGGCACCAGCGATAATCAGTGGGATCGCCACGTTATCGTAGAc                                                  <  1:365926/149‑1 (MQ=255)
                                                                                     cagCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  2:478086/148‑1 (MQ=255)
                                                                                     acGGTGTTGTTCACCGCCGGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCTGTAGCCCGCCTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  1:261571/150‑1 (MQ=255)
                                                                                     acGCTGTTGTTGACAGCCCGAAAGCTGAATAGGGAAGTTCTTCGCCTTGTCAAGTAGCCCGACTTCATCCAGTGCCGCCGACACCCGGCGACGAATAGCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  1:162335/150‑1 (MQ=255)
                                                                                     acGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACCCCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  1:388356/150‑1 (MQ=255)
                                                                                     acGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  2:291454/150‑1 (MQ=255)
                                                                                     acGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  2:374683/150‑1 (MQ=255)
                                                                                     acGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  2:355248/150‑1 (MQ=255)
                                                                                     acGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  2:581009/150‑1 (MQ=255)
                                                                                     acGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  1:100095/150‑1 (MQ=255)
                                                                                     acGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  2:202641/150‑1 (MQ=255)
                                                                                     acGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  2:517594/150‑1 (MQ=255)
                                                                                     acGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  1:476276/150‑1 (MQ=255)
                                                                                     acGCGGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGATTTGTCCAGTAGCCCGGCGTTATCCAGTGCCGCCGACACCCGGCGCCGAATATCGTCACCGCCGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc                                                  <  1:447024/150‑1 (MQ=255)
                                                                                        cTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCCCCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGt                                               >  1:340829/1‑150 (MQ=255)
                                                                                        cTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAcctg                                               >  1:476733/1‑147 (MQ=255)
                                                                                        cTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGt                                               >  2:287133/1‑150 (MQ=255)
                                                                                          gttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAcctgtc                                             >  2:428826/1‑145 (MQ=255)
                                                                                          gttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTAc                                             <  2:341805/150‑1 (MQ=255)
                                                                                          gttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAcctgtc                                             >  1:213063/1‑145 (MQ=255)
                                                                                          gttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAcctgtc                                             >  2:350671/1‑145 (MQ=255)
                                                                                          gttgttCACCGCCCGAAAGCTGAATAGGGAAGCTCTTCGCTTTGTCCAGTAGCCCGACGTTACCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCCGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGGAc                                             <  1:536254/150‑1 (MQ=255)
                                                                                             gttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGt                                          <  2:70024/150‑1 (MQ=255)
                                                                                               tCACCGCCCGAAAGCTGAATAGGGTAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGCCGACGATTATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTcc                                        <  1:243148/150‑1 (MQ=255)
                                                                                               tCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTcc                                        <  2:306523/150‑1 (MQ=255)
                                                                                                     cccGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCTAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGt                                  >  2:496339/1‑150 (MQ=255)
                                                                                                     cccGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGt                                  <  2:362238/150‑1 (MQ=255)
                                                                                                       cGAAAGCTGAATAGGGAAGTTCTTCGCTGTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTTCATCAGTAg                                <  1:49424/150‑1 (MQ=255)
                                                                                                         aaaGCGGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGAt                              <  2:494608/150‑1 (MQ=255)
                                                                                                           ggccgAATAGGGAGGTTCTTGCCCTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCGGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATgg                            <  1:683261/146‑1 (MQ=255)
                                                                                                               gAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGAt                        <  1:668965/150‑1 (MQ=255)
                                                                                                               gAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGAt                        >  2:284862/1‑150 (MQ=255)
                                                                                                                 aTAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATcc                      >  1:298233/1‑150 (MQ=255)
                                                                                                                           tcttcGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCa            >  2:13639/1‑150 (MQ=255)
                                                                                                                                     gTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTg  >  2:517316/1‑150 (MQ=255)
                                                                                                                                     gTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTg  >  2:651265/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                      
CCGACAACGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTG  >  REL606/3532913‑3533195

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: