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 F1 I2 R1 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,349,220 (G)7→9 intergenic (‑157/+211) fabI ← / ← ycjD enoyl‑[acyl‑carrier‑protein] reductase (NADH)/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,349,2131.G90.4% 226.1 / 20.1 83intergenic (‑150/+218)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase (NADH)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/6);  new base G (17/58);  total (19/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01
*NC_0009131,349,2132.G81.9% 203.8 / 46.9 83intergenic (‑150/+218)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase (NADH)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (5/10);  new base G (14/54);  total (19/64)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 115 total aligned reads displayed.

TAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCA  >  NC_000913/1349122‑1349292
                                                                                                 ||                                                                                        
tAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g                                                                                         <  1:6056945/93‑1 (MQ=255)
 aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑g                                                                                      >  1:2415666/1‑93 (MQ=255)
 aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑g                                                                                      >  2:3387342/1‑93 (MQ=255)
 aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑g                                                                                      >  2:5388314/1‑93 (MQ=255)
   tAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg                                                                                    <  1:4704316/93‑1 (MQ=255)
   tAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg                                                                                    <  1:1221315/93‑1 (MQ=255)
   tAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg                                                                                    >  1:783208/1‑93 (MQ=255)
    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑g                                                                                 <  1:4597534/93‑1 (MQ=255)
    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑g                                                                                 <  2:3373784/93‑1 (MQ=255)
    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑g                                                                                 >  1:1006415/1‑93 (MQ=255)
     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gg                                                                                <  2:2121958/93‑1 (MQ=255)
     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gg                                                                                <  2:1236626/93‑1 (MQ=255)
     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gg                                                                                >  2:5626174/1‑93 (MQ=255)
     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gg                                                                                >  2:5998091/1‑93 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               <  1:3449201/93‑1 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               <  1:616245/93‑1 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               <  1:1135277/93‑1 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               <  2:5545302/93‑1 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               >  2:736367/1‑93 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCGGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               >  1:1177168/1‑93 (MQ=255)
       cAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gggg                                                                              >  1:4029118/1‑92 (MQ=255)
        aGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggggg                                                                             >  1:4044465/1‑91 (MQ=255)
        aGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggggg                                                                             >  1:5832396/1‑91 (MQ=255)
         gTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gggggg                                                                            >  2:5249082/1‑90 (MQ=255)
         gTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGg‑‑ggc                                                                          <  2:2429656/93‑1 (MQ=255)
         gTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGg‑‑ggc                                                                          >  2:642843/1‑93 (MQ=255)
           aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGg‑‑ggcgg                                                                        >  1:29413/1‑91 (MQ=255)
             gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGg‑‑ggcgctt                                                                      >  2:237707/1‑89 (MQ=255)
             gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGg‑‑ggcgctt                                                                      >  2:4956103/1‑89 (MQ=255)
              aaCGGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGG‑‑GGG‑‑GGG‑‑CGCtttt                                                                     >  1:942791/1‑93 (MQ=255)
                aaGATAAACGGTTATTATAATCACCCTGGCTGTCAGTAGTTATAGTTGCCAGGTCGGACCGGAGCAGGATGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTc                                                                    <  2:2853589/92‑1 (MQ=255)
                  gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTccc                                                                  >  1:4095644/1‑93 (MQ=255)
                   aTAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTtccc                                                                 >  1:979596/1‑93 (MQ=255)
                    tAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCt                                                                <  1:1352558/93‑1 (MQ=255)
                     aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCTc                                                               >  2:771180/1‑93 (MQ=255)
                      aaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTccccccc                                                              >  1:2040629/1‑90 (MQ=255)
                       aCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCTCAc                                                             >  2:1813759/1‑93 (MQ=255)
                        cGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCCCAcc                                                            >  2:756075/1‑93 (MQ=255)
                         ggTTATTATAGTCACCCTGGCTGTGAGTAGTTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCAccc                                                           <  1:4288771/93‑1 (MQ=255)
                         ggTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCAccc                                                           <  2:5645069/93‑1 (MQ=255)
                         ggTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCt                                                          >  1:4233108/1‑93 (MQ=255)
                          gTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCccccccccc                                                          >  1:5659498/1‑92 (MQ=255)
                           tttttaaAATCAACCTGGGCTGTAATAGCTAAAGTTCCCCGGTCCGGCCGGGACCGGCCTGGG‑‑GAGG‑‑G‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTa                                                         <  2:3362634/90‑1 (MQ=255)
                           ttattaTAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTa                                                         <  1:2971150/93‑1 (MQ=255)
                            tattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCccccccccccc                                                        >  1:1538658/1‑90 (MQ=255)
                             attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAAc                                                     <  1:4472311/93‑1 (MQ=255)
                             attatAATCAACCTGACTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAAc                                                     <  2:5962353/93‑1 (MQ=255)
                               tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCct                                                  >  1:3177827/1‑93 (MQ=255)
                               tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCct                                                  <  2:3024143/93‑1 (MQ=255)
                                 taATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCctc                                                 <  2:3142863/93‑1 (MQ=255)
                                   aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCtctc                                               <  1:4562490/93‑1 (MQ=39)
                                    tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTcc                                              <  2:4425908/93‑1 (MQ=39)
                                     cAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTccc                                             >  1:2143565/1‑93 (MQ=39)
                                         cTGGCTGTGATTAGCTATAGTTCCCGGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCaga                                         <  2:2759138/93‑1 (MQ=37)
                                         cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTTCCCCCCCCCCACCCCCTccccccc                                         >  1:3820194/1‑90 (MQ=255)
                                         cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCCCCCCCCAAACCCCTCCCCcga                                         >  2:3030242/1‑93 (MQ=255)
                                         cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCaga                                         <  2:2504697/93‑1 (MQ=38)
                                            ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgg                                       <  1:3120308/93‑1 (MQ=38)
                                            ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgg                                       <  1:5911470/93‑1 (MQ=38)
                                             gCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAggg                                      <  1:3390771/93‑1 (MQ=38)
                                              cTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgggg                                     <  1:778632/93‑1 (MQ=37)
                                              cTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgggg                                     <  2:3580948/93‑1 (MQ=37)
                                              cTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgggg                                     <  1:5869653/93‑1 (MQ=37)
                                               tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCCCCCCCCTAACCCTCCTCCCACCgggg                                     >  2:5825638/1‑93 (MQ=255)
                                               tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCCCACCCCCACCCCCTCCCCCCAggggg                                    >  2:2987704/1‑92 (MQ=255)
                                                gtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCg                                   <  1:4909907/93‑1 (MQ=37)
                                                 tgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                  <  2:5278214/93‑1 (MQ=37)
                                                 tgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                  <  1:2081604/93‑1 (MQ=37)
                                                 tgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                  <  2:3472245/93‑1 (MQ=37)
                                                  gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                 <  1:6144376/93‑1 (MQ=37)
                                                  gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                 <  1:1256274/93‑1 (MQ=37)
                                                  gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                 <  1:1196156/93‑1 (MQ=37)
                                                  gAGTAGCTATAGTTGACAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                 <  2:2944012/93‑1 (MQ=25)
                                                   aGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgg                                <  1:4738844/93‑1 (MQ=37)
                                                   aGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgg                                <  1:1943569/93‑1 (MQ=37)
                                                    gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCGGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAggg                               <  2:3151092/93‑1 (MQ=25)
                                                     tggttataTTTGCCAGTCCCGACGGGGGCAGCCTGCGCCGGG‑‑G‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgggg                              <  2:2480520/89‑1 (MQ=9)
                                                     tAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgggg                              <  2:5475066/93‑1 (MQ=25)
                                                      aGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGa                             <  2:4054790/93‑1 (MQ=25)
                                                       gCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCCCCCCCCCAACCCCTCCCCCCCGgggggggggggg                              >  1:5031659/1‑93 (MQ=255)
                                                       gCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGAc                            <  2:768334/93‑1 (MQ=25)
                                                        gtataGTTGCCAGGGCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGAcc                           <  2:5496529/92‑1 (MQ=21)
                                                        cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGAcc                           <  1:1589656/93‑1 (MQ=25)
                                                          ataGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCCCCCCCCCACCCCCCCCCCCCAGAGggggggggggcgc                           >  1:2789774/1‑90 (MQ=255)
                                                            aGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAt                       <  1:2362574/93‑1 (MQ=25)
                                                             gTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAtt                      <  2:818272/93‑1 (MQ=25)
                                                              ttGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtg                     <  1:3725250/93‑1 (MQ=25)
                                                               tGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtgt                    <  2:326072/93‑1 (MQ=25)
                                                                      tCCGACCGGAGCAGGCTGCG‑‑GCAGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAatat            <  2:3909769/93‑1 (MQ=255)
                                                                      cgcGACCGGAGCAGGGTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCGCCCTAACCCCGCCCCCAGAGGGGCGAGGGGACCTGACTGGGCAAata             <  2:5790339/91‑1 (MQ=16)
                                                                       ccGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAatat            <  1:334463/93‑1 (MQ=25)
                                                                       ccGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAatat            <  1:1620376/93‑1 (MQ=25)
                                                                         caCCGGAGCAGTCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTg          <  2:6308122/92‑1 (MQ=21)
                                                                           ccGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGGGTGGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGtt        <  1:342607/93‑1 (MQ=25)
                                                                               gAGCAGGCTGCG‑‑GCA‑‑GGAGGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAcc     <  2:3517896/93‑1 (MQ=25)
                                                                                 gCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAcccc   <  2:84701/93‑1 (MQ=25)
                                                                                  cAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTATCCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa  <  2:3026622/93‑1 (MQ=25)
                                                                                  cAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa  <  2:5783191/93‑1 (MQ=25)
                                                                                  cAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa  <  1:5772779/93‑1 (MQ=25)
                                                                                  cAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGGGCCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa  <  2:5436536/93‑1 (MQ=25)
                                                                                                 ||                                                                                        
TAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCA  >  NC_000913/1349122‑1349292

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAG  >  NC_000913/1349127‑1349313
                                                                                                   |                                                                                             
AACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCGGCGGC‑‑A‑‑GGGG‑‑GGG                                                                                               >  SRR4302533.1184233/1‑94 (MQ=60)
 ACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGG                                                                                              <  SRR4302533.1142028/94‑1 (MQ=60)
 ACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGG                                                                                              <  SRR4302533.3470496/94‑1 (MQ=60)
 ACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGG                                                                                              >  SRR4302533.4053640/1‑94 (MQ=60)
 ACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGG                                                                                              <  SRR4302533.619864/94‑1 (MQ=60)
  CAGTACGAGCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGGG                                                                                             >  SRR4302533.754473/1‑94 (MQ=60)
  CAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGGG                                                                                             >  SRR4302533.5867645/1‑94 (MQ=60)
  CAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGGG                                                                                             >  SRR4302533.4069082/1‑94 (MQ=60)
     TACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggggc                                                                                          >  SRR4302533.2682121/1‑89 (MQ=60)
     TACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgg                                                                                          >  SRR4302533.29612/1‑89 (MQ=60)
        GAACGGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGG‑‑GGGCGCTTTT                                                                                        >  SRR4302533.948337/1‑94 (MQ=60)
        GAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTT                                                                                         >  SRR4302533.5259953/1‑94 (MQ=60)
        GAACAGATAAACGGTCATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTT                                                                                         >  SRR4302533.4004132/1‑94 (MQ=60)
         AACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTT                                                                                        <  SRR4302533.4374813/94‑1 (MQ=60)
          ACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTC                                                                                       <  SRR4302533.1538041/94‑1 (MQ=60)
          ACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTC                                                                                       <  SRR4302533.1304795/94‑1 (MQ=60)
           CAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCCTTTCC                                                                                      >  SRR4302533.3957815/1‑94 (MQ=60)
            AGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCC                                                                                     >  SRR4302533.4120578/1‑94 (MQ=60)
             GATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTTCCC                                                                                    >  SRR4302533.985339/1‑94 (MQ=60)
               TAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTC                                                                                  <  SRR4302533.1360692/94‑1 (MQ=60)
               TAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgcttttttcgcc                                                                                >  SRR4302533.2233314/1‑79 (MQ=60)
                AAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCccc                                                                                 >  SRR4302533.2053207/1‑91 (MQ=60)
                 AACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCGGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgcttttttcccccc                                                                              >  SRR4302533.3632991/1‑77 (MQ=60)
                 AACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCAC                                                                                <  SRR4302533.5245325/94‑1 (MQ=60)
                   CGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGG‑‑GGGCGCTTTTCCCCTCACCCT                                                                             >  SRR4302533.4258689/1‑94 (MQ=60)
                   CGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgctttttttccgcccc                                                                            >  SRR4302533.5399259/1‑75 (MQ=60)
                    GGTTATTATAGTCACCCTGGCTGTGAGTAGTTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCT                                                                             <  SRR4302533.4314561/94‑1 (MQ=60)
                    GGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgctttttcccccccccc                                                                           >  SRR4302533.5693699/1‑74 (MQ=60)
                      TTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAA                                                                           <  SRR4302533.2989320/94‑1 (MQ=60)
                      TTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgctttttttcctcctcccc                                                                         >  SRR4302533.3292247/1‑72 (MQ=60)
                      TTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgctttttcccccccccccc                                                                         >  SRR4302533.1547923/1‑72 (MQ=60)
                        ATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACC                                                                         <  SRR4302533.4499363/94‑1 (MQ=60)
                         TTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCT                                                                       >  SRR4302533.3197372/1‑94 (MQ=60)
                              ATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCC                                                                   <  SRR4302533.4590159/94‑1 (MQ=60)
                               TCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCC                                                                  >  SRR4302533.2156772/1‑94 (MQ=60)
                                   CCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgcttttttccccccccccaccccctccccccc                                                            >  SRR4302533.3843554/1‑59 (MQ=60)
                                    CTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGgcgctttttcccgccccccccccccccccccccccc                                                           >  SRR4302533.2357250/1‑58 (MQ=60)
                                     TGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGG                                                            <  SRR4302533.4600772/94‑1 (MQ=60)
                                     TGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGG                                                            <  SRR4302533.3161181/94‑1 (MQ=60)
                                      GGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGG                                                           <  SRR4302533.3139497/94‑1 (MQ=60)
                                      GGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGG                                                           <  SRR4302533.5947185/94‑1 (MQ=60)
                                      GGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgctttttttctccccctcccacacgcccccccccc                                                         >  SRR4302533.5351689/1‑56 (MQ=60)
                                      GGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgctttttcccccgcccccccacccccccccccccg                                                         >  SRR4302533.452215/1‑56 (MQ=60)
                                       GCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGG                                                          <  SRR4302533.3411711/94‑1 (MQ=60)
                                        CTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGC                                                         <  SRR4302533.5905131/94‑1 (MQ=60)
                                        CTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGC                                                         <  SRR4302533.783215/94‑1 (MQ=60)
                                          GTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGA                                                       <  SRR4302533.4939650/94‑1 (MQ=60)
                                           TGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAG                                                      <  SRR4302533.2094414/94‑1 (MQ=60)
                                           TGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGG                                                     <  SRR4302533.2209829/94‑1 (MQ=60)
                                            GAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGG                                                     <  SRR4302533.6181372/94‑1 (MQ=60)
                                            GAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGG                                                     <  SRR4302533.1263784/94‑1 (MQ=60)
                                            GAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGG                                                     <  SRR4302533.1203344/94‑1 (MQ=60)
                                             AGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGG                                                    <  SRR4302533.1955426/94‑1 (MQ=60)
                                             AGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGG                                                    <  SRR4302533.1428717/94‑1 (MQ=60)
                                             AGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGG                                                    <  SRR4302533.4767577/94‑1 (MQ=60)
                                               TAGCTATAGTTGCCAGGTCCGACCGGACCAGGCTGCGGCCGG‑‑GGGG‑‑GGGCGCTTTTCCCCctgcccccgcccctcacccgagccgggtcggggg                                                  >  SRR4302533.2039659/1‑60 (MQ=52)
                                                AGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgcttttttgcccccccccccccccccccccccccccccggggggg                                               >  SRR4302533.2215940/1‑46 (MQ=60)
                                                AGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgctttttcccccccccccaacccctcccccccggggggggggggg                                               >  SRR4302533.5062124/1‑46 (MQ=60)
                                                  CTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCG                                               <  SRR4302533.1599270/94‑1 (MQ=60)
                                                   TATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgcttttttcggccccgcaccctcccccccccccccccccgcgccgggc                                            >  SRR4302533.4314040/1‑43 (MQ=60)
                                                   TATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgctttttcccccccccccaccccccccccccagagggggggggggcgc                                            >  SRR4302533.2806930/1‑43 (MQ=60)
                                                      AGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATT                                           <  SRR4302533.2377193/94‑1 (MQ=60)
                                                      AGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgctttttgccccccccccccactaccccccccggcgcggggcggcggcggt                                         >  SRR4302533.2544911/1‑40 (MQ=60)
                                                       GTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgctttttcccccccccccccaccccccccccccgggggggggggggcggggt                                        >  SRR4302533.1493646/1‑39 (MQ=60)
                                                        TTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGT                                         <  SRR4302533.3748104/94‑1 (MQ=60)
                                                         TGCCAGGTCCGACTGAAGCAGGCTGCGGG‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTG                                        <  SRR4302533.3962750/94‑1 (MQ=49)
                                                          GCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgcttctttctcccccccccaccctccccccacgccgcgaggggaccgggcgtccg                                     >  SRR4302533.1111681/1‑36 (MQ=60)
                                                           CCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgcttttttgccccccccccacccactcccccccctgcggcggggggagacgttgta                                    >  SRR4302533.4791703/1‑35 (MQ=60)
                                                               GTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGggcgcttttcaccccccccccccccccccccccccccaaggaggggggagggagttttgaaca                                >  SRR4302533.5640030/1‑31 (MQ=60)
                                                                 CCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATT                                <  SRR4302533.1630192/94‑1 (MQ=60)
                                                                 CCGACCGGAGCAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATT                                <  SRR4302533.336507/94‑1 (MQ=60)
                                                                     CCGGAGCAGGCTGCGGC‑‑A‑‑GGGGGTGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTA                            <  SRR4302533.344699/94‑1 (MQ=60)
                                                                           CAGGCTGCGGC‑‑AGGGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAG                      <  SRR4302533.5807640/94‑1 (MQ=60)
                                                                                       caggGGGG‑‑GGGCGCTTTTCCCCTCTCCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAG            <  SRR4302533.519249/90‑1 (MQ=60)
                                                                                        G‑‑GGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGTGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCT         <  SRR4302533.4465857/94‑1 (MQ=60)
                                                                                         ggGGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGC          <  SRR4302533.1973002/92‑1 (MQ=60)
                                                                                              G‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATA     <  SRR4302533.753449/94‑1 (MQ=60)
                                                                                                  GGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACA   <  SRR4302533.4800219/94‑1 (MQ=60)
                                                                                                   GCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAG  <  SRR4302533.5950301/94‑1 (MQ=60)
                                                                                                   GCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAG  <  SRR4302533.2799924/94‑1 (MQ=60)
                                                                                                   GCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAG  <  SRR4302533.4712690/94‑1 (MQ=60)
                                                                                                   |                                                                                             
AACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGC‑‑A‑‑GGGG‑‑GGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAG  >  NC_000913/1349127‑1349313

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