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 3,934,995 +C intergenic (‑308/‑173) yieP ← / → BW25113_RS19520 FadR/GntR family transcriptional regulator/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,934,9911.C96.7% 326.8 / 2.0 92intergenic (‑304/‑177)yieP/BW25113_RS19520FadR/GntR family transcriptional regulator/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (49/40);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGAT  >  NZ_CP009273/3934902‑3935075
                                                                                           |                                                                                    
ccgccgGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGcaca                                                                                       >  2:2910153/1‑90 (MQ=33)
     gggTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACAccccc                                                                                 >  1:6462/1‑89 (MQ=18)
          aGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGcgcc                                                                             >  2:981191/1‑90 (MQ=255)
          aGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGcgcc                                                                             >  1:598761/1‑90 (MQ=255)
             ggggTTCTCCTGAGAACTCCGGCAGAGAAAGCAGAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCt                                                                          <  1:646763/90‑1 (MQ=11)
             ggggTTCTCCTGAGAACTACGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTGGCGGGAAAGCGTATTATGCACACCCCCCCGCCGCt                                                                          <  1:2699355/90‑1 (MQ=11)
                gTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTgag                                                                       >  1:384049/1‑90 (MQ=255)
                gTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTgag                                                                       >  1:2451229/1‑90 (MQ=255)
                       tGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagc                                                                <  1:971154/90‑1 (MQ=255)
                       tGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagc                                                                <  1:2439398/90‑1 (MQ=255)
                        gagaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagca                                                               >  1:893594/1‑90 (MQ=255)
                        gagaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagca                                                               >  1:2375898/1‑90 (MQ=255)
                          gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  2:1910774/1‑90 (MQ=255)
                          gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  2:1406528/1‑90 (MQ=255)
                          gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  2:2516079/1‑90 (MQ=255)
                          gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  2:2777866/1‑90 (MQ=255)
                          gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  1:256821/1‑90 (MQ=255)
                          gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  2:2777867/1‑90 (MQ=255)
                          gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  2:1564347/1‑90 (MQ=255)
                          gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  1:2148922/1‑90 (MQ=255)
                          gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             <  1:2085889/90‑1 (MQ=255)
                          gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAaaagcaaag                                                            >  1:2757934/1‑90 (MQ=38)
                            aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaagc                                                           >  2:1848478/1‑90 (MQ=255)
                            aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaagc                                                           <  1:20884/90‑1 (MQ=255)
                             cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCg                                                          >  2:1305099/1‑90 (MQ=255)
                             cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCg                                                          >  2:651167/1‑90 (MQ=255)
                             cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCg                                                          >  1:874197/1‑90 (MQ=255)
                             cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCg                                                          >  1:1971813/1‑90 (MQ=255)
                               ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCt                                                                          >  2:1006564/1‑72 (MQ=255)
                               ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCt                                                                          <  1:1006564/72‑1 (MQ=255)
                               ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGc                                                        >  2:2347914/1‑90 (MQ=255)
                               ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGc                                                        >  1:639893/1‑90 (MQ=255)
                               ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGc                                                        >  2:88894/1‑90 (MQ=255)
                               ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGc                                                        >  1:1485406/1‑90 (MQ=255)
                                  gCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt                                                     <  1:634137/90‑1 (MQ=255)
                                  gCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt                                                     <  1:1848478/90‑1 (MQ=255)
                                  gCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt                                                     <  2:2018294/90‑1 (MQ=255)
                                    agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGc                                                   <  2:2668352/90‑1 (MQ=255)
                                    agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGc                                                   <  2:1633326/90‑1 (MQ=255)
                                     gagaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGct                                                  <  1:88894/90‑1 (MQ=255)
                                       gaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGctct                                                >  1:1877695/1‑90 (MQ=255)
                                       gaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGctct                                                >  1:698545/1‑90 (MQ=255)
                                         aaGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCttt                                              >  2:2134643/1‑90 (MQ=255)
                                            cAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAAc                                           <  2:3010825/90‑1 (MQ=255)
                                              aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACCCCCCCGCGCCGCTGAGATAAAGCAAAGCGGCACTGCTCTTTAACaa                                         <  2:2900866/90‑1 (MQ=11)
                                              aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACaa                                         <  2:384049/90‑1 (MQ=255)
                                              aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACaa                                         <  1:2140060/90‑1 (MQ=255)
                                              aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACaa                                         <  2:13687/90‑1 (MQ=255)
                                                 ataaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAAttt                                      >  2:306548/1‑90 (MQ=255)
                                                    aatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATc                                   <  1:2239401/90‑1 (MQ=255)
                                                      tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                 <  2:2148922/90‑1 (MQ=255)
                                                      tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                 <  2:1109729/90‑1 (MQ=255)
                                                      tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                 <  1:1910774/90‑1 (MQ=255)
                                                        cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                               <  1:2652006/90‑1 (MQ=255)
                                                        cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                               <  2:518963/90‑1 (MQ=255)
                                                          tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACaa                             <  2:639893/90‑1 (MQ=255)
                                                          tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACaa                             <  1:981191/90‑1 (MQ=255)
                                                           gACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAAt                            >  1:224965/1‑90 (MQ=255)
                                                            aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                           <  2:1476156/90‑1 (MQ=255)
                                                            aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                           <  2:2451229/90‑1 (MQ=255)
                                                            aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                           <  2:2554150/90‑1 (MQ=255)
                                                            aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                           <  2:519709/90‑1 (MQ=255)
                                                             ctctGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCt                          <  2:2365567/90‑1 (MQ=255)
                                                               ctGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                           >  1:2387514/1‑87 (MQ=255)
                                                               ctGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                           <  2:2387514/87‑1 (MQ=255)
                                                               ctGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgt                        >  1:1392454/1‑90 (MQ=255)
                                                                 gTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgtgt                      >  2:1945330/1‑90 (MQ=255)
                                                                 gTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgtgt                      >  2:2729482/1‑90 (MQ=255)
                                                                 gTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgtgt                      >  2:2342977/1‑90 (MQ=255)
                                                                     cGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGc                  <  2:1997021/90‑1 (MQ=255)
                                                                     cGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGc                  <  1:701280/90‑1 (MQ=255)
                                                                     cGGGAAAGCGTATTATGCACACCCCCCCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGc                  <  2:1971813/90‑1 (MQ=11)
                                                                      gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                 >  2:2996772/1‑90 (MQ=255)
                                                                      gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                 >  1:1042282/1‑90 (MQ=255)
                                                                      gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                 >  2:1216787/1‑90 (MQ=255)
                                                                      gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                 >  2:1800199/1‑90 (MQ=255)
                                                                      gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                 >  1:877386/1‑90 (MQ=255)
                                                                       ggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCAc                <  1:350693/90‑1 (MQ=255)
                                                                         aaaGCGTATTATGCACACCCCCGCGCCGCTGGGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc              >  2:2258396/1‑90 (MQ=11)
                                                                         aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc              >  1:448228/1‑90 (MQ=255)
                                                                         aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc              >  2:748143/1‑90 (MQ=255)
                                                                         aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc              >  1:1462736/1‑90 (MQ=255)
                                                                         aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc              >  2:2251388/1‑90 (MQ=255)
                                                                          aaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCg             <  2:1090929/90‑1 (MQ=255)
                                                                           aGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGa            >  2:2617655/1‑90 (MQ=255)
                                                                             cGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAg          <  1:25145/90‑1 (MQ=255)
                                                                              gTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTTGGCACTCGAAGa         >  1:1796324/1‑90 (MQ=11)
                                                                               tattatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGAt        <  1:2342977/90‑1 (MQ=255)
                                                                                attatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATa       >  1:609751/1‑90 (MQ=255)
                                                                                 ttatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATAc      <  2:1744530/90‑1 (MQ=255)
                                                                                 ttatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg     <  1:518099/90‑1 (MQ=33)
                                                                                  tatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg     >  2:1947203/1‑90 (MQ=255)
                                                                                  tatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg     <  1:306548/90‑1 (MQ=255)
                                                                                     gCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGAt  >  2:2268550/1‑90 (MQ=255)
                                                                                           |                                                                                    
CCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGAT  >  NZ_CP009273/3934902‑3935075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ 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
ACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAG  >  NZ_CP009273/3934897‑3935090
                                                                                                     |                                                                                               
ACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCC                                                                                                >  SRR3722030.6564/1‑100 (MQ=8)
    gtcagatgtgtataagagacagCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGC                                                                                             <  SRR3722030.1418985/78‑1 (MQ=9)
    cagatgtgtataagagacagCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGC                                                                                             <  SRR3722030.2671886/80‑1 (MQ=8)
    GCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGC                                                                                             <  SRR3722030.3037737/100‑1 (MQ=7)
    GCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGC                                                                                             <  SRR3722030.2392041/100‑1 (MQ=7)
    GCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGC                                                                                             <  SRR3722030.2538786/100‑1 (MQ=7)
     CCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCC                                                                                            >  SRR3722030.608033/1‑100 (MQ=7)
           GGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAG                                                                                      >  SRR3722030.389800/1‑100 (MQ=7)
           GGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAG                                                                                      >  SRR3722030.2495998/1‑100 (MQ=7)
              cggagtGGTTCTTCTGAGAACTCCGGCAGAGAAAGCAAAACTAAATGCTTGGCTCTGTAGCGGGAAAGCCAATTATGCAATCCCCCGCGCCGCTGAGAAA                                                                                   <  SRR3722030.818961/94‑1 (MQ=4)
                  agagacagCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGC                                                                               <  SRR3722030.1021791/92‑1 (MQ=7)
                  GGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAGAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGC                                                                               <  SRR3722030.656791/100‑1 (MQ=6)
                  GGGGTTCTCCTGAGAACTACGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTGGCGGGAAAGCGTATTATGCA‑‑CACCCCCCCGCCGCTGAGAAAAAGC                                                                               <  SRR3722030.2749197/100‑1 (MQ=6)
                   GGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCA                                                                              >  SRR3722030.2419103/1‑100 (MQ=9)
                   GGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCA                                                                              >  SRR3722030.907307/1‑100 (MQ=9)
                     GTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAA                                                                            >  SRR3722030.260528/1‑100 (MQ=19)
                     GTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAA                                                                            >  SRR3722030.2187361/1‑100 (MQ=19)
                     GTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCGCGCCGCTGAGAAAAAGCAAAG                                                                           >  SRR3722030.2808873/1‑100 (MQ=23)
                        CTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCG                                                                         >  SRR3722030.2006342/1‑100 (MQ=24)
                        CTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCG                                                                         >  SRR3722030.887608/1‑100 (MQ=24)
                          CCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGC                                                                       >  SRR3722030.649824/1‑100 (MQ=24)
                          CCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGC                                                                       >  SRR3722030.1509189/1‑100 (MQ=24)
                            TGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCAC                                                                     <  SRR3722030.985952/100‑1 (MQ=24)
                            TGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCAC                                                                     <  SRR3722030.2483897/100‑1 (MQ=24)
                               GAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGC                                                                  <  SRR3722030.2122985/100‑1 (MQ=24)
                                 gatctacactggccatgtcgtcggcagcgtcagatgtgtataagagacagGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTC                                                                <  SRR3722030.473617/50‑1 (MQ=10)
                                 gagatctacactggccatgtcgtcggcagcgtcagatgtgtataagagacagATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTC                                                                <  SRR3722030.2915501/48‑1 (MQ=20)
                                 ACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTC                                                                <  SRR3722030.21225/100‑1 (MQ=24)
                                  CTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCT                                                               >  SRR3722030.1909956/1‑100 (MQ=24)
                                  CTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCT                                                               >  SRR3722030.709349/1‑100 (MQ=24)
                                       GCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAAC                                                          <  SRR3722030.643983/100‑1 (MQ=24)
                                       GCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAAC                                                          <  SRR3722030.1879993/100‑1 (MQ=24)
                                          GAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAAT                                                       <  SRR3722030.90223/100‑1 (MQ=24)
                                                   AAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAC                                              <  SRR3722030.2178331/100‑1 (MQ=24)
                                                     cactggccatgtcgtcggcagcgtcagatgtgtataagagacagcCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAC                                              <  SRR3722030.2534160/55‑1 (MQ=11)
                                                      ATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAAT                                           >  SRR3722030.228177/1‑100 (MQ=24)
                                                         AATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTG                                        <  SRR3722030.2279807/100‑1 (MQ=24)
                                                          ATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGT                                       >  SRR3722030.2430943/1‑100 (MQ=24)
                                                          ATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGT                                       >  SRR3722030.1414340/1‑100 (MQ=24)
                                                          ATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCcgcctgtctcttatacacatctgacgctgccgacgacatggccagtgtagatctcggtg                                      >  SRR3722030.529260/1‑41 (MQ=24)
                                                           atgtcgtcggcagcgtcagatgtgtataagagacagcacCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTG                                        <  SRR3722030.945457/61‑1 (MQ=10)
                                                           TGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTG                                      <  SRR3722030.1943854/100‑1 (MQ=24)
                                                             CTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTG                                    <  SRR3722030.2700833/100‑1 (MQ=24)
                                                               gtcgtcggcagcgtcagatgtgtataagagacagcCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTG                                    <  SRR3722030.446583/65‑1 (MQ=10)
                                                               TGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGG                                  <  SRR3722030.996112/100‑1 (MQ=24)
                                                                 ACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCA                                >  SRR3722030.1057979/1‑100 (MQ=24)
                                                                 ACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCA                                >  SRR3722030.890844/1‑100 (MQ=24)
                                                                    CTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTC                             >  SRR3722030.455041/1‑100 (MQ=24)
                                                                    CTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTC                             >  SRR3722030.1486054/1‑100 (MQ=24)
                                                                         GCGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTTGGCACTCGAAGA                        >  SRR3722030.1826539/1‑100 (MQ=23)
                                                                          CGGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGAT                       <  SRR3722030.712120/100‑1 (MQ=24)
                                                                           GGGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATA                      >  SRR3722030.619204/1‑100 (MQ=24)
                                                                            tgtataagagacagatgcacacCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGAT                       <  SRR3722030.1453041/78‑1 (MQ=8)
                                                                            GGAAAGCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATAC                     <  SRR3722030.355889/100‑1 (MQ=24)
                                                                                 agGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATT                <  SRR3722030.2163197/98‑1 (MQ=6)
                                                                                 GCGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATT                >  SRR3722030.1275917/1‑100 (MQ=6)
                                                                                  agagacagatgcacacCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGAT                 <  SRR3722030.671734/84‑1 (MQ=8)
                                                                                  CGTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTC               <  SRR3722030.25546/100‑1 (MQ=6)
                                                                                   GTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGctgtctcttatacacatctgacgctgccgacgac              >  SRR3722030.2413322/1‑66 (MQ=8)
                                                                                   GTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCT              >  SRR3722030.1317728/1‑100 (MQ=6)
                                                                                   GTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCT              >  SRR3722030.2429758/1‑100 (MQ=6)
                                                                                   GTATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGGGGGCACTCGAAGATACGGATTCT              >  SRR3722030.2475412/1‑100 (MQ=5)
                                                                                    TATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTT             <  SRR3722030.2385560/100‑1 (MQ=6)
                                                                                     ATTATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACctgtctcttatacacatctgacgctgccgaccacatggccagtgta            >  SRR3722030.3038091/1‑54 (MQ=9)
                                                                                      TTATGCA‑‑CACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAAC          <  SRR3722030.526057/100‑1 (MQ=7)
                                                                                       TATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAAC          <  SRR3722030.311012/100‑1 (MQ=6)
                                                                                        ATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACG         >  SRR3722030.2147837/1‑100 (MQ=6)
                                                                                        ATGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACG         >  SRR3722030.315556/1‑100 (MQ=6)
                                                                                         TGCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGT        <  SRR3722030.305771/100‑1 (MQ=6)
                                                                                          GCA‑‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTC       >  SRR3722030.2910738/1‑100 (MQ=6)
                                                                                           CA‑‑CCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCA    >  SRR3722030.2881236/1‑100 (MQ=6)
                                                                                                CCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAG  >  SRR3722030.2046521/1‑100 (MQ=7)
                                                                                                CCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAG  >  SRR3722030.2937868/1‑100 (MQ=7)
                                                                                                CCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAG  >  SRR3722030.1647576/1‑100 (MQ=7)
                                                                                                CCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAG  >  SRR3722030.1225057/1‑100 (MQ=7)
                                                                                                     |                                                                                               
ACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCA‑‑CACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAG  >  NZ_CP009273/3934897‑3935090

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