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
A1 F3 I150 R1 9 47.7 3300572 95.9% 3165248 84.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,754,231 G→T R127S (CGC→AGC)  selA ← L‑seryl‑tRNA(Sec) selenium transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,754,2310GT97.4% 299.0 / ‑2.1 76R127S (CGC→AGC) selAL‑seryl‑tRNA(Sec) selenium transferase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (34/40);  total (35/41)
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.99e-01

GCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754149‑3754316
                                                                                  |                                                                                     
gcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcg                                                                                >  1:473930/1‑90 (MQ=255)
gcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcg                                                                                >  2:1303060/1‑90 (MQ=255)
gcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcg                                                                                <  2:1228575/90‑1 (MQ=255)
 cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATCCGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:1629329/90‑1 (MQ=255)
 cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:639387/90‑1 (MQ=255)
 cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:686575/90‑1 (MQ=255)
 cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:419758/90‑1 (MQ=255)
 cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:163959/90‑1 (MQ=255)
 cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:1567688/90‑1 (MQ=255)
      ttATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:946174/90‑1 (MQ=255)
      ttATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:441123/90‑1 (MQ=255)
      ttATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  1:997617/90‑1 (MQ=255)
       tATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCa                                                                         >  1:382630/1‑90 (MQ=255)
         ttGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGa                                                                       <  1:550551/90‑1 (MQ=255)
          tGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                                      >  1:1160929/1‑90 (MQ=255)
          tGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                                      <  2:1298654/90‑1 (MQ=255)
            aCGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCa                                                                    <  2:1327763/90‑1 (MQ=255)
             cGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCAt                                                                   >  2:949167/1‑90 (MQ=255)
               aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               >  2:1149083/1‑76 (MQ=255)
               aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:1149083/76‑1 (MQ=255)
               aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCa                                                                         >  2:1534994/1‑82 (MQ=255)
               aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCa                                                                         <  1:1534994/82‑1 (MQ=255)
               aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  2:1182548/1‑90 (MQ=255)
               aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:386038/1‑90 (MQ=255)
               aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  2:438957/1‑90 (MQ=255)
               aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  2:956592/1‑90 (MQ=255)
               aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATAc                                                                 >  1:805873/1‑90 (MQ=255)
                   aGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGa                                                             <  1:297017/90‑1 (MQ=255)
                    ggCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAg                                                            <  1:1041657/90‑1 (MQ=255)
                    ggCATCTTCCACCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAg                                                            <  2:136086/90‑1 (MQ=255)
                       aTCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGt                                                         >  1:761677/1‑90 (MQ=255)
                        tCTTCCGCCCCCGTAATACGGCACACCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTc                                                        <  1:1297496/90‑1 (MQ=255)
                         cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                                      <  1:461352/75‑1 (MQ=255)
                         cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                                      >  2:461352/1‑75 (MQ=255)
                         cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCa                                                       >  1:345930/1‑90 (MQ=255)
                         cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCa                                                       >  1:713022/1‑90 (MQ=255)
                         cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCa                                                       >  2:645532/1‑90 (MQ=255)
                          ttCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                                      <  2:1629060/90‑1 (MQ=255)
                          ttCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                                      <  2:382630/90‑1 (MQ=255)
                              gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGc                                                  >  2:395958/1‑90 (MQ=255)
                              gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGCCACTGGc                                                  >  1:453182/1‑90 (MQ=255)
                                ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCGCTGGCGa                                                <  1:1164878/90‑1 (MQ=255)
                                ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:105547/90‑1 (MQ=255)
                                ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:440916/90‑1 (MQ=255)
                                ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:1209068/90‑1 (MQ=255)
                                ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:228213/90‑1 (MQ=255)
                                ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:187152/90‑1 (MQ=255)
                                ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:588334/90‑1 (MQ=255)
                                   cGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:418187/87‑1 (MQ=255)
                                   cGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                >  2:418187/1‑87 (MQ=255)
                                   cGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACg                                             >  2:999576/1‑90 (MQ=255)
                                         aCGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGc                                       >  1:1463808/1‑90 (MQ=255)
                                            gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:552650/52‑1 (MQ=255)
                                            gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          >  1:552650/1‑52 (MQ=255)
                                            gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:904685/78‑1 (MQ=255)
                                            gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                >  1:904685/1‑78 (MQ=255)
                                            gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg                                    >  2:366152/1‑90 (MQ=255)
                                            gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg                                    >  2:1219846/1‑90 (MQ=255)
                                              acaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:617452/45‑1 (MQ=255)
                                              acaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               >  2:617452/1‑45 (MQ=255)
                                               cagcagCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgcgc                                 >  1:1614309/1‑90 (MQ=255)
                                               cagcagCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgcgc                                 >  2:925289/1‑90 (MQ=255)
                                               cagcagCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgcgc                                 >  1:168251/1‑90 (MQ=255)
                                                 gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                               <  2:824979/90‑1 (MQ=255)
                                                 gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                               <  1:1057772/90‑1 (MQ=255)
                                                   agcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACg                             >  1:311109/1‑90 (MQ=255)
                                                      tGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt                          >  1:1166365/1‑90 (MQ=255)
                                                       gcgcCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgc                                   <  2:1029448/80‑1 (MQ=255)
                                                       gcgcCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgc                                   >  1:1029448/1‑80 (MQ=255)
                                                              cgcCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACcgc                  <  2:280993/90‑1 (MQ=255)
                                                                   gATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCtt             <  1:438957/90‑1 (MQ=255)
                                                                      cgcgATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCg          <  1:1594111/90‑1 (MQ=255)
                                                                       gcgATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGc         <  2:1076807/90‑1 (MQ=255)
                                                                            gTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                               <  1:298288/63‑1 (MQ=255)
                                                                            gTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                               >  2:298288/1‑63 (MQ=255)
                                                                              ccGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTaaa  >  2:90054/1‑90 (MQ=255)
                                                                                  |                                                                                     
GCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754149‑3754316

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGGGAGTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACGTGAT  >  NZ_CP009273/4267355‑4267552
                                                                                                   |                                                                                                  
GGGGAGTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAAT                                                                                                    >  SRR3722019.1639758/1‑100 (MQ=60)
  GGAGTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGC                                                                                                  >  SRR3722019.1641965/1‑100 (MQ=60)
  GGAGTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGC                                                                                                  <  SRR3722019.1219607/100‑1 (MQ=60)
  GGAGTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGC                                                                                                  >  SRR3722019.997592/1‑100 (MQ=60)
  GGAGTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGC                                                                                                  <  SRR3722019.962828/100‑1 (MQ=60)
   GAGTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTTTGGAAAAGAAATTACCCCGCATTAAAGCGTGGTTAACCCCCGGCGAGGGGGCGAAATGCA                                                                                                 >  SRR3722019.22704/1‑100 (MQ=60)
   GAGTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCA                                                                                                 >  SRR3722019.633766/1‑100 (MQ=60)
     GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGC                                                                                               >  SRR3722019.429230/1‑100 (MQ=60)
     GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGC                                                                                               >  SRR3722019.1573998/1‑100 (MQ=60)
     GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGC                                                                                               >  SRR3722019.612283/1‑100 (MQ=60)
     GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGC                                                                                               >  SRR3722019.115046/1‑100 (MQ=60)
     GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGC                                                                                               >  SRR3722019.11356/1‑100 (MQ=60)
     GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGC                                                                                               >  SRR3722019.1100057/1‑100 (MQ=60)
       ATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGG                                                                                             >  SRR3722019.500776/1‑100 (MQ=60)
       ATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGG                                                                                             <  SRR3722019.594065/100‑1 (MQ=60)
       ATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGG                                                                                             <  SRR3722019.1168168/100‑1 (MQ=60)
       ATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGG                                                                                             <  SRR3722019.62212/100‑1 (MQ=60)
               TCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGG                                                                                     <  SRR3722019.615549/100‑1 (MQ=60)
               TCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGG                                                                                     <  SRR3722019.1113183/100‑1 (MQ=60)
                       ACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGctgtctct                                                                             >  SRR3722019.1537719/1‑92 (MQ=60)
                        CTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGC                                                                            <  SRR3722019.489824/100‑1 (MQ=60)
                                  AGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTA                                                                  <  SRR3722019.1012818/100‑1 (MQ=60)
                                  AGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTA                                                                  <  SRR3722019.806044/100‑1 (MQ=60)
                                        TTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAG                                                            <  SRR3722019.1375684/100‑1 (MQ=60)
                                             GAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTTTGAAAGTAAAG                                                       >  SRR3722019.712744/1‑100 (MQ=60)
                                               AAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGG                                                     <  SRR3722019.1127828/100‑1 (MQ=60)
                                                cacggaactcctcgtcggcagcgtcagatgtgtataagagacagGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGT                                                    <  SRR3722019.1214841/56‑1 (MQ=60)
                                                AAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGT                                                    <  SRR3722019.873248/100‑1 (MQ=60)
                                                  GAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCcggatttcttttagaacaaaaggttg                                                  >  SRR3722019.315737/1‑74 (MQ=60)
                                                  GAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTTTGAAAGGAAAGGGGTG                                                  >  SRR3722019.591622/1‑100 (MQ=60)
                                                        ACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACC                                            <  SRR3722019.1492926/100‑1 (MQ=60)
                                                        ACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACC                                            <  SRR3722019.808656/100‑1 (MQ=60)
                                                          tgtgtataagagacagGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAG                                          <  SRR3722019.1553860/84‑1 (MQ=60)
                                                          CCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGAGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAG                                          <  SRR3722019.1204409/100‑1 (MQ=60)
                                                            CGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTA                                        <  SRR3722019.1040994/100‑1 (MQ=60)
                                                             GCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGAAAGTAAAGGGGTGATTACCAGTAT                                       >  SRR3722019.394930/1‑100 (MQ=60)
                                                             GCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTAT                                       >  SRR3722019.167002/1‑100 (MQ=60)
                                                             GCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGAAAAGGGGTGATTACCAGTAT                                       >  SRR3722019.6878/1‑100 (MQ=60)
                                                             GCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCCGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTAT                                       >  SRR3722019.1568224/1‑100 (MQ=60)
                                                              CATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCCTTTCTATGAAAATAAAGGGGTGATTACCctgtct                                      >  SRR3722019.1179473/1‑94 (MQ=60)
                                                                  AAAGCGCTGCTAACCCCCGGCGAAGTGGTGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTA                                  <  SRR3722019.380446/100‑1 (MQ=60)
                                                                  AAAGCGCTGCTAACCCCCGGCGAAGTGGCGAGATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTA                                  <  SRR3722019.294446/100‑1 (MQ=60)
                                                                    AGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAAC                                <  SRR3722019.985202/100‑1 (MQ=60)
                                                                    AGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAAC                                <  SRR3722019.526515/100‑1 (MQ=60)
                                                                    AGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAAC                                <  SRR3722019.1343302/100‑1 (MQ=60)
                                                                      CGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTAGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAG                              <  SRR3722019.83540/100‑1 (MQ=60)
                                                                        CTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCCTTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCG                            >  SRR3722019.142589/1‑100 (MQ=60)
                                                                        CTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTATCCGTAACAGCG                            >  SRR3722019.641257/1‑100 (MQ=60)
                                                                          GCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGGTTACCAGTATCCGTAACAGCGGC                          >  SRR3722019.276215/1‑100 (MQ=60)
                                                                               CCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGCATGGGCGCTGCAGTTTTATGAAAGAAAAGGGGTGATTACCAGTATCCGTAACctgtctcttat                     >  SRR3722019.1160404/1‑89 (MQ=60)
                                                                                        AAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCCTTTCTATGAAAATAAAGGGGTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATA            >  SRR3722019.433719/1‑100 (MQ=60)
                                                                                             GCGAAATGCAGCGGTGTGGCGGTATCGGCGTCGCATTTGTATGAAAGAAAAGGGGTGTTTACCAATATCCGTAAAAGAGGCAATCAGCGGCGATATAAAA       >  SRR3722019.1657172/1‑100 (MQ=50)
                                                                                             GCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGGTATAAAA       >  SRR3722019.253408/1‑100 (MQ=60)
                                                                                              CGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACG      <  SRR3722019.807855/100‑1 (MQ=60)
                                                                                              CGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACG      <  SRR3722019.1324985/100‑1 (MQ=60)
                                                                                              CGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACG      <  SRR3722019.453507/100‑1 (MQ=60)
                                                                                                  ATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTTTGAAAATAAAGGGTTGATTACCAGGATCCGTAAAAGCGGGAATCAGCGGCGATATAAAAGTGAT  >  SRR3722019.1189369/1‑100 (MQ=60)
                                                                                                  ATGCAGCGGTGTGGCGGGATCGGCGCTGCATTTCTATGAAAATAAAAGGGTGGTTTCCAGTATTCGTAACAGCGGCAAACAGCGGCGATATAAAAATGAT  >  SRR3722019.944821/1‑100 (MQ=56)
                                                                                                   |                                                                                                  
GGGGAGTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACGTGAT  >  NZ_CP009273/4267355‑4267552

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: