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 I160 R1 220 104.1 6051714 95.9% 5803593 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,117,071 T→C F44S (TTC→TCC)  rpmE → 50S ribosomal protein L31

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,117,0710TC100.0% 388.2 / NA 102F44S (TTC→TCC) rpmE50S ribosomal protein L31
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (48/54);  total (48/54)

Only 100 of 108 total aligned reads displayed.

GCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGAC  >  NZ_CP009273/4116984‑4117156
                                                                                       |                                                                                     
 ctgctCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAc                                                                                    <  2:1515664/90‑1 (MQ=255)
 ctgctCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAc                                                                                    <  1:77161/90‑1 (MQ=255)
 ctgctCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAc                                                                                    <  1:2680781/90‑1 (MQ=255)
 ctgctCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAc                                                                                    <  2:525882/90‑1 (MQ=255)
   gctCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTg                                                                                  >  1:2475049/1‑90 (MQ=255)
      ctTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCa                                                                               <  2:51647/90‑1 (MQ=255)
       ttGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCaa                                                                              <  1:310314/90‑1 (MQ=255)
         gCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAAc                                                                            <  2:2680246/90‑1 (MQ=255)
         gCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAAc                                                                            <  2:2761290/90‑1 (MQ=255)
              aaCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGt                                                                       >  1:2492195/1‑90 (MQ=255)
              aaCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGt                                                                       >  2:2968446/1‑90 (MQ=255)
               aCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTg                                                                      >  1:2533199/1‑90 (MQ=255)
               aCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTg                                                                      <  2:1922668/90‑1 (MQ=255)
                cGTAATGAAAATCCGCTCCACCGTTGGTCATGCCCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGa                                                                     <  1:1335086/90‑1 (MQ=255)
                cGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGa                                                                     <  1:422236/90‑1 (MQ=255)
                cGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGa                                                                     <  1:2738352/90‑1 (MQ=255)
                cGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGa                                                                     <  2:2212520/90‑1 (MQ=255)
                  tAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATg                                                                   <  1:1173857/90‑1 (MQ=255)
                  tAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATg                                                                   <  2:181978/90‑1 (MQ=255)
                  tAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATg                                                                   <  2:2843559/90‑1 (MQ=255)
                      gAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGc                                                               <  2:2492195/90‑1 (MQ=255)
                       aaaaTCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCt                                                              <  1:2538996/90‑1 (MQ=255)
                        aaaTCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCGAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTa                                                             <  2:1831171/90‑1 (MQ=255)
                        aaaTCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAAc                                                                            <  1:53733/75‑1 (MQ=255)
                        aaaTCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAAc                                                                            >  2:53733/1‑75 (MQ=255)
                          aTCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTAcc                                                           >  1:27864/1‑90 (MQ=255)
                           tCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGt                                                                  >  1:1309514/1‑82 (MQ=255)
                           tCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGt                                                                  <  2:1309514/82‑1 (MQ=255)
                           tCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCg                                                          <  1:2917224/90‑1 (MQ=255)
                           tCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCg                                                          <  2:2787571/90‑1 (MQ=255)
                              gCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTg                                                       >  1:1447441/1‑90 (MQ=255)
                              gCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTg                                                       >  2:2739623/1‑90 (MQ=255)
                                tCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGc                                                     >  2:1590439/1‑90 (MQ=255)
                                     cGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGc                                                     <  2:2284253/85‑1 (MQ=255)
                                     cGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGc                                                     >  1:2284253/1‑85 (MQ=255)
                                        tGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGa                                             >  2:611958/1‑90 (MQ=255)
                                        tGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGa                                             >  2:278972/1‑90 (MQ=255)
                                         ggTCATGACCTGAACTTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGAc                                            <  1:370564/90‑1 (MQ=255)
                                             aTGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCt                                        >  2:1737021/1‑90 (MQ=255)
                                             aTGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCt                                        >  1:2842940/1‑90 (MQ=255)
                                              tGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCtt                                       >  1:2693753/1‑90 (MQ=255)
                                              tGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCtt                                       >  2:518479/1‑90 (MQ=255)
                                              tGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCtt                                       >  1:2441311/1‑90 (MQ=255)
                                              tGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCtt                                       >  1:2458098/1‑90 (MQ=255)
                                               gACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTc                                      <  2:152137/90‑1 (MQ=255)
                                                   tGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAc                                                                                    >  1:1608372/1‑40 (MQ=255)
                                                   tGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAc                                                                                    <  2:1608372/40‑1 (MQ=255)
                                                   tGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATg                                                                   <  2:1592020/57‑1 (MQ=255)
                                                   tGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAATGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaaca                                  >  2:2257470/1‑90 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCGCCCGTTCTCCACTTGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 <  1:1487579/90‑1 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGc                                                               <  1:1196701/60‑1 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGc                                                               >  2:1196701/1‑60 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 >  1:1058992/1‑90 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 >  2:2357960/1‑90 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 <  2:27864/90‑1 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 <  2:316367/90‑1 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 >  2:307261/1‑90 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 <  2:309858/90‑1 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 <  1:1857272/90‑1 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 >  2:120040/1‑90 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 <  1:2325255/90‑1 (MQ=255)
                                                    gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                 >  1:1058364/1‑90 (MQ=255)
                                                      aCCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGc                               <  2:1336509/90‑1 (MQ=255)
                                                      aCCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGc                               <  2:2714845/90‑1 (MQ=255)
                                                        cTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGt                             <  2:504235/90‑1 (MQ=255)
                                                        cTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGt                             <  1:635464/90‑1 (MQ=255)
                                                         tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt                            <  1:541347/90‑1 (MQ=255)
                                                         tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt                            >  1:1693787/1‑90 (MQ=255)
                                                         tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt                            >  1:1708387/1‑90 (MQ=255)
                                                         tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt                            >  1:2848864/1‑90 (MQ=255)
                                                         tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt                            >  2:13620/1‑90 (MQ=255)
                                                           gACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGTCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTc                          <  2:1805180/90‑1 (MQ=255)
                                                             cGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCaa                        <  1:611958/90‑1 (MQ=255)
                                                             cGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCaa                        <  2:2408237/90‑1 (MQ=255)
                                                              gtgtGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAAc                       >  1:977950/1‑90 (MQ=255)
                                                                gtgCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTg                                                                      >  1:2423687/1‑41 (MQ=255)
                                                                gtgCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTg                                                                      <  2:2423687/41‑1 (MQ=255)
                                                                gtgCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACAt                     >  1:1520880/1‑90 (MQ=255)
                                                                gtgCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACAt                     >  2:2794511/1‑90 (MQ=255)
                                                                gtgCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACAt                     >  2:1452680/1‑90 (MQ=255)
                                                                gtgCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACAt                     >  2:2275511/1‑90 (MQ=255)
                                                                  gcagcaAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATcc                   >  2:3004964/1‑90 (MQ=255)
                                                                    agcaAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCg                 <  1:1512519/90‑1 (MQ=255)
                                                                     gcaAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCgg                >  2:1552620/1‑90 (MQ=255)
                                                                       aaGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgc              <  1:2243263/90‑1 (MQ=255)
                                                                       aaGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgc              <  1:2118876/90‑1 (MQ=255)
                                                                       aaGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgc              >  2:1209313/1‑90 (MQ=255)
                                                                         gTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcag            <  2:1885041/90‑1 (MQ=255)
                                                                          tGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcagc           >  1:754709/1‑90 (MQ=255)
                                                                          tGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcagc           >  1:2117374/1‑90 (MQ=255)
                                                                          tGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcagc           >  1:2095449/1‑90 (MQ=255)
                                                                          tGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcagc           <  2:1702481/90‑1 (MQ=255)
                                                                            ccACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCaa         <  1:2488185/90‑1 (MQ=255)
                                                                                 cGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAg    <  1:1351992/90‑1 (MQ=255)
                                                                                   ttctccACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGAc  >  1:656808/1‑90 (MQ=255)
                                                                                   ttctccACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGAc  >  1:1319477/1‑90 (MQ=255)
                                                                                   ttctccACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGAc  >  2:1883118/1‑90 (MQ=255)
                                                                                   gtctccACTGGCAAACAGCGGGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGAc  >  1:1529631/2‑90 (MQ=255)
                                                                                    tctccACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaaca                                  <  1:668667/57‑1 (MQ=255)
                                                                                    tctccACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaaca                                  >  2:668667/1‑57 (MQ=255)
                                                                                       |                                                                                     
GCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGAC  >  NZ_CP009273/4116984‑4117156

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 32 ≤ 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
AAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAA  >  NZ_CP009273/4116972‑4117170
                                                                                                   |                                                                                                   
AAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTC                                                                                                     >  SRR3722030.824443/1‑100 (MQ=60)
AAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTC                                                                                                     >  SRR3722030.468183/1‑100 (MQ=60)
AAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTC                                                                                                     <  SRR3722030.2427048/100‑1 (MQ=60)
 AATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCC                                                                                                    >  SRR3722030.1375998/1‑100 (MQ=60)
   cggcgaccaccgagatctacactggccatgtcgtcggcagcgtcagatgtgtataagagacagACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAC                                                                                                  <  SRR3722030.2699408/37‑1 (MQ=60)
   TTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAC                                                                                                  <  SRR3722030.228132/100‑1 (MQ=60)
   TTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAC                                                                                                  <  SRR3722030.1940834/100‑1 (MQ=60)
   TTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCAC                                                                                                  <  SRR3722030.1396343/100‑1 (MQ=60)
    TACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACT                                                                                                 <  SRR3722030.1731438/100‑1 (MQ=60)
     ACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTG                                                                                                >  SRR3722030.2520268/1‑100 (MQ=60)
        GCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCA                                                                                             <  SRR3722030.2632005/100‑1 (MQ=60)
        GCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCA                                                                                             <  SRR3722030.2809811/100‑1 (MQ=60)
         CTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAA                                                                                            <  SRR3722030.1432672/100‑1 (MQ=60)
             tacactggccatgtcgtcggcagcgtcagatgtgtataagagacagGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAG                                                                                        <  SRR3722030.174371/54‑1 (MQ=60)
             cggcgaccaccgagatctacactggccatgtcgtcggcagcgtcagatgtgtataagagacagGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAG                                                                                        <  SRR3722030.2243427/37‑1 (MQ=60)
             CTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAG                                                                                        <  SRR3722030.78356/100‑1 (MQ=60)
             CTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAG                                                                                        <  SRR3722030.2730233/100‑1 (MQ=60)
                CTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGT                                                                                     >  SRR3722030.2537778/1‑100 (MQ=60)
                 TCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTG                                                                                    >  SRR3722030.2579593/1‑100 (MQ=60)
                   TTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGAT                                                                                  <  SRR3722030.314833/100‑1 (MQ=60)
                     gacagAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGT                                                                                <  SRR3722030.54559/95‑1 (MQ=60)
                            CGTAATGAAAATCCGCTCCACCGTTGGTCATGCCCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACC                                                                         <  SRR3722030.1355654/100‑1 (MQ=60)
                            CGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACC                                                                         <  SRR3722030.2788928/100‑1 (MQ=60)
                            CGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACC                                                                         <  SRR3722030.428626/100‑1 (MQ=60)
                            CGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACC                                                                         >  SRR3722030.28292/1‑100 (MQ=60)
                             GTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCG                                                                        >  SRR3722030.1329473/1‑100 (MQ=60)
                              TAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGG                                                                       <  SRR3722030.1191261/100‑1 (MQ=60)
                                ATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTG                                                                     >  SRR3722030.1470454/1‑100 (MQ=60)
                                  cagatgtgtataagagacaGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGC                                                                   <  SRR3722030.1214414/81‑1 (MQ=60)
                                   AAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCC                                                                  <  SRR3722030.2585500/100‑1 (MQ=60)
                                       TCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGT                                                              <  SRR3722030.2971019/100‑1 (MQ=60)
                                       TCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGT                                                              >  SRR3722030.2325619/1‑100 (MQ=60)
                                               ACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCT                                                      >  SRR3722030.2895449/1‑100 (MQ=60)
                                                CCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTT                                                     >  SRR3722030.2503022/1‑100 (MQ=60)
                                                CCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTT                                                     >  SRR3722030.2743477/1‑100 (MQ=60)
                                                CCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTT                                                     >  SRR3722030.2485856/1‑100 (MQ=60)
                                                     GGTCATGACCTGAACTTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACA                                                <  SRR3722030.376083/100‑1 (MQ=60)
                                                     GGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCtgtctcttatacacatctgacgctgccgacgacatggcc                                                >  SRR3722030.1634458/1‑61 (MQ=60)
                                                     GGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGctgtctcttatacacatctgacg                                                >  SRR3722030.1617841/1‑77 (MQ=60)
                                                      GTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAA                                               >  SRR3722030.314371/1‑100 (MQ=60)
                                                      GTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAA                                               >  SRR3722030.1074267/1‑100 (MQ=60)
                                                      GTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAA                                               >  SRR3722030.1074900/1‑100 (MQ=60)
                                                           GACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTT                                          >  SRR3722030.1736266/1‑100 (MQ=60)
                                                           GACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTT                                          >  SRR3722030.1721299/1‑100 (MQ=60)
                                                           GACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGGGGCCGGGTTGACCGCTTCAACAAGCGTT                                          >  SRR3722030.2901472/1‑100 (MQ=60)
                                                               cgtcagatgtgtataagagacAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAA                                      <  SRR3722030.679030/79‑1 (MQ=60)
                                                                tcggcagcgtcagatgtgtataagagacagGTTCTCCACTGGCAAACAGCGTGATATTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAAC                                     <  SRR3722030.845029/70‑1 (MQ=60)
                                                                GAACCTCGACGTGTGCAGCAAGTGCCGCCCGTTCTCCACTTGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAAC                                     <  SRR3722030.1511402/100‑1 (MQ=60)
                                                                GAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAAC                                     >  SRR3722030.992837/1‑100 (MQ=60)
                                                                GAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAAC                                     <  SRR3722030.2367460/100‑1 (MQ=60)
                                                                GAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAAC                                     <  SRR3722030.1889006/100‑1 (MQ=60)
                                                                  ACCTTGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACAT                                   >  SRR3722030.1545401/1‑100 (MQ=60)
                                                                  ACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCtgtctcttatacacatctgacgctgccgacgacatggc                                   >  SRR3722030.2467873/1‑62 (MQ=60)
                                                                    CTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCC                                 <  SRR3722030.645336/100‑1 (MQ=60)
                                                                     TCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCC                                <  SRR3722030.549700/100‑1 (MQ=60)
                                                                         CGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGC                            <  SRR3722030.621437/100‑1 (MQ=60)
                                                                           TGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAG                          <  SRR3722030.2497047/100‑1 (MQ=60)
                                                                            GTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCtgtctcttatacacatctgacgctgccgacgacatggccagtgtagatctcggtggtcgccg                         >  SRR3722030.73104/1‑38 (MQ=60)
                                                                            GTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGC                         >  SRR3722030.2155154/1‑100 (MQ=60)
                                                                            GTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGC                         >  SRR3722030.2132770/1‑100 (MQ=60)
                                                                            GTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGC                         >  SRR3722030.766336/1‑100 (MQ=60)
                                                                                AGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAAT                     <  SRR3722030.1536870/100‑1 (MQ=60)
                                                                                   AAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAG                  <  SRR3722030.2156681/100‑1 (MQ=60)
                                                                                   AAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAG                  <  SRR3722030.2283762/100‑1 (MQ=60)
                                                                                     GTGTCACCCGGTCTCCACTGGCAAACAGCGGGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGAC                >  SRR3722030.1554307/1‑100 (MQ=60)
                                                                                     GTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGAC                >  SRR3722030.1339680/1‑100 (MQ=60)
                                                                                     GTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGAC                >  SRR3722030.666978/1‑100 (MQ=60)
                                                                                        CCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGG             <  SRR3722030.2533693/100‑1 (MQ=60)
                                                                                             CGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTC        >  SRR3722030.2952200/1‑100 (MQ=60)
                                                                                             CGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTC        <  SRR3722030.1372981/100‑1 (MQ=60)
                                                                                              GTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGATTCAACATCCCGGGCAGCAAATAAGACCGGATTTCC       >  SRR3722030.1644255/1‑100 (MQ=60)
                                                                                                 CTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCtgtctcttata    >  SRR3722030.2672015/1‑89 (MQ=60)
                                                                                                 CTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAA    >  SRR3722030.779779/1‑100 (MQ=60)
                                                                                                 CTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAA    >  SRR3722030.938517/1‑100 (MQ=60)
                                                                                                   CCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAA  >  SRR3722030.1395216/1‑100 (MQ=60)
                                                                                                   |                                                                                                   
AAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAA  >  NZ_CP009273/4116972‑4117170

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