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,870,363 T→C Q293R (CAG→CGG)  dgoA ← 2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,870,3630TC95.4% 740.5 / 17.9 240Q293R (CAG→CGG) dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (7/4);  new base C (123/106);  total (130/110)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

Only 100 of 250 total aligned reads displayed.

ACTGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCGC  >  NC_000913/3870272‑3870452
                                                                                           |                                                                                         
 cTGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCgg                                                                                         >  2:2303455/1‑93 (MQ=255)
  tGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACGGCGGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCggc                                                                                        >  2:2221196/1‑92 (MQ=255)
  tGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGt                                                                                        >  2:2742517/1‑93 (MQ=255)
     aGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCTCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCaa                                                                                     <  1:3451801/93‑1 (MQ=255)
     aGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCaa                                                                                     <  1:5737776/93‑1 (MQ=255)
      ggCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAAt                                                                                    <  2:5191303/93‑1 (MQ=255)
       gCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATc                                                                                   >  1:5038968/1‑93 (MQ=255)
       gCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATc                                                                                   >  1:3717851/1‑93 (MQ=255)
        cTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATcc                                                                                  >  1:1395429/1‑93 (MQ=255)
         tttAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATccc                                                                                 >  2:6439017/1‑93 (MQ=255)
         tttAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATccc                                                                                 >  2:4440098/1‑93 (MQ=255)
             aTTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCg                                                                             >  2:2550011/1‑93 (MQ=255)
             aTTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATACCGGCg                                                                             >  2:1267293/1‑93 (MQ=255)
               tttGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGAt                                                                           >  2:956482/1‑93 (MQ=255)
                ttGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATg                                                                          >  2:4963695/1‑93 (MQ=255)
                gtGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATg                                                                          >  2:6206719/2‑93 (MQ=255)
                 tGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATgg                                                                         <  1:459148/93‑1 (MQ=255)
                 tGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATgg                                                                         <  2:734944/93‑1 (MQ=255)
                  gTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGgc                                                                        >  1:291640/1‑93 (MQ=255)
                  gTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGgc                                                                        <  2:1907886/93‑1 (MQ=255)
                   tCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGgcg                                                                       <  1:920632/93‑1 (MQ=255)
                     gcgAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCt                                                                     >  2:3127616/1‑93 (MQ=255)
                     gcgAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGAGCt                                                                     >  2:2489167/1‑93 (MQ=255)
                          ccAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCtcatca                                                                >  2:5859056/1‑93 (MQ=255)
                          ccAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCtcatca                                                                >  1:4119498/1‑93 (MQ=255)
                            aGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAgg                                                              >  2:4644712/1‑93 (MQ=255)
                             ggCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGa                                                             <  2:2936397/93‑1 (MQ=255)
                             ggCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGa                                                             <  2:4230494/93‑1 (MQ=255)
                             ggCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGa                                                             >  1:2170388/1‑93 (MQ=255)
                             ggCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATAAGGa                                                             >  1:1541876/1‑93 (MQ=255)
                               ccgccCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATc                                                           <  1:1119351/93‑1 (MQ=255)
                                cgccCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCg                                                          <  1:676310/93‑1 (MQ=255)
                                cgccCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCg                                                          <  2:2163968/93‑1 (MQ=255)
                                cgccCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCg                                                          <  1:618109/93‑1 (MQ=255)
                                  cccATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGa                                                        <  2:2869895/93‑1 (MQ=255)
                                    cATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGAcc                                                      >  1:3147396/1‑93 (MQ=255)
                                    cATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGAcc                                                      <  1:4583197/93‑1 (MQ=255)
                                    cATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGAcc                                                      >  2:5760733/1‑93 (MQ=255)
                                    cATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGAcc                                                      >  1:1820364/1‑93 (MQ=255)
                                    cATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGAcc                                                      <  1:1422107/93‑1 (MQ=255)
                                       cAGTTGCCAGACCGGCGCATTCAGCGCTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGAccgcc                                                   <  2:3983994/93‑1 (MQ=255)
                                         gTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACcgccgc                                                 >  2:6004460/1‑93 (MQ=255)
                                         gTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACcgccgc                                                 >  2:1026158/1‑93 (MQ=255)
                                           tGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGa                                               <  2:75998/93‑1 (MQ=255)
                                            gCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGaa                                              >  1:4932552/1‑92 (MQ=255)
                                            gCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGAt                                              >  1:1889493/1‑93 (MQ=255)
                                             ccAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATa                                             >  1:1156625/1‑93 (MQ=255)
                                              cAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAg                                            <  2:174190/93‑1 (MQ=255)
                                              cAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGAACGCCGCGATAg                                            >  2:1963621/1‑93 (MQ=255)
                                               aGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGa                                           >  1:513586/1‑93 (MQ=255)
                                                gACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGaa                                          <  1:387316/93‑1 (MQ=255)
                                                gACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGaa                                          >  1:4833891/1‑93 (MQ=255)
                                                 aCCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAAt                                         <  1:4527244/93‑1 (MQ=255)
                                                  ccGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATc                                        <  1:4432789/93‑1 (MQ=255)
                                                   cGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATcc                                       <  1:2390849/93‑1 (MQ=255)
                                                    ggCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCg                                      <  1:883225/93‑1 (MQ=255)
                                                    cgcgcATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCg                                      <  2:1955278/92‑1 (MQ=255)
                                                     gcgcATTCAGCACTTTTCCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCgg                                     <  2:1180194/93‑1 (MQ=255)
                                                      cgcATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGTATCGGACCGCCTCGATAGAATCCGgc                                    <  2:1351994/93‑1 (MQ=255)
                                                         aTTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGGCGCGATAGAATCCGGcgcg                                 >  1:4312450/1‑93 (MQ=255)
                                                         aTTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGcgcg                                 >  2:3672965/1‑93 (MQ=255)
                                                         aTTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGcgcg                                 >  2:3825596/1‑93 (MQ=255)
                                                         aTTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGcgcg                                 >  1:5416932/1‑93 (MQ=255)
                                                          ttCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGa                                >  2:3162644/1‑93 (MQ=255)
                                                           tCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGAt                               >  1:5347985/1‑93 (MQ=255)
                                                            cAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATa                              >  2:2378523/1‑93 (MQ=255)
                                                            cAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATa                              <  2:4524476/93‑1 (MQ=255)
                                                             aGCTCTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATAc                             <  2:4014902/93‑1 (MQ=255)
                                                             aGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATAc                             <  1:4315779/93‑1 (MQ=255)
                                                              gCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACa                            <  1:6017166/93‑1 (MQ=255)
                                                              gCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACa                            >  2:3878734/1‑93 (MQ=255)
                                                                aCTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATc                          >  2:637188/1‑93 (MQ=255)
                                                                aCTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATc                          >  1:3593664/1‑93 (MQ=255)
                                                                aCTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATc                          <  1:1661625/93‑1 (MQ=255)
                                                                aCTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATc                          >  1:2113585/1‑93 (MQ=255)
                                                                  tttACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCAc                        >  1:6117601/1‑93 (MQ=255)
                                                                  tttACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCAc                        <  1:4453867/93‑1 (MQ=255)
                                                                   ttACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACt                       >  1:2274221/1‑93 (MQ=255)
                                                                   ttACCATTGATATCCCATAACTCACGGTCACTCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACt                       <  1:3373541/93‑1 (MQ=255)
                                                                    tACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACtt                      <  2:1456843/93‑1 (MQ=255)
                                                                     aCCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTg                     >  1:2484842/1‑93 (MQ=255)
                                                                           gATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATaa               >  1:6069461/1‑93 (MQ=255)
                                                                              ataccATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTc            >  2:415747/4‑93 (MQ=255)
                                                                               tCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCa           <  1:2861300/93‑1 (MQ=255)
                                                                               tCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCa           >  2:5587782/1‑93 (MQ=255)
                                                                                cccATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCAt          >  2:5448251/1‑93 (MQ=255)
                                                                                cccATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCAt          <  1:6395784/93‑1 (MQ=255)
                                                                                cccATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCAt          <  1:1319154/93‑1 (MQ=255)
                                                                                 ccATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCAtt         <  2:1806638/93‑1 (MQ=255)
                                                                                 ccATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCAtt         <  1:168609/93‑1 (MQ=255)
                                                                                 ccATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCAtt         <  2:856466/93‑1 (MQ=255)
                                                                                  cATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTg        >  1:5212294/1‑93 (MQ=255)
                                                                                  cATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTg        >  1:2327979/1‑93 (MQ=255)
                                                                                    tAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGAt      <  2:5130539/93‑1 (MQ=255)
                                                                                    tAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGAt      >  2:1855347/1‑93 (MQ=255)
                                                                                    tAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGAt      <  2:2732672/93‑1 (MQ=255)
                                                                                      aCGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATgc    <  2:1510190/93‑1 (MQ=255)
                                                                                      aCGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATgc    >  1:5940049/1‑93 (MQ=255)
                                                                                        gCCCGGTCAATCCCGGCGATGGCGCTCGTCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATgcgc  <  1:4070529/93‑1 (MQ=255)
                                                                                        gCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATgcgc  >  2:3496469/1‑93 (MQ=255)
                                                                                           |                                                                                         
ACTGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCGC  >  NC_000913/3870272‑3870452

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ 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
CAACTGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCGCGATG  >  NC_000913/3870270‑3870456
                                                                                             |                                                                                             
CAACTGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCC                                                                                               >  SRR4302537.5868248/1‑94 (MQ=60)
CAACTGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCC                                                                                               >  SRR4302537.1181713/1‑94 (MQ=60)
 AACTGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCG                                                                                              >  SRR4302537.3442198/1‑94 (MQ=60)
 AACTGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCG                                                                                              <  SRR4302537.2846995/94‑1 (MQ=60)
     GTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCA                                                                                          >  SRR4302537.3792274/1‑94 (MQ=60)
     GTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCA                                                                                          >  SRR4302537.6494346/1‑94 (MQ=60)
       AGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCTCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAAT                                                                                        <  SRR4302537.3476213/94‑1 (MQ=60)
       AGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAAT                                                                                        <  SRR4302537.5777059/94‑1 (MQ=60)
       AGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAAT                                                                                        <  SRR4302537.2414901/94‑1 (MQ=60)
        GGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATC                                                                                       >  SRR4302537.5073658/1‑94 (MQ=60)
        GGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATC                                                                                       >  SRR4302537.2003016/1‑94 (MQ=60)
        GGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATC                                                                                       >  SRR4302537.3743981/1‑94 (MQ=60)
         GCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCC                                                                                      >  SRR4302537.1405159/1‑94 (MQ=60)
         GCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCC                                                                                      <  SRR4302537.3375093/94‑1 (MQ=60)
         GCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCC                                                                                      >  SRR4302537.4627784/1‑94 (MQ=60)
          CTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCC                                                                                     <  SRR4302537.6034738/94‑1 (MQ=60)
          CTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCC                                                                                     >  SRR4302537.2255330/1‑94 (MQ=60)
          CTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCC                                                                                     <  SRR4302537.96883/94‑1 (MQ=60)
          CTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCC                                                                                     <  SRR4302537.5135451/94‑1 (MQ=60)
           TTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCG                                                                                    >  SRR4302537.4619216/1‑94 (MQ=60)
            TTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGG                                                                                   >  SRR4302537.3654372/1‑94 (MQ=60)
                TTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTCCCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGAT                                                                               <  SRR4302537.3626201/94‑1 (MQ=60)
                TTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGAT                                                                               <  SRR4302537.2384262/94‑1 (MQ=60)
                 TTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATG                                                                              <  SRR4302537.51732/94‑1 (MQ=60)
                   TGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGC                                                                            >  SRR4302537.293669/1‑94 (MQ=60)
                   TGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGC                                                                            <  SRR4302537.462310/94‑1 (MQ=60)
                   TGTCGCGAACCAGGCCGCCCATCAGTAGCCAGACCGGCGCATTCAGCACTTTACCTTTGATAACCCATAACGCCCCGTCAATCCCGGCGATGGC                                                                            <  SRR4302537.3590181/94‑1 (MQ=60)
                     TCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGC                                                                          <  SRR4302537.926985/94‑1 (MQ=60)
                      CGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCT                                                                         >  SRR4302537.1333101/1‑94 (MQ=60)
                        CGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCA                                                                       >  SRR4302537.1793468/1‑94 (MQ=60)
                           ACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCA                                                                    >  SRR4302537.4148263/1‑94 (MQ=60)
                            CCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAG                                                                   <  SRR4302537.5982412/94‑1 (MQ=60)
                             CAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGG                                                                  >  SRR4302537.2367083/1‑94 (MQ=60)
                             CAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGG                                                                  <  SRR4302537.4482622/94‑1 (MQ=60)
                              AGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGA                                                                 <  SRR4302537.704482/94‑1 (MQ=60)
                              AGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGA                                                                 >  SRR4302537.2185683/1‑94 (MQ=60)
                              AGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATAAGGA                                                                 >  SRR4302537.1552638/1‑94 (MQ=60)
                               GGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGAT                                                                >  SRR4302537.1693369/1‑94 (MQ=60)
                               GGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGAT                                                                >  SRR4302537.4585728/1‑94 (MQ=60)
                                GCCGCCCATCAGTTTCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATC                                                               <  SRR4302537.3091688/94‑1 (MQ=60)
                                GCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATC                                                               <  SRR4302537.417698/94‑1 (MQ=60)
                                 CCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGTATCG                                                              >  SRR4302537.6169996/1‑94 (MQ=60)
                                 CCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCG                                                              <  SRR4302537.1127099/94‑1 (MQ=60)
                                  CGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGG                                                             >  SRR4302537.6436946/1‑94 (MQ=60)
                                  CGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGG                                                             <  SRR4302537.622315/94‑1 (MQ=60)
                                  CGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGG                                                             <  SRR4302537.6158393/94‑1 (MQ=60)
                                  CGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGG                                                             >  SRR4302537.2172560/1‑94 (MQ=60)
                                  CGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGG                                                             <  SRR4302537.680943/94‑1 (MQ=60)
                                   CCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGA                                                            <  SRR4302537.910034/94‑1 (MQ=60)
                                    CCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGAC                                                           >  SRR4302537.539006/1‑94 (MQ=60)
                                     CCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACC                                                          >  SRR4302537.3169475/1‑94 (MQ=60)
                                     CCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACC                                                          >  SRR4302537.1833050/1‑94 (MQ=60)
                                     CCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACC                                                          >  SRR4302537.5552585/1‑94 (MQ=60)
                                     CCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACC                                                          <  SRR4302537.1556446/94‑1 (MQ=60)
                                      CATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCG                                                         <  SRR4302537.1432010/94‑1 (MQ=60)
                                      CATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCG                                                         <  SRR4302537.4614898/94‑1 (MQ=60)
                                             TGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGAT                                                  >  SRR4302537.4409191/1‑94 (MQ=60)
                                             TGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGAT                                                  >  SRR4302537.1902626/1‑94 (MQ=60)
                                             TGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGAT                                                  >  SRR4302537.4890455/1‑94 (MQ=60)
                                             TGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGAA                                                  >  SRR4302537.4966550/1‑94 (MQ=60)
                                              GCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATA                                                 >  SRR4302537.1164671/1‑94 (MQ=60)
                                                CAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGA                                               >  SRR4302537.517110/1‑94 (MQ=60)
                                                 AGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAA                                              >  SRR4302537.4867269/1‑94 (MQ=60)
                                                  GACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAAT                                             <  SRR4302537.4194469/94‑1 (MQ=60)
                                                  GACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAAT                                             <  SRR4302537.389991/94‑1 (MQ=60)
                                                   ACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATC                                            <  SRR4302537.4558591/94‑1 (MQ=60)
                                                   ACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATC                                            >  SRR4302537.1203242/1‑94 (MQ=60)
                                                    CCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCC                                           <  SRR4302537.4463400/94‑1 (MQ=60)
                                                     CGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCG                                          >  SRR4302537.3928928/1‑94 (MQ=60)
                                                     CGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCG                                          <  SRR4302537.6145573/94‑1 (MQ=60)
                                                     CGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCG                                          <  SRR4302537.2407685/94‑1 (MQ=60)
                                                      GGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGG                                         <  SRR4302537.889302/94‑1 (MQ=60)
                                                      GGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGG                                         <  SRR4302537.1500470/94‑1 (MQ=60)
                                                       GCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGC                                        <  SRR4302537.1179463/94‑1 (MQ=60)
                                                        CGCATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCG                                       >  SRR4302537.2344102/1‑94 (MQ=60)
                                                          CATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGGCGCGATAGAATCCGGCGCG                                     >  SRR4302537.4342420/1‑94 (MQ=60)
                                                          CATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCG                                     >  SRR4302537.5453859/1‑94 (MQ=60)
                                                          CATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCG                                     <  SRR4302537.4619731/94‑1 (MQ=60)
                                                          CATTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCG                                     <  SRR4302537.5270801/94‑1 (MQ=60)
                                                            TTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGAT                                   <  SRR4302537.5741325/94‑1 (MQ=60)
                                                            TTCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGAT                                   >  SRR4302537.5384627/1‑94 (MQ=60)
                                                             TCAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATA                                  >  SRR4302537.6253728/1‑94 (MQ=60)
                                                              CAGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATAC                                 >  SRR4302537.1698233/1‑94 (MQ=60)
                                                               AGCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACA                                <  SRR4302537.4345779/94‑1 (MQ=60)
                                                                GCACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACAT                               <  SRR4302537.6058505/94‑1 (MQ=60)
                                                                 CACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATC                              >  SRR4302537.3504778/1‑94 (MQ=60)
                                                                 CACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATC                              >  SRR4302537.2128478/1‑94 (MQ=60)
                                                                 CACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATC                              >  SRR4302537.929037/1‑94 (MQ=60)
                                                                 CACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATC                              >  SRR4302537.1585349/1‑94 (MQ=60)
                                                                 CACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATC                              >  SRR4302537.3618999/1‑94 (MQ=60)
                                                                  ACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCA                             >  SRR4302537.4536653/1‑94 (MQ=60)
                                                                  ACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCA                             <  SRR4302537.882684/94‑1 (MQ=60)
                                                                  ACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCA                             >  SRR4302537.1430575/1‑94 (MQ=60)
                                                                  ACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCA                             <  SRR4302537.4343618/94‑1 (MQ=60)
                                                                  ACTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCA                             <  SRR4302537.1673268/94‑1 (MQ=60)
                                                                   CTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCAC                            >  SRR4302537.4921110/1‑94 (MQ=60)
                                                                   CTTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCAC                            >  SRR4302537.6159646/1‑94 (MQ=60)
                                                                    TTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACT                           <  SRR4302537.1136545/94‑1 (MQ=60)
                                                                    TTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACT                           <  SRR4302537.4484626/94‑1 (MQ=60)
                                                                    TTTACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACT                           >  SRR4302537.2290279/1‑94 (MQ=60)
                                                                     TTACCATTGATATCCCATAACTCACGGTCACTCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTT                          <  SRR4302537.3397449/94‑1 (MQ=60)
                                                                      TACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTG                         >  SRR4302537.2900340/1‑94 (MQ=60)
                                                                      TACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTG                         >  SRR4302537.2502335/1‑94 (MQ=60)
                                                                       ACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGC                        >  SRR4302537.1926711/1‑94 (MQ=60)
                                                                       ACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGC                        >  SRR4302537.5881488/1‑94 (MQ=60)
                                                                       ACCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGC                        >  SRR4302537.340303/1‑94 (MQ=60)
                                                                        CCTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCC                       >  SRR4302537.2961982/1‑94 (MQ=60)
                                                                         CTTTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCA                      >  SRR4302537.392596/1‑94 (MQ=60)
                                                                           TTGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATA                    >  SRR4302537.6086232/1‑94 (MQ=60)
                                                                            tgcgATCCCATAACGCCCGGTCCATCTCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGCTAGAATCCGcagaggtcactcacttgccctaa                   >  SRR4302537.6111158/5‑71 (MQ=60)
                                                                            TGATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAA                   >  SRR4302537.6111162/1‑94 (MQ=60)
                                                                             GATATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAG                  >  SRR4302537.1877871/1‑94 (MQ=60)
                                                                               TATCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTC                <  SRR4302537.527331/94‑1 (MQ=60)
                                                                                 TCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCAT              <  SRR4302537.2881363/94‑1 (MQ=60)
                                                                                 TCCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCAT              <  SRR4302537.1430395/94‑1 (MQ=60)
                                                                                  CCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATT             <  SRR4302537.6439339/94‑1 (MQ=60)
                                                                                  CCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATT             >  SRR4302537.2249525/1‑94 (MQ=60)
                                                                                  CCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATT             <  SRR4302537.1328349/94‑1 (MQ=60)
                                                                                  CCCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATT             <  SRR4302537.4099328/94‑1 (MQ=60)
                                                                                   CCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTG            >  SRR4302537.5248133/1‑94 (MQ=60)
                                                                                   CCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTG            <  SRR4302537.2810157/94‑1 (MQ=60)
                                                                                   CCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTG            >  SRR4302537.2344396/1‑94 (MQ=60)
                                                                                   CCATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTG            <  SRR4302537.169812/94‑1 (MQ=60)
                                                                                     ATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGAT          >  SRR4302537.2598450/1‑94 (MQ=60)
                                                                                     ATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGAT          <  SRR4302537.4997421/94‑1 (MQ=60)
                                                                                     ATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGAT          >  SRR4302537.451897/1‑94 (MQ=60)
                                                                                     ATAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGAT          <  SRR4302537.5213547/94‑1 (MQ=60)
                                                                                      TAACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATG         <  SRR4302537.3413136/94‑1 (MQ=60)
                                                                                       AACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGC        >  SRR4302537.5980819/1‑94 (MQ=60)
                                                                                        ACGCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCG       >  SRR4302537.4352011/1‑94 (MQ=60)
                                                                                          GCCCGGTCAATCCCGGCGATGGCGCTCGTCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCGCG     <  SRR4302537.4098956/94‑1 (MQ=60)
                                                                                          GCCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCGCG     >  SRR4302537.2800961/1‑94 (MQ=60)
                                                                                           CCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCGCGA    >  SRR4302537.2331979/1‑94 (MQ=60)
                                                                                           CCCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCGCGA    >  SRR4302537.1704937/1‑94 (MQ=60)
                                                                                            CCGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCGCGAT   <  SRR4302537.3268618/94‑1 (MQ=60)
                                                                                             CGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCGCGATG  <  SRR4302537.5169472/94‑1 (MQ=60)
                                                                                             |                                                                                             
CAACTGTAGGCTTTAATTTTGTCGCGAACCAGGCCGCCCATCAGTTGCCAGACCGGCGCATTCAGCACTTTACCTTTGATATCCCATAACGCCTGGTCAATCCCGGCGATGGCGCTCATCAGGATCGGACCGCCGCGATAGAATCCGGCGCGATACATCACTTGCCATAAGTCATTGATGCGCGATG  >  NC_000913/3870270‑3870456

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