Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A2 F5 I175 R1 328 64.7 2582366 96.2% 2484236 120.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 734,002 T→C pseudogene (455/759 nt) ybfD → ISAs1 family transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273734,0020TC93.2% 174.5 / 3.7 59pseudogene (455/759 nt)ybfDISAs1 family transposase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (29/26);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAG  >  NZ_CP009273/733870‑734136
                                                                                                                                    |                                                                                                                                      
ttGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGc                                                                                                                                  <  1:248393/139‑1 (MQ=255)
        aaaaTACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCt                                                                                                                           >  1:1136116/1‑138 (MQ=255)
        aaaaTACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCt                                                                                                                           >  2:277742/1‑138 (MQ=255)
        aaaaTACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTg                                                                                                                          <  1:296178/139‑1 (MQ=255)
          aaTACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGaa                                                                                                                        <  1:205251/139‑1 (MQ=255)
                 cGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAActtctt                                                                                                                  >  1:1285327/1‑138 (MQ=255)
                   gCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAAc                                                                                                               <  1:724020/139‑1 (MQ=255)
                    cATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACa                                                                                                              <  1:870180/139‑1 (MQ=255)
                    cATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACa                                                                                                              >  2:757671/1‑139 (MQ=255)
                      ttCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATg                                                                                                            <  2:1077972/139‑1 (MQ=38)
                      ttCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATg                                                                                                            >  1:682941/1‑139 (MQ=38)
                      ttCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATg                                                                                                            <  2:921606/139‑1 (MQ=38)
                           aTGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCt                                                                                                          <  1:277742/136‑1 (MQ=37)
                           aTGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCt                                                                                                          <  2:682941/136‑1 (MQ=37)
                              aCAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGatat                                                                                                    >  2:828107/1‑139 (MQ=37)
                                       gCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGaa                                                                                            >  2:182025/1‑138 (MQ=37)
                                         cGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGaaaaa                                                                                         >  1:600247/1‑139 (MQ=37)
                                                   cGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATg                                                                                                            <  2:89388/110‑1 (MQ=32)
                                                   cGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATg                                                                                                            >  1:89388/1‑110 (MQ=32)
                                                    gATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAAcc                                                                              <  2:568866/139‑1 (MQ=25)
                                                     aTTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACcg                                                                             <  1:744218/139‑1 (MQ=37)
                                                         aTGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACcgatgc                                                                         <  2:393604/139‑1 (MQ=37)
                                                            tCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACcgatgcga                                                                       <  2:1136116/138‑1 (MQ=37)
                                                                     cGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGc                                                               >  2:1072717/1‑137 (MQ=37)
                                                                       ttCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAg                                                            >  1:176128/1‑138 (MQ=37)
                                                                        tctcAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGa                                                           >  1:356108/1‑138 (MQ=37)
                                                                         ctcAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGaa                                                          <  2:718255/138‑1 (MQ=37)
                                                                                 tGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATggg                                                                   <  1:1129925/121‑1 (MQ=38)
                                                                                 tGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATggg                                                                   >  2:1129925/1‑121 (MQ=38)
                                                                                 tGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTg                                                  >  1:1166109/1‑138 (MQ=37)
                                                                                     caGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCagag                                             <  1:1264322/139‑1 (MQ=37)
                                                                                      aGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCagag                                             >  1:1139788/1‑138 (MQ=37)
                                                                                        tCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCagag                                             >  1:53155/1‑136 (MQ=37)
                                                                                             tCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACcg                                                                             <  2:995341/99‑1 (MQ=255)
                                                                                             tCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACcg                                                                             >  1:995341/1‑99 (MQ=255)
                                                                                                 cGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACaaaaa                                 >  1:741400/1‑139 (MQ=37)
                                                                                                  ggACAGATCAAGCCGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACaaaaa                                 >  1:450917/1‑138 (MQ=255)
                                                                                                  ggACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCt                                                                                                                           <  2:974594/48‑1 (MQ=255)
                                                                                                  ggACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCt                                                                                                                           >  1:974594/1‑48 (MQ=255)
                                                                                                  ggACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACaaaaa                                 >  2:715682/1‑138 (MQ=37)
                                                                                                  ggACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACaaaaa                                 <  2:1067548/138‑1 (MQ=37)
                                                                                                     cAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAgg                             <  2:1285327/139‑1 (MQ=21)
                                                                                                      aGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAgg                             >  2:566287/1‑138 (MQ=21)
                                                                                                       gATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGa                           <  1:865668/139‑1 (MQ=34)
                                                                                                       gATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAgg                             >  1:931961/1‑137 (MQ=21)
                                                                                                                  aTAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACaaaaa                                 <  1:636647/122‑1 (MQ=25)
                                                                                                                  aTAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACaaaaa                                 >  2:636647/1‑122 (MQ=25)
                                                                                                                  aTAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGattatt                 >  2:362729/1‑138 (MQ=21)
                                                                                                                  aTAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGattatt                 >  2:213352/1‑138 (MQ=21)
                                                                                                                  aTAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGattatt                 >  2:1003222/1‑138 (MQ=21)
                                                                                                                  aTAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGattatt                 >  1:521354/1‑138 (MQ=21)
                                                                                                                  aTAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGattatt                 >  2:722660/1‑138 (MQ=21)
                                                                                                                     aGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATtatttatt             <  2:600247/139‑1 (MQ=21)
                                                                                                                     aGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATtatttatt             <  2:741400/139‑1 (MQ=21)
                                                                                                                     aGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATtatttatt             <  1:722660/139‑1 (MQ=21)
                                                                                                                            cAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGt       <  2:1139788/138‑1 (MQ=21)
                                                                                                                              aTGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTaa     >  1:144246/1‑138 (MQ=34)
                                                                                                                               tGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTaa     >  1:329453/1‑137 (MQ=21)
                                                                                                                                 agaCCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAg  <  2:1161668/138‑1 (MQ=21)
                                                                                                                                    |                                                                                                                                      
TTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAG  >  NZ_CP009273/733870‑734136

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCT  >  NZ_CP009273/733867‑734152
                                                                                                                                       |                                                                                                                                                      
CAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTG                                                                                                                                          >  SRR3722129.1151930/1‑150 (MQ=60)
  ATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAAC                                                                                                                                       <  SRR3722129.251909/151‑1 (MQ=60)
         GAAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTA                                                                                                                                 >  SRR3722129.1303137/1‑150 (MQ=60)
          AAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAAC                                                                                                                               <  SRR3722129.300417/151‑1 (MQ=60)
            AAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACAT                                                                                                                             <  SRR3722129.208117/151‑1 (MQ=60)
              ATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGC                                                                                                                           >  SRR3722129.693022/1‑151 (MQ=52)
                     GGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATAT                                                                                                                    <  SRR3722129.734629/151‑1 (MQ=42)
                      GCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATT                                                                                                                   <  SRR3722129.882601/151‑1 (MQ=40)
                           CTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCctgtctcttatacacatctccgagcccacgagact                                                                                                               >  SRR3722129.1222482/1‑115 (MQ=51)
                             TATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAG                                                                                                               <  SRR3722129.281729/148‑1 (MQ=32)
                                 ACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAAT                                                                                                        >  SRR3722129.609203/1‑151 (MQ=31)
                                           CCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAActgtctcttatacacatct                                                                                              >  SRR3722129.90555/1‑132 (MQ=25)
                                                       GATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGT                                                                                  <  SRR3722129.755099/151‑1 (MQ=31)
                                                               TCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGA                                                                           >  SRR3722129.178545/1‑150 (MQ=32)
                                                                CATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAA                                                                          >  SRR3722129.361247/1‑150 (MQ=32)
                                                                 agagacagGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAG                                                                        <  SRR3722129.1145664/143‑1 (MQ=32)
                                                                         GTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgtaaaaaaaaaaaaca                                                                 >  SRR3722129.268439/1‑69 (MQ=11)
                                                                         GTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGC                                                                 >  SRR3722129.1182365/1‑150 (MQ=32)
                                                                              CAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGA                                                            >  SRR3722129.1155643/1‑150 (MQ=32)
                                                                                ACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGA                                                            >  SRR3722129.53861/1‑148 (MQ=32)
                                                                                     GCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTctgtctcttatacacatctccgagcccacg                                                    >  SRR3722129.1009378/1‑121 (MQ=34)
                                                                                       ACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAA                                                  <  SRR3722129.1281900/151‑1 (MQ=31)
                                                                                         AGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAAC                                                >  SRR3722129.752247/1‑151 (MQ=31)
                                                                                          GTCTGGTCATCGGACAGATCAAGCCGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAAC                                                >  SRR3722129.457436/1‑150 (MQ=31)
                                                                                          GTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACtgtctcttatacacatctccgagcccacgagactaaggcgaatatcgtatgccgtcttctgcttgaaaaaaaaaaaaac                                                >  SRR3722129.988392/1‑71 (MQ=11)
                                                                                               GTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGG                                            >  SRR3722129.945183/1‑149 (MQ=32)
                                                                                                   gagacagGATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGA                                      <  SRR3722129.646107/144‑1 (MQ=21)
                                                                                                         AGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTT                                <  SRR3722129.877999/151‑1 (MQ=22)
                                                                                                          GATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAActgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaa                                >  SRR3722129.292686/1‑69 (MQ=0)
                                                                                                          GATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATctgtctcttatacacatctccgacc                               >  SRR3722129.1058199/1‑126 (MQ=22)
                                                                                                          GATCAAGACGGATAAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTT                                >  SRR3722129.528970/1‑150 (MQ=21)
                                                                                                                      TAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAA                    >  SRR3722129.146266/1‑150 (MQ=22)
                                                                                                                       AAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAG                  <  SRR3722129.733252/151‑1 (MQ=21)
                                                                                                                       AAGAAATCCAATGAGACCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAA                    >  SRR3722129.334201/1‑149 (MQ=21)
                                                                                                                                       CCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCT  >  SRR3722129.295179/1‑151 (MQ=21)
                                                                                                                                       |                                                                                                                                                      
CAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCT  >  NZ_CP009273/733867‑734152

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: