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 F5 I97 R1 275 44.0 1733530 96.5% 1672856 122.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,545,596 T→C E181G (GAG→GGG)  yjiC ← YjfZ/YjiC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,545,5960TC96.6% 211.4 / ‑3.4 59E181G (GAG→GGG) yjiCYjfZ/YjiC family protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (27/30);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATC  >  NZ_CP009273/4545460‑4545726
                                                                                                                                        |                                                                                                                                  
ggCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCa                                                                                                                                  >  2:135363/1‑139 (MQ=255)
ggCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCa                                                                                                                                  >  2:503304/1‑139 (MQ=255)
  cTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCAtt                                                                                                                                >  1:140556/1‑139 (MQ=255)
      tataCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTc                                                                                                                            >  2:764011/1‑139 (MQ=255)
                attGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTg                                                                                                                  >  2:835666/1‑139 (MQ=255)
                     gTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCt                                                                                                               >  2:406197/1‑137 (MQ=255)
                        aTTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgag                                                                                                            >  1:863905/1‑137 (MQ=255)
                             ccAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTgg                                                                                                     <  2:863905/139‑1 (MQ=255)
                              cAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGa                                                                                                    <  2:35626/139‑1 (MQ=255)
                                    aaTATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTg                                                                                                                               >  1:699656/1‑106 (MQ=255)
                                    aaTATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTg                                                                                                                               <  2:699656/106‑1 (MQ=255)
                                    aaTATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGtt                                                                                              <  2:412246/139‑1 (MQ=255)
                                    aaTATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGt                                                                                               <  2:421652/138‑1 (MQ=255)
                                     atatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTg                                                                                             <  1:503304/139‑1 (MQ=255)
                                            cATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAatata                                                                                      >  2:792660/1‑139 (MQ=255)
                                              ttAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                    <  1:782562/139‑1 (MQ=255)
                                                         ggCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                         <  1:221294/139‑1 (MQ=255)
                                                               ggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCa                                                                   <  2:140556/139‑1 (MQ=255)
                                                               ggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCa                                                                   <  1:135363/139‑1 (MQ=255)
                                                                gCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgag                                                                                                            <  1:10502/97‑1 (MQ=255)
                                                                gCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgag                                                                                                            >  2:10502/1‑97 (MQ=255)
                                                                       cGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                                                           <  1:835666/139‑1 (MQ=255)
                                                                                   cTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAcc                                               <  1:771165/139‑1 (MQ=255)
                                                                                                cGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGa                                                                                                    <  1:750902/73‑1 (MQ=255)
                                                                                                cGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGa                                                                                                    >  2:750902/1‑73 (MQ=255)
                                                                                                  tCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAt                                >  2:582046/1‑139 (MQ=255)
                                                                                                  tCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAt                                >  2:150718/1‑139 (MQ=255)
                                                                                                    aTTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAttt                              >  1:452721/1‑139 (MQ=255)
                                                                                                     ttATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTc                             <  1:609201/139‑1 (MQ=255)
                                                                                                        tCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCAt                           <  1:150718/138‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCGATCGTGGTTCGGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        <  1:830803/137‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                         <  1:609174/88‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                         >  2:609174/1‑88 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        <  1:406197/137‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        <  1:582046/137‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        <  1:792660/137‑1 (MQ=255)
                                                                                                               aaTATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                   <  1:207890/139‑1 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCt            >  2:490622/1‑139 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgaga                                                                                                           <  2:266367/44‑1 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgaga                                                                                                           >  1:266367/1‑44 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                                           <  1:417304/60‑1 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                                           >  2:417304/1‑60 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg                                                    >  1:622013/1‑99 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg                                                    <  2:622013/99‑1 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                   >  1:775368/1‑132 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                   <  2:775368/132‑1 (MQ=255)
                                                                                                                        atCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAggg                                                                                                <  2:418886/53‑1 (MQ=255)
                                                                                                                        atCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAggg                                                                                                >  1:418886/1‑53 (MQ=255)
                                                                                                                         tCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        <  2:175989/124‑1 (MQ=255)
                                                                                                                         tCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        >  1:175989/1‑124 (MQ=255)
                                                                                                                           aaTCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAcc                                               <  1:401006/99‑1 (MQ=255)
                                                                                                                           aaTCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAcc                                               >  2:401006/1‑99 (MQ=255)
                                                                                                                           aaTCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCaa         >  1:247439/1‑137 (MQ=255)
                                                                                                                                ttGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATg                          <  2:494431/115‑1 (MQ=255)
                                                                                                                                ttGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATg                          >  1:494431/1‑115 (MQ=255)
                                                                                                                                ttGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                   <  1:628428/122‑1 (MQ=255)
                                                                                                                                ttGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                   >  2:628428/1‑122 (MQ=255)
                                                                                                                                ttGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATc  >  2:639748/1‑139 (MQ=255)
                                                                                                                                 tGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATc  >  1:691289/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                                  
GGCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATC  >  NZ_CP009273/4545460‑4545726

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACCGGGATAGGCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGC  >  NZ_CP009273/4545451‑4545734
                                                                                                                                                 |                                                                                                                                          
ACCGGGATAGGCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTG                                                                                                                                       >  SRR3722120.142457/1‑151 (MQ=60)
 agtctcgtgggctcggagatgtgtataagagacagACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGG                                                                                                                                      <  SRR3722120.161343/116‑1 (MQ=60)
                      ATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGA                                                                                                                   >  SRR3722120.875442/1‑149 (MQ=60)
                                cgtgggctcggagatgtgtataagagacagATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGG                                                                                                         <  SRR3722120.10630/119‑1 (MQ=60)
                                  TTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTctgtctcttatacacatctccg                                                                                                      >  SRR3722120.709205/1‑128 (MQ=60)
                                      acggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAAT                                                                                                 <  SRR3722120.761061/95‑1 (MQ=60)
                                             AATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCT                                                                                          <  SRR3722120.510505/151‑1 (MQ=60)
                                                ttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTAT                                                                                        <  SRR3722120.423426/82‑1 (MQ=60)
                                                      ATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTT                                                                                 <  SRR3722120.793104/151‑1 (MQ=60)
                                                                 cgccttagtctcgtgggctcggagatgtgtataagagacagGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACC                                                                      <  SRR3722120.617617/110‑1 (MQ=60)
                                                                 CGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACC                                                                      <  SRR3722120.224390/151‑1 (MQ=60)
                                                                       GGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAAC                                                                <  SRR3722120.137197/151‑1 (MQ=60)
                                                                               CCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAC                                                        <  SRR3722120.846924/151‑1 (MQ=60)
                                                                                           cgtgggctcggagatgtgtataagagacagATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATG                                            <  SRR3722120.406947/121‑1 (MQ=60)
                                                                                           GCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATG                                            <  SRR3722120.781585/151‑1 (MQ=60)
                                                                                                  CTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTC                                     >  SRR3722120.459322/1‑151 (MQ=60)
                                                                                                             ATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTAT                          <  SRR3722120.617644/151‑1 (MQ=60)
                                                                                                                ATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAG                        <  SRR3722120.152751/150‑1 (MQ=60)
                                                                                                                    GATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCGATCGTGGTTCGGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCC                     <  SRR3722120.841967/149‑1 (MQ=60)
                                                                                                                    GATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaacaaa                    >  SRR3722120.270192/1‑66 (MQ=60)
                                                                                                                    GATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCctgtctcttatacacatctccgagcccacg                   >  SRR3722120.630659/1‑121 (MQ=60)
                                                                                                                    GATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTC                   >  SRR3722120.785838/1‑151 (MQ=60)
                                                                                                                    GATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCC                     <  SRR3722120.412196/149‑1 (MQ=60)
                                                                                                                    GATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCC                     <  SRR3722120.590157/149‑1 (MQ=60)
                                                                                                                    GATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCC                     <  SRR3722120.803317/149‑1 (MQ=60)
                                                                                                                      TCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaa                  >  SRR3722120.425027/1‑76 (MQ=60)
                                                                                                                       gagacagGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAAC                <  SRR3722120.637132/144‑1 (MQ=60)
                                                                                                                       CAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCctgtc                >  SRR3722120.178388/1‑146 (MQ=60)
                                                                                                                       CAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAAC                <  SRR3722120.210759/151‑1 (MQ=60)
                                                                                                                         ATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAAC                >  SRR3722120.251027/1‑149 (MQ=60)
                                                                                                                              GATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCtgtctcttatac          >  SRR3722120.501547/1‑138 (MQ=60)
                                                                                                                               ATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCG         >  SRR3722120.700739/1‑150 (MQ=60)
                                                                                                                                     ATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaac  >  SRR3722120.718459/1‑74 (MQ=60)
                                                                                                                                     ATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctg  >  SRR3722120.848994/1‑89 (MQ=60)
                                                                                                                                                 |                                                                                                                                          
ACCGGGATAGGCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGC  >  NZ_CP009273/4545451‑4545734

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: