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 882 254.4 12957588 98.0% 12698436 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,906 A→G T259A (ACG→GCG)  ilvD → dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,9060AG97.1% 918.2 / 9.2 277T259A (ACG→GCG) ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (5/3);  new base G (133/136);  total (138/139)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 287 total aligned reads displayed.

GTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGT  >  NC_000913/3951816‑3951995
                                                                                           |                                                                                         
gttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCg                                                                                         >  1:6393974/1‑93 (MQ=255)
 ttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGgcgc                                                                                        >  1:6121965/1‑93 (MQ=255)
 ttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGgcgc                                                                                        >  2:5492293/1‑93 (MQ=255)
 ttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGgcgc                                                                                        >  1:2824833/1‑93 (MQ=255)
 ttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGc                                                                                        >  2:730812/1‑93 (MQ=255)
 ttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCCCCAGTAAGGCGGCGTTTGAAAACGCCATGgcgc                                                                                        <  2:5205755/93‑1 (MQ=255)
   gAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTg                                                                                      >  2:4113013/1‑93 (MQ=255)
   gAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTg                                                                                      <  1:3289705/93‑1 (MQ=255)
    aaTTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTgg                                                                                     >  1:5161428/1‑93 (MQ=255)
    aaTTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTgg                                                                                     >  1:4497913/1‑93 (MQ=255)
     aTTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGa                                                                                    <  1:4066028/93‑1 (MQ=255)
       tGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGata                                                                                  >  2:186574/1‑93 (MQ=255)
        gACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGatat                                                                                 >  2:5032885/1‑93 (MQ=255)
        gACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGatat                                                                                 <  2:6334048/93‑1 (MQ=255)
         aCCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATc                                                                                >  2:1587158/1‑93 (MQ=255)
           cAAACGTTATTACGAGCAAAGCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATcgc                                                                              <  2:3371355/93‑1 (MQ=255)
           cAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATcgc                                                                              >  1:4988106/1‑93 (MQ=255)
            aaaCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATcgcg                                                                             >  1:4697705/1‑93 (MQ=255)
             aaCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATcgcgt                                                                            >  2:3905769/1‑92 (MQ=255)
             aaCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGa                                                                            >  1:92293/1‑93 (MQ=255)
              aCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAGGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGAt                                                                           <  2:2331413/93‑1 (MQ=255)
              aCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGAt                                                                           <  1:3910807/93‑1 (MQ=255)
              aCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGAt                                                                           >  1:4320489/1‑93 (MQ=255)
              aCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGAt                                                                           <  2:1927326/93‑1 (MQ=255)
              aCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGAt                                                                           >  1:1977509/1‑93 (MQ=255)
              aCGTTAATACGAACAAAACGGAGAAAGGTCCCTTGCGGGTAATATCGCCAGTAAGGCGGCGTTTTGAAAACGCCATGGCGCTGGATATCGCGa                                                                            <  2:3133793/93‑1 (MQ=255)
                gTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATgg                                                                         <  2:2971030/93‑1 (MQ=255)
                gTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATgg                                                                         <  2:1764108/93‑1 (MQ=255)
                gTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATgg                                                                         >  1:1683969/1‑93 (MQ=255)
                 ttattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATggg                                                                        >  1:2635251/1‑93 (MQ=255)
                 ttattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATggg                                                                        <  2:4448155/93‑1 (MQ=255)
                 ttattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATggg                                                                        >  2:3105251/1‑93 (MQ=255)
                 ttattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATggg                                                                        <  1:3115332/93‑1 (MQ=255)
                  tattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGt                                                                       >  1:6201501/1‑93 (MQ=255)
                   attaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTg                                                                      >  2:3149250/1‑93 (MQ=255)
                   attaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTg                                                                      >  1:5809608/1‑93 (MQ=255)
                    ttaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTgg                                                                     <  2:5082198/93‑1 (MQ=255)
                     taCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGa                                                                    >  1:3053402/1‑93 (MQ=255)
                      aCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGAt                                                                   <  1:5905402/93‑1 (MQ=255)
                      aCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGAt                                                                   <  2:1288141/93‑1 (MQ=255)
                            aaaaCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTa                                                             >  1:2092391/1‑93 (MQ=255)
                             acacgacgAAAGGGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAACCGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTaa                                                            <  2:4616343/91‑1 (MQ=255)
                             aaaCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTaa                                                            >  1:5928986/1‑93 (MQ=255)
                               acgacgAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAaca                                                          >  2:5773940/1‑93 (MQ=255)
                               acgacgAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAaca                                                          >  2:3727337/1‑93 (MQ=255)
                                cgacgaAAGTGCTCTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAacac                                                         <  2:2539327/93‑1 (MQ=255)
                                cgacgaAAGTGCACTGCCACGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAacac                                                         <  2:2224388/93‑1 (MQ=255)
                                    gaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTa                                                     <  1:3604818/93‑1 (MQ=255)
                                    gaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTa                                                     >  2:3572289/1‑93 (MQ=255)
                                    gaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTa                                                     >  1:1059438/1‑93 (MQ=255)
                                     aaaGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTAc                                                    <  2:4501142/93‑1 (MQ=255)
                                        gTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTc                                                 <  1:5475307/93‑1 (MQ=255)
                                          gCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAc                                               <  2:4695607/93‑1 (MQ=255)
                                          gCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAc                                               <  2:2043666/93‑1 (MQ=255)
                                           cACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAcc                                              <  2:4601567/93‑1 (MQ=255)
                                           cACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAcc                                              <  2:1568815/93‑1 (MQ=255)
                                           cACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAcc                                              <  1:606802/93‑1 (MQ=255)
                                           cACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAcc                                              <  1:1432762/93‑1 (MQ=255)
                                            aCTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACct                                             <  2:6071861/93‑1 (MQ=255)
                                             cTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGCGGATCGACTAACACCGTACTTCACctg                                            >  1:4838289/1‑93 (MQ=255)
                                               gcagcgTAATATCGCCAGTAAGGCGGCGTTTGAACACGCCATTGCGCTTGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACctgct                                          >  2:6455836/4‑93 (MQ=255)
                                                ccGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACctgctg                                         <  1:293213/93‑1 (MQ=255)
                                                ccGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACctgctg                                         <  2:1421225/93‑1 (MQ=255)
                                                ccGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACctgctg                                         <  1:5625082/93‑1 (MQ=255)
                                                 cgcgTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTgg                                        >  1:679819/1‑93 (MQ=255)
                                                  gcgTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGCCTAACACCGTACTTCACCTGCTggc                                       <  1:4229670/93‑1 (MQ=255)
                                                   cgTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcg                                      >  1:4413366/1‑93 (MQ=255)
                                                      aaTATCGCCAGTGAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcggcg                                   <  1:383750/93‑1 (MQ=255)
                                                      aaTATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcggcg                                   <  2:2503567/93‑1 (MQ=255)
                                                      aaTATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcggcg                                   <  1:3723773/93‑1 (MQ=255)
                                                      aaTATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcggcg                                   <  1:5967662/93‑1 (MQ=255)
                                                       atatCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcggcgg                                  >  2:6158545/1‑93 (MQ=255)
                                                          tCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGgcgc                               >  2:148545/1‑93 (MQ=255)
                                                          tCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGgcgc                               <  1:3213618/93‑1 (MQ=255)
                                                           cGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCa                              <  2:3239032/93‑1 (MQ=255)
                                                            gCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCTCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAg                             >  1:3453612/1‑93 (MQ=255)
                                                            gCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAg                             <  1:6308391/93‑1 (MQ=255)
                                                              cAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGa                           >  2:2832026/1‑93 (MQ=255)
                                                               aGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCAACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGaa                          <  1:5034172/93‑1 (MQ=255)
                                                                gTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAg                         >  1:6458583/1‑93 (MQ=255)
                                                                 tAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGc                        <  1:4574684/93‑1 (MQ=255)
                                                                    ggcggcGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGa                     >  1:2480876/1‑93 (MQ=255)
                                                                     gcggcgTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGaa                    <  2:1156130/93‑1 (MQ=255)
                                                                       ggcgTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAAt                  <  1:1088913/93‑1 (MQ=255)
                                                                       ggcgTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAAt                  <  1:724176/93‑1 (MQ=255)
                                                                            tttGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGTCGCAGGAAGCGGAAATCGAc              >  2:2116548/1‑93 (MQ=255)
                                                                            tttGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGAc              >  2:5633913/1‑93 (MQ=255)
                                                                                aaaaCGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCa          >  1:244290/1‑93 (MQ=255)
                                                                                 aaaCGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCAc         >  1:69889/1‑93 (MQ=255)
                                                                                 aaaCGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCAc         >  1:479918/1‑93 (MQ=255)
                                                                                  aaCGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCAcc        <  1:1349103/93‑1 (MQ=255)
                                                                                  aaCGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGCACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCAcc        >  2:5894737/1‑93 (MQ=255)
                                                                                   aCGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCa       <  1:829606/93‑1 (MQ=255)
                                                                                   aCGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCa       >  1:4159009/1‑93 (MQ=255)
                                                                                    cGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCAt      <  2:1427174/93‑1 (MQ=255)
                                                                                     gCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATg     <  2:4293682/93‑1 (MQ=255)
                                                                                      ccATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa    <  2:3516029/93‑1 (MQ=255)
                                                                                      ccATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa    >  2:3013073/1‑93 (MQ=255)
                                                                                        aTGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGt  >  2:5397500/1‑93 (MQ=255)
                                                                                        aTGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGt  >  2:2667502/1‑93 (MQ=255)
                                                                                           |                                                                                         
GTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGT  >  NC_000913/3951816‑3951995

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGAT  >  NC_000913/3951813‑3951998
                                                                                             |                                                                                            
ATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGG                                                                                              >  SRR4302537.4897910/1‑94 (MQ=60)
ATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGG                                                                                              >  SRR4302537.3104399/1‑94 (MQ=60)
ATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGG                                                                                              >  SRR4302537.2008277/1‑94 (MQ=60)
ATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGG                                                                                              >  SRR4302537.3013619/1‑94 (MQ=60)
  TGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCG                                                                                            >  SRR4302537.5054821/1‑94 (MQ=60)
  TGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCG                                                                                            <  SRR4302537.115137/94‑1 (MQ=60)
  TGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCG                                                                                            >  SRR4302537.6437523/1‑94 (MQ=60)
  TGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCG                                                                                            >  SRR4302537.4101269/1‑94 (MQ=60)
  TGTAGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCG                                                                                            <  SRR4302537.2839017/94‑1 (MQ=60)
   GTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGC                                                                                           >  SRR4302537.6164044/1‑94 (MQ=60)
   GTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGC                                                                                           >  SRR4302537.2844578/1‑94 (MQ=60)
     TGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTG                                                                                         <  SRR4302537.3203049/94‑1 (MQ=60)
      GAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGG                                                                                        <  SRR4302537.1739200/94‑1 (MQ=60)
      GAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGG                                                                                        >  SRR4302537.5197011/1‑94 (MQ=60)
      GAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGG                                                                                        >  SRR4302537.4529025/1‑94 (MQ=60)
      GAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGG                                                                                        <  SRR4302537.6144191/94‑1 (MQ=60)
      GAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGG                                                                                        <  SRR4302537.3312954/94‑1 (MQ=60)
        ATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGAT                                                                                      <  SRR4302537.4094438/94‑1 (MQ=60)
        ATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGAT                                                                                      <  SRR4302537.5392404/94‑1 (MQ=60)
         TTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATA                                                                                     >  SRR4302537.4787685/1‑94 (MQ=60)
          TGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATAT                                                                                    <  SRR4302537.5550859/94‑1 (MQ=60)
          TGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATAT                                                                                    >  SRR4302537.4089087/1‑94 (MQ=60)
           GACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATC                                                                                   <  SRR4302537.68925/94‑1 (MQ=60)
            ACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCG                                                                                  <  SRR4302537.6374363/94‑1 (MQ=60)
            ACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTAGATATCG                                                                                  >  SRR4302537.2100643/1‑94 (MQ=60)
             CCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGC                                                                                 >  SRR4302537.5022478/1‑94 (MQ=60)
              CAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCG                                                                                >  SRR4302537.3413029/1‑94 (MQ=60)
              CAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCG                                                                                >  SRR4302537.4730204/1‑94 (MQ=60)
               AAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGA                                                                               >  SRR4302537.92961/1‑94 (MQ=60)
               AAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGA                                                                               <  SRR4302537.4977193/94‑1 (MQ=60)
               AAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGA                                                                               >  SRR4302537.6043153/1‑94 (MQ=60)
                AACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGAT                                                                              <  SRR4302537.4064853/94‑1 (MQ=60)
                AACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGAT                                                                              >  SRR4302537.1991251/1‑94 (MQ=60)
                AACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGAT                                                                              >  SRR4302537.4350523/1‑94 (MQ=60)
                 ACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATG                                                                             <  SRR4302537.4246933/94‑1 (MQ=60)
                 ACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATG                                                                             <  SRR4302537.3938144/94‑1 (MQ=60)
                  CGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGG                                                                            >  SRR4302537.1695759/1‑94 (MQ=60)
                  CGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGG                                                                            >  SRR4302537.6159140/1‑94 (MQ=60)
                  CGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGG                                                                            <  SRR4302537.6229639/94‑1 (MQ=60)
                  CGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGG                                                                            >  SRR4302537.4867181/1‑94 (MQ=60)
                   GTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGG                                                                           <  SRR4302537.4525513/94‑1 (MQ=60)
                   GTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGG                                                                           <  SRR4302537.6230816/94‑1 (MQ=60)
                   GTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGG                                                                           >  SRR4302537.2653693/1‑94 (MQ=60)
                    TTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGT                                                                          <  SRR4302537.4427769/94‑1 (MQ=60)
                    TTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGT                                                                          >  SRR4302537.6244239/1‑94 (MQ=60)
                    TTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGT                                                                          <  SRR4302537.3137089/94‑1 (MQ=60)
                     TATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTG                                                                         >  SRR4302537.5849440/1‑94 (MQ=60)
                     TATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTG                                                                         <  SRR4302537.6481942/94‑1 (MQ=60)
                      ATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGG                                                                        <  SRR4302537.956445/94‑1 (MQ=60)
                       TTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGA                                                                       >  SRR4302537.3074682/1‑94 (MQ=60)
                        TACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGAT                                                                      <  SRR4302537.1066895/94‑1 (MQ=60)
                         ACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATC                                                                     <  SRR4302537.5945925/94‑1 (MQ=60)
                          CGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCG                                                                    <  SRR4302537.2848603/94‑1 (MQ=60)
                             GCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACT                                                                 <  SRR4302537.6392158/94‑1 (MQ=60)
                              CAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTA                                                                >  SRR4302537.2107148/1‑94 (MQ=60)
                               AAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAGACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAA                                                               <  SRR4302537.4702360/94‑1 (MQ=60)
                               AAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAA                                                               >  SRR4302537.5969671/1‑94 (MQ=60)
                               AAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAA                                                               >  SRR4302537.5036353/1‑94 (MQ=60)
                                AAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAAC                                                              <  SRR4302537.5876777/94‑1 (MQ=60)
                                 AACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACA                                                             >  SRR4302537.3500467/1‑94 (MQ=60)
                                   CGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACC                                                           <  SRR4302537.4377319/94‑1 (MQ=60)
                                    GACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCG                                                          >  SRR4302537.3904327/1‑94 (MQ=60)
                                    GACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCG                                                          >  SRR4302537.4789331/1‑94 (MQ=60)
                                      CGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTA                                                        >  SRR4302537.1066741/1‑94 (MQ=60)
                                      CGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTA                                                        >  SRR4302537.318153/1‑94 (MQ=60)
                                       GAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTAC                                                       <  SRR4302537.4864445/94‑1 (MQ=60)
                                       GAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTAC                                                       <  SRR4302537.3630224/94‑1 (MQ=60)
                                           GTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCA                                                   <  SRR4302537.5512678/94‑1 (MQ=60)
                                           GTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCA                                                   <  SRR4302537.1650612/94‑1 (MQ=60)
                                            TGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAC                                                  >  SRR4302537.2005046/1‑94 (MQ=60)
                                              CACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCT                                                <  SRR4302537.610950/94‑1 (MQ=60)
                                              CACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCT                                                <  SRR4302537.1442731/94‑1 (MQ=60)
                                               ACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTG                                               >  SRR4302537.5117509/1‑94 (MQ=60)
                                               ACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTG                                               >  SRR4302537.6105988/1‑94 (MQ=60)
                                               ACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTG                                               <  SRR4302537.249556/94‑1 (MQ=60)
                                               ACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTG                                               >  SRR4302537.5999927/1‑94 (MQ=60)
                                               ACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTG                                               >  SRR4302537.3030467/1‑94 (MQ=60)
                                               ACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGCGGATCGACTAACACCGTACTTCACCTG                                               >  SRR4302537.4871689/1‑94 (MQ=60)
                                                CTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGC                                              >  SRR4302537.3858982/1‑94 (MQ=60)
                                                 TGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCT                                             <  SRR4302537.1290196/94‑1 (MQ=60)
                                                  GCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTG                                            <  SRR4302537.6259887/94‑1 (MQ=60)
                                                   CGGCGTAATATCGCCAGTAAGGCAGCGTTTGAAAACGCCAGGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGG                                           <  SRR4302537.5515309/94‑1 (MQ=60)
                                                   CCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGG                                           <  SRR4302537.292825/94‑1 (MQ=60)
                                                   CCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGG                                           <  SRR4302537.3557380/94‑1 (MQ=60)
                                                   CCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGG                                           <  SRR4302537.753045/94‑1 (MQ=60)
                                                   CCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGG                                           <  SRR4302537.295256/94‑1 (MQ=60)
                                                   CCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGG                                           >  SRR4302537.684480/1‑94 (MQ=60)
                                                   CCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGG                                           <  SRR4302537.5663548/94‑1 (MQ=60)
                                                    CGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGCTGGGTGGATCGACTAACACCGTACTTCACCTGCTGGG                                          >  SRR4302537.681578/1‑94 (MQ=60)
                                                    CGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGC                                          <  SRR4302537.51946/94‑1 (MQ=60)
                                                    CGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGC                                          >  SRR4302537.706522/1‑94 (MQ=60)
                                                     GCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGCCTAACACCGTACTTCACCTGCTGGCG                                         <  SRR4302537.4259118/94‑1 (MQ=60)
                                                     GCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCG                                         >  SRR4302537.186767/1‑94 (MQ=60)
                                                     GCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCG                                         <  SRR4302537.2007715/94‑1 (MQ=60)
                                                     GCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCG                                         >  SRR4302537.4443825/1‑94 (MQ=60)
                                                      CGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGG                                        >  SRR4302537.3842346/1‑94 (MQ=60)
                                                      CGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGG                                        >  SRR4302537.391646/1‑94 (MQ=60)
                                                       GTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGC                                       >  SRR4302537.5739971/1‑94 (MQ=60)
                                                       GTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGC                                       <  SRR4302537.4257103/94‑1 (MQ=60)
                                                       GTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACATGCTGGCGGG                                       >  SRR4302537.1270556/1‑94 (MQ=60)
                                                        TAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCG                                      >  SRR4302537.3257341/1‑94 (MQ=60)
                                                        TAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCG                                      >  SRR4302537.5867910/1‑94 (MQ=60)
                                                         AATATCGCCAGTGAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGG                                     <  SRR4302537.386408/94‑1 (MQ=60)
                                                         AATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGG                                     <  SRR4302537.2458736/94‑1 (MQ=60)
                                                         AATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGG                                     <  SRR4302537.6008620/94‑1 (MQ=60)
                                                         AATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGG                                     <  SRR4302537.3749939/94‑1 (MQ=60)
                                                          ATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGC                                    >  SRR4302537.3095363/1‑94 (MQ=60)
                                                           TATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCG                                   >  SRR4302537.1413754/1‑94 (MQ=60)
                                                           TATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCG                                   >  SRR4302537.94812/1‑94 (MQ=60)
                                                           TATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCG                                   <  SRR4302537.4169604/94‑1 (MQ=60)
                                                             TCGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCA                                 <  SRR4302537.3236233/94‑1 (MQ=60)
                                                              CGCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCTCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAG                                >  SRR4302537.3478035/1‑94 (MQ=60)
                                                               GCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTTGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGG                               <  SRR4302537.1413481/94‑1 (MQ=60)
                                                               GCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGG                               <  SRR4302537.5848923/94‑1 (MQ=60)
                                                               GCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGG                               <  SRR4302537.6351586/94‑1 (MQ=60)
                                                               GCCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACCAACACCGTACTTCACCTGCTGGCGGCGGCGCAGG                               <  SRR4302537.6351415/94‑1 (MQ=60)
                                                                CCAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGA                              >  SRR4302537.5959951/1‑94 (MQ=60)
                                                                 CAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAA                             <  SRR4302537.3296893/94‑1 (MQ=60)
                                                                 CAGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAA                             <  SRR4302537.4672600/94‑1 (MQ=60)
                                                                  AGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAG                            >  SRR4302537.6502352/1‑94 (MQ=60)
                                                                  AGTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCAACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAG                            <  SRR4302537.5068830/94‑1 (MQ=60)
                                                                   GTAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGC                           >  SRR4302537.6491589/1‑94 (MQ=60)
                                                                    TAAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCG                          <  SRR4302537.4606346/94‑1 (MQ=60)
                                                                     AAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGG                         <  SRR4302537.5328478/94‑1 (MQ=60)
                                                                     AAGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGG                         <  SRR4302537.1512055/94‑1 (MQ=60)
                                                                      AGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGA                        >  SRR4302537.5946816/1‑94 (MQ=60)
                                                                      AGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGA                        <  SRR4302537.727304/94‑1 (MQ=60)
                                                                      AGGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGA                        >  SRR4302537.2498314/1‑94 (MQ=60)
                                                                       GGCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAA                       <  SRR4302537.2703563/94‑1 (MQ=60)
                                                                        GCGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAA                      <  SRR4302537.4563490/94‑1 (MQ=60)
                                                                         CGGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAAT                     >  SRR4302537.3667676/1‑94 (MQ=60)
                                                                         CGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAAT                     <  SRR4302537.769769/94‑1 (MQ=60)
                                                                          GGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATC                    >  SRR4302537.3088340/1‑94 (MQ=60)
                                                                          GGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATC                    <  SRR4302537.1096436/94‑1 (MQ=60)
                                                                          GGCGTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATC                    <  SRR4302537.729118/94‑1 (MQ=60)
                                                                             GTTTGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGAC                 >  SRR4302537.279275/1‑94 (MQ=60)
                                                                                TGAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTC              >  SRR4302537.4562256/1‑94 (MQ=60)
                                                                                 GAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCA             >  SRR4302537.246021/1‑94 (MQ=60)
                                                                                 GAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCA             <  SRR4302537.4325169/94‑1 (MQ=60)
                                                                                 GAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCA             <  SRR4302537.6502979/94‑1 (MQ=60)
                                                                                 GAAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCA             <  SRR4302537.3738764/94‑1 (MQ=60)
                                                                                  AAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCAC            >  SRR4302537.70401/1‑94 (MQ=60)
                                                                                  AAAACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCAC            >  SRR4302537.483209/1‑94 (MQ=60)
                                                                                    AACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCA          <  SRR4302537.1358505/94‑1 (MQ=60)
                                                                                    AACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCA          >  SRR4302537.4188052/1‑94 (MQ=60)
                                                                                    AACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCA          <  SRR4302537.451801/94‑1 (MQ=60)
                                                                                     ACGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCAT         <  SRR4302537.835282/94‑1 (MQ=60)
                                                                                      CGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATG        >  SRR4302537.2170652/1‑94 (MQ=60)
                                                                                      CGCCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATG        <  SRR4302537.881955/94‑1 (MQ=60)
                                                                                        CCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAG      <  SRR4302537.2977303/94‑1 (MQ=60)
                                                                                        CCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAG      >  SRR4302537.2580430/1‑94 (MQ=60)
                                                                                        CCATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAG      <  SRR4302537.2515835/94‑1 (MQ=60)
                                                                                          ATGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTG    >  SRR4302537.1828915/1‑94 (MQ=60)
                                                                                          ATGGCGCTCGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTG    <  SRR4302537.89929/94‑1 (MQ=60)
                                                                                           TGGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGA   <  SRR4302537.4340460/94‑1 (MQ=60)
                                                                                            GGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGAT  <  SRR4302537.4304631/94‑1 (MQ=60)
                                                                                            GGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGAT  <  SRR4302537.225449/94‑1 (MQ=60)
                                                                                            GGCGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGAT  <  SRR4302537.3474726/94‑1 (MQ=60)
                                                                                             |                                                                                            
ATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGAT  >  NC_000913/3951813‑3951998

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