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 2,579,447 T→C S598P (TCG→CCG)  tktB → transketolase 2 isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,579,4470TC99.5% 596.9 / NA 190S598P (TCG→CCG) tktBtransketolase 2 isozyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (91/98);  minor base A (1/0);  total (92/98)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGC  >  NC_000913/2579357‑2579535
                                                                                          |                                                                                        
gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCg                                                                                        <  1:121682/93‑1 (MQ=255)
  gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                      <  2:3919188/93‑1 (MQ=255)
  gaAATGAGCAGGAGAAGGTCGCAATGTACGCTTAGTTTCCCTGCCCTCGACCGATATTTTCGCCGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                      <  2:3182913/93‑1 (MQ=255)
  gaAATGAGCAGGAGAAGGTCGAAATGTACGCGTAGTTTCCCTGCCCGCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                      <  2:4535115/93‑1 (MQ=255)
   aaaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtg                                                                                     <  1:5511442/93‑1 (MQ=255)
    aaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtgt                                                                                    >  1:3550371/1‑93 (MQ=255)
    aaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgggt                                                                                    >  1:1828451/1‑93 (MQ=255)
      ttAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTg                                                                                  <  1:5322626/93‑1 (MQ=255)
        aGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGcc                                                                                >  1:6189712/1‑93 (MQ=255)
           aGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTc                                                                             <  1:5862962/93‑1 (MQ=255)
            ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCt                                                                            >  2:2944096/1‑93 (MQ=255)
            ggAGAAGGTCGCAATGTACGCGTAGTTCCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCt                                                                            >  2:5812837/1‑93 (MQ=255)
               gaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAAc                                                                         <  2:4834033/93‑1 (MQ=255)
                aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACg                                                                        <  1:1028829/93‑1 (MQ=255)
                 aGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGt                                                                       >  2:1425625/1‑93 (MQ=255)
                  ggTCGGAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGtt                                                                      <  1:4675992/93‑1 (MQ=255)
                  ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGtt                                                                      >  1:22779/1‑93 (MQ=255)
                   gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTg                                                                     <  2:3868273/93‑1 (MQ=255)
                     cGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCg                                                                   >  1:974620/1‑93 (MQ=255)
                      gCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCgg                                                                  >  1:251965/1‑93 (MQ=255)
                      gCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCgg                                                                  <  2:5063386/93‑1 (MQ=255)
                       cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                 >  2:4153323/1‑93 (MQ=255)
                       cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                 >  2:3544267/1‑93 (MQ=255)
                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                <  1:6248700/93‑1 (MQ=255)
                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                <  2:2838188/93‑1 (MQ=255)
                         aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTc                                                               <  2:5187141/93‑1 (MQ=255)
                            tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgcg                                                            <  1:4115909/93‑1 (MQ=255)
                            tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgcg                                                            >  1:3702149/1‑93 (MQ=255)
                            tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgcg                                                            <  1:1386410/93‑1 (MQ=255)
                              cgcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTg                                                          >  1:5341341/1‑93 (MQ=255)
                              cgcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTg                                                          <  2:2181411/93‑1 (MQ=255)
                               gcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTgg                                                         >  1:5190532/1‑93 (MQ=255)
                                cgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGc                                                        <  2:4738071/93‑1 (MQ=255)
                                cgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGc                                                        <  1:394255/93‑1 (MQ=255)
                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTgg                                                   >  2:586533/1‑93 (MQ=255)
                                      tccccTGCCCTCGACCGATATTTTCTACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGa                                                  <  2:4269616/91‑1 (MQ=255)
                                       tCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGaa                                                 <  1:2655854/93‑1 (MQ=255)
                                       tCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGaa                                                 <  2:5492644/93‑1 (MQ=255)
                                        cccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAg                                                <  1:4611796/93‑1 (MQ=255)
                                         ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAgg                                               >  2:1872371/1‑92 (MQ=255)
                                         ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGc                                               >  2:208148/1‑93 (MQ=255)
                                           tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGGGGGAGCAg                                             >  2:2551907/1‑93 (MQ=255)
                                            gCCCTCGACCGATATTTTCGACGGCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAgg                                            <  2:4689862/93‑1 (MQ=255)
                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCGTGGCGGGGGGAGAAgggtt                                         >  1:1617662/1‑90 (MQ=255)
                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAtt                                        <  2:5398759/93‑1 (MQ=255)
                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAtt                                        <  1:1117751/93‑1 (MQ=255)
                                                 cGACTGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTg                                       <  1:4507465/93‑1 (MQ=255)
                                                  gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGc                                      <  1:5963482/93‑1 (MQ=255)
                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGcc                                     <  1:2144224/93‑1 (MQ=255)
                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGTGGGAGCAGGGATTGcc                                     >  2:3347501/1‑93 (MQ=255)
                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGACTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAACAGGTATTGcg                                     >  1:4254773/1‑92 (MQ=255)
                                                    ccGATTTTTTCTACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCg                                    <  2:4474243/93‑1 (MQ=255)
                                                    ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCg                                    <  1:3935379/93‑1 (MQ=255)
                                                    ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGCCg                                    >  2:129860/1‑93 (MQ=255)
                                                    ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCTCGGTGGTCGACGGGGGGATTGCCg                                    >  2:5587859/1‑93 (MQ=255)
                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGa                                   <  2:144605/93‑1 (MQ=255)
                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGa                                   >  1:4638003/1‑93 (MQ=255)
                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGCCGa                                   >  2:3761214/1‑93 (MQ=255)
                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAGGTATTGCCGa                                   >  2:752535/1‑93 (MQ=255)
                                                       atatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGAtt                                 <  1:3082347/93‑1 (MQ=255)
                                                        tatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTa                                <  2:3039460/93‑1 (MQ=255)
                                                        tatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTa                                <  1:1472778/93‑1 (MQ=255)
                                                         atTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGGGGGGGGAGGAGGGGTTGCCGGtttt                               >  2:4502405/1‑91 (MQ=255)
                                                           tttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTg                             >  2:2747294/1‑93 (MQ=255)
                                                            ttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTgg                            <  2:4804074/93‑1 (MQ=255)
                                                             tCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGt                           <  1:5142668/93‑1 (MQ=255)
                                                             tCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGt                           <  2:2288254/93‑1 (MQ=255)
                                                               gACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTAc                         <  2:4973119/93‑1 (MQ=255)
                                                                  gCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCCTGGCGGTGGAAGCAGGTGTTGCCGCTTACTGGTACAAg                      >  1:4817379/1‑93 (MQ=255)
                                                                     cAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGTACAAGTAt                   >  2:17017/1‑93 (MQ=255)
                                                                     cAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGTACAAGTAt                   >  1:4999268/1‑93 (MQ=255)
                                                                      aGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGGACAAATATg                  >  2:2659400/1‑93 (MQ=255)
                                                                      aGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATg                  >  1:5889701/1‑93 (MQ=255)
                                                                       ggATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                 <  2:1637798/93‑1 (MQ=255)
                                                                       ggATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGt                 >  1:3037503/1‑93 (MQ=255)
                                                                       ggATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGt                 >  1:5691698/1‑93 (MQ=255)
                                                                        gATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGtt                <  2:2818710/93‑1 (MQ=255)
                                                                        gATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGCGGCGGTGGGAGCAGGGATTGCCGCTTACTGGGCCCAGTATGtt                >  2:2497370/1‑93 (MQ=255)
                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg               <  2:5356651/93‑1 (MQ=255)
                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAAGTATTGCCGATTACTGGGACAAGTATGGTg               >  2:6286539/1‑93 (MQ=255)
                                                                          tGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGGACAAGTATGTTgg              >  2:361841/1‑93 (MQ=255)
                                                                           gAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGCATTGCCGGTTACTGGTACCAGGATGTTTGt             >  2:4635415/1‑93 (MQ=255)
                                                                            aGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc            <  1:4159880/93‑1 (MQ=255)
                                                                             ggAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCt           <  1:5163025/93‑1 (MQ=255)
                                                                             ggAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCt           <  1:473285/93‑1 (MQ=255)
                                                                              gAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGGACACGTATGTTGGTCTg          >  2:527375/1‑93 (MQ=255)
                                                                              gAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTg          <  1:3851325/93‑1 (MQ=255)
                                                                               aaTATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGTTGGTCTGa         >  1:2311815/1‑93 (MQ=255)
                                                                                atatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaa        >  2:1931679/1‑93 (MQ=255)
                                                                                atatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaa        <  1:5443535/93‑1 (MQ=255)
                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa       <  2:243749/93‑1 (MQ=255)
                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa       <  2:4252942/93‑1 (MQ=255)
                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa       >  2:3475278/1‑93 (MQ=255)
                                                                                  atCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg      <  1:1031944/93‑1 (MQ=255)
                                                                                  atCGGGAGCCGGTGTGGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg      <  1:1066555/93‑1 (MQ=255)
                                                                                   tCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgg     <  1:6057548/93‑1 (MQ=255)
                                                                                   tCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGTATTGCCGATTACTTGTACAAGTATGTTGGTCTGAAAgg     >  2:3277211/1‑93 (MQ=38)
                                                                                      ggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGc  >  2:3539704/1‑93 (MQ=255)
                                                                                      ggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGc  >  1:6094806/1‑93 (MQ=255)
                                                                                      ggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGACGGTGGAAGCATGTATTGCCGATTACTGGTACAAGTATGTAGGTCTGAAAGGGGc  <  1:3195938/93‑1 (MQ=38)
                                                                                          |                                                                                        
GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGC  >  NC_000913/2579357‑2579535

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ 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
CAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTG  >  NC_000913/2579355‑2579540
                                                                                            |                                                                                             
CAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCC                                                                                              >  SRR4302533.2013868/1‑94 (MQ=60)
CAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCC                                                                                              >  SRR4302533.2909089/1‑94 (MQ=60)
 AGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGACG                                                                                             >  SRR4302533.5757185/1‑94 (MQ=60)
  GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGG                                                                                            <  SRR4302533.122449/94‑1 (MQ=60)
    GAAATTAGCATGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTG                                                                                          <  SRR4302533.4261599/94‑1 (MQ=60)
     AAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGT                                                                                         >  SRR4302533.3572245/1‑94 (MQ=60)
     AAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGT                                                                                         <  SRR4302533.5544716/94‑1 (MQ=60)
     AAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGGGT                                                                                         >  SRR4302533.1839617/1‑94 (MQ=60)
      AATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTT                                                                                        >  SRR4302533.4171528/1‑94 (MQ=60)
       ATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTG                                                                                       <  SRR4302533.924678/94‑1 (MQ=60)
        TTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGC                                                                                      <  SRR4302533.5354524/94‑1 (MQ=60)
         TAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCC                                                                                     >  SRR4302533.6226856/1‑94 (MQ=60)
            CAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTC                                                                                  >  SRR4302533.6090413/1‑94 (MQ=60)
            CAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTC                                                                                  >  SRR4302533.2573381/1‑94 (MQ=60)
            CAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTC                                                                                  <  SRR4302533.3940221/94‑1 (MQ=60)
             AGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCT                                                                                 >  SRR4302533.1532447/1‑94 (MQ=60)
             AGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCT                                                                                 <  SRR4302533.5898408/94‑1 (MQ=60)
                  AAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGT                                                                            <  SRR4302533.1034871/94‑1 (MQ=60)
                  AAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGT                                                                            >  SRR4302533.2314782/1‑94 (MQ=60)
                   AGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTT                                                                           >  SRR4302533.22931/1‑94 (MQ=60)
                    GGTCGGAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTG                                                                          <  SRR4302533.4704330/94‑1 (MQ=60)
                    GGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTG                                                                          >  SRR4302533.5355012/1‑94 (MQ=60)
                    GGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTG                                                                          <  SRR4302533.3113827/94‑1 (MQ=60)
                    GGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTG                                                                          >  SRR4302533.3271531/1‑94 (MQ=60)
                    GGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTG                                                                          >  SRR4302533.140664/1‑94 (MQ=60)
                      tcgcaagtgcggcttgtttcCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCG                                                                        <  SRR4302533.1137854/74‑1 (MQ=60)
                      TCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCG                                                                        >  SRR4302533.980341/1‑94 (MQ=60)
                      TCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGAGGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCG                                                                        <  SRR4302533.4687485/94‑1 (MQ=60)
                       CGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGG                                                                       >  SRR4302533.253516/1‑94 (MQ=60)
                       CGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGG                                                                       >  SRR4302533.2775537/1‑94 (MQ=60)
                       CGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGG                                                                       >  SRR4302533.116770/1‑94 (MQ=60)
                          AATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTC                                                                    <  SRR4302533.6286088/94‑1 (MQ=60)
                           ATGTCCGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCG                                                                   <  SRR4302533.1273108/94‑1 (MQ=60)
                           ATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCG                                                                   <  SRR4302533.1818153/94‑1 (MQ=60)
                             GTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCG                                                                 >  SRR4302533.3724880/1‑94 (MQ=60)
                              TACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGT                                                                <  SRR4302533.1394752/94‑1 (MQ=60)
                              TACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGT                                                                <  SRR4302533.4140915/94‑1 (MQ=60)
                               ACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTG                                                               >  SRR4302533.5373374/1‑94 (MQ=60)
                                CGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGG                                                              >  SRR4302533.5178750/1‑94 (MQ=60)
                                CGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGG                                                              >  SRR4302533.5221849/1‑94 (MQ=60)
                                 GCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGC                                                             >  SRR4302533.6272412/1‑94 (MQ=60)
                                  CGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCG                                                            <  SRR4302533.396646/94‑1 (MQ=60)
                                  CGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCGGTGGCG                                                            >  SRR4302533.4708630/1‑94 (MQ=60)
                                        TTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAA                                                      <  SRR4302533.5494871/94‑1 (MQ=60)
                                         TCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAG                                                     <  SRR4302533.2672167/94‑1 (MQ=60)
                                         TCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAG                                                     >  SRR4302533.3602592/1‑94 (MQ=60)
                                          CCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGC                                                    <  SRR4302533.4639738/94‑1 (MQ=60)
                                          CCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGC                                                    >  SRR4302533.4060398/1‑94 (MQ=60)
                                            CTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAG                                                  <  SRR4302533.4659635/94‑1 (MQ=60)
                                            CTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGTGGGAGCAG                                                  >  SRR4302533.2163455/1‑94 (MQ=60)
                                                CCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCGTGGCGGggggagaagggtt                                              >  SRR4302533.1627463/1‑81 (MQ=60)
                                                 CTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGGAGGGATT                                             >  SRR4302533.5211487/1‑94 (MQ=60)
                                                  TCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTG                                            <  SRR4302533.1124383/94‑1 (MQ=60)
                                                   CGACTGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGC                                           <  SRR4302533.4534765/94‑1 (MQ=60)
                                                    GACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTtgcgggggggggggagcacgtagtcgc                                          >  SRR4302533.2476401/1‑67 (MQ=60)
                                                    GACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCC                                          <  SRR4302533.5999560/94‑1 (MQ=60)
                                                    GACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCC                                          <  SRR4302533.6130293/94‑1 (MQ=60)
                                                    GACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGGGGAAGCAGGTATTGCC                                          >  SRR4302533.4776681/1‑94 (MQ=60)
                                                    GACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGACTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAACAGGTATTGCG                                          >  SRR4302533.4280429/1‑94 (MQ=60)
                                                     ACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCG                                         <  SRR4302533.2157433/94‑1 (MQ=60)
                                                      CCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGA                                        <  SRR4302533.3959405/94‑1 (MQ=60)
                                                      CCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGA                                        >  SRR4302533.4666106/1‑94 (MQ=60)
                                                        GATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATT                                      <  SRR4302533.3083424/94‑1 (MQ=60)
                                                        GATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATT                                      <  SRR4302533.5995416/94‑1 (MQ=60)
                                                        GATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATT                                      <  SRR4302533.4604076/94‑1 (MQ=60)
                                                         ATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTA                                     <  SRR4302533.3101257/94‑1 (MQ=60)
                                                          TATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTAC                                    <  SRR4302533.1481678/94‑1 (MQ=60)
                                                            TTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTG                                  <  SRR4302533.5279216/94‑1 (MQ=60)
                                                               TCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTA                               <  SRR4302533.5173749/94‑1 (MQ=60)
                                                                   CGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAG                           <  SRR4302533.5167166/94‑1 (MQ=60)
                                                                   CGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCCTGGCGGTGGAAGCAGGTGTTGCCGCTTACTGGTACAAG                           >  SRR4302533.4846579/1‑94 (MQ=60)
                                                                    GCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTGGCGGCTCGCGTGGTGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGT                          <  SRR4302533.6198112/94‑1 (MQ=60)
                                                                     CCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCCGGTATTGCCGATTCCTGGTACAAGCA                         >  SRR4302533.2438290/1‑94 (MQ=60)
                                                                      CCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGTACAAGTAT                        >  SRR4302533.5029572/1‑94 (MQ=60)
                                                                      CCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTAT                        <  SRR4302533.546200/94‑1 (MQ=60)
                                                                       CAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATG                       <  SRR4302533.111029/94‑1 (MQ=60)
                                                                       CAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATG                       >  SRR4302533.5925283/1‑94 (MQ=60)
                                                                       CAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATG                       >  SRR4302533.1049562/1‑94 (MQ=60)
                                                                        AGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGT                      >  SRR4302533.5726073/1‑94 (MQ=60)
                                                                        AGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGT                      >  SRR4302533.3056088/1‑94 (MQ=60)
                                                                           ATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGG                   <  SRR4302533.5553783/94‑1 (MQ=60)
                                                                              AGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCT                <  SRR4302533.4185092/94‑1 (MQ=60)
                                                                              AGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCT                <  SRR4302533.3859280/94‑1 (MQ=60)
                                                                               GGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTG               <  SRR4302533.476114/94‑1 (MQ=60)
                                                                               GGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTG               <  SRR4302533.5194200/94‑1 (MQ=60)
                                                                                GAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGA              <  SRR4302533.3874844/94‑1 (MQ=60)
                                                                                GAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGTTGGTCTGA              >  SRR4302533.2326137/1‑94 (MQ=60)
                                                                                GAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGTTGGTCTGA              >  SRR4302533.2060725/1‑94 (MQ=60)
                                                                                  ATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAA            <  SRR4302533.5476333/94‑1 (MQ=60)
                                                                                   TATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAG           <  SRR4302533.4359222/94‑1 (MQ=60)
                                                                                    ATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGG          <  SRR4302533.1038012/94‑1 (MQ=60)
                                                                                    ATCGGGAGCCGGTGTGGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGG          <  SRR4302533.1072822/94‑1 (MQ=60)
                                                                                     TCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGG         <  SRR4302533.6094224/94‑1 (MQ=60)
                                                                                     TCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGG         <  SRR4302533.4939324/94‑1 (MQ=60)
                                                                                       GGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGC       >  SRR4302533.6131626/1‑94 (MQ=60)
                                                                                        GGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGACGGTGGAAGCATGTATTGCCGATTACTGGTACAAGTATGTAGGTCTGAAAGGGGCA      <  SRR4302533.3215602/94‑1 (MQ=60)
                                                                                           GCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATT   <  SRR4302533.2563065/94‑1 (MQ=60)
                                                                                           GCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATT   >  SRR4302533.3694887/1‑94 (MQ=60)
                                                                                           GCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATT   <  SRR4302533.3446762/94‑1 (MQ=60)
                                                                                            CCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTT  >  SRR4302533.1493971/1‑94 (MQ=60)
                                                                                            CCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTG  >  SRR4302533.381298/1‑94 (MQ=60)
                                                                                            CCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTG  <  SRR4302533.3117003/94‑1 (MQ=60)
                                                                                            CCGGTGTTGCCTTCTAACGTTGCGGCTCGCGCGGCGGGGGAAGCAGGCATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTG  >  SRR4302533.2757382/1‑94 (MQ=60)
                                                                                            |                                                                                             
CAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTG  >  NC_000913/2579355‑2579540

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