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 3,873,493 T→C C124C (TGT→TGC)  yidX → putative replicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,873,4930TC99.6% 918.0 / ‑5.1 281C124C (TGT→TGC) yidXputative replicase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (132/147);  total (134/147)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 285 total aligned reads displayed.

CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAT  >  NC_000913/3873403‑3873581
                                                                                          |                                                                                        
 ggCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCAcc                                                                                       <  1:1216672/93‑1 (MQ=255)
  gCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCa                                                                                      >  2:4795915/1‑93 (MQ=255)
   cacaATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAg                                                                                     <  2:6245476/93‑1 (MQ=255)
     caATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGt                                                                                   >  1:5712947/1‑93 (MQ=255)
       atatCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCa                                                                                 >  2:3902194/1‑93 (MQ=255)
        tatCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCAc                                                                                >  2:3654354/1‑93 (MQ=255)
           cGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTgg                                                                             <  2:277557/93‑1 (MQ=255)
            gCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGa                                                                            >  2:2908243/1‑93 (MQ=255)
             cTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaa                                                                           >  1:3575674/1‑93 (MQ=255)
              tttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGaaa                                                                          >  2:5064340/1‑93 (MQ=255)
              tttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaa                                                                          >  1:2460951/1‑93 (MQ=255)
               ttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaaa                                                                         <  2:4280749/93‑1 (MQ=255)
               ttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaaa                                                                         >  1:4506536/1‑93 (MQ=255)
                tactactacTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAAc                                                                        <  2:2297702/93‑1 (MQ=255)
                 actactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAAcc                                                                       >  1:535534/1‑93 (MQ=255)
                  ctactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCt                                                                      >  2:151960/1‑93 (MQ=255)
                   tactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTa                                                                     >  2:3310516/1‑93 (MQ=255)
                   tactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTa                                                                     <  2:4372987/93‑1 (MQ=255)
                   tactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTa                                                                     >  2:5821533/1‑93 (MQ=255)
                   tactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTa                                                                     <  2:909541/93‑1 (MQ=255)
                    actactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaa                                                                    <  2:3137535/93‑1 (MQ=255)
                    actactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaa                                                                    >  1:4841985/1‑93 (MQ=255)
                    actactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaa                                                                    <  1:3015263/93‑1 (MQ=255)
                     ctactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaa                                                                   <  1:4151909/93‑1 (MQ=255)
                      tactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaaa                                                                  <  2:3112105/93‑1 (MQ=255)
                      tactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaaa                                                                  <  1:3790880/93‑1 (MQ=255)
                       actGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAg                                                                 <  2:31938/93‑1 (MQ=255)
                       actGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAg                                                                 <  1:4972388/93‑1 (MQ=255)
                         tGGAATAGCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCGTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                               <  2:2967860/93‑1 (MQ=255)
                         tGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                               <  1:5964561/93‑1 (MQ=255)
                         tGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                               >  2:5453711/1‑93 (MQ=255)
                         tGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                               >  2:4445564/1‑93 (MQ=255)
                          ggAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGca                                                              >  2:1462753/1‑93 (MQ=255)
                            aaTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGcaca                                                            <  2:5657541/93‑1 (MQ=255)
                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                           >  2:5386376/1‑93 (MQ=255)
                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                           <  1:2573526/93‑1 (MQ=255)
                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                           >  2:2320103/1‑93 (MQ=255)
                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                           >  2:1987142/1‑93 (MQ=255)
                               aaCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATc                                                         >  2:3557446/1‑93 (MQ=255)
                                 cGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCa                                                       <  1:1367095/93‑1 (MQ=255)
                                 cGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCa                                                       <  2:553554/93‑1 (MQ=255)
                                 cGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCa                                                       >  1:3250798/1‑93 (MQ=255)
                                  gaagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAACCCTAAAAGGCACAATCCAt                                                      <  2:6081728/93‑1 (MQ=255)
                                  gaagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAt                                                      <  2:3409588/93‑1 (MQ=255)
                                  gaagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAt                                                      <  1:1006908/93‑1 (MQ=255)
                                   aagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATg                                                     <  2:4491087/93‑1 (MQ=255)
                                   aagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATg                                                     <  2:4331098/93‑1 (MQ=255)
                                    agaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCGCCAGGTCACTGGAAAACCTAAAAGGCACAATCCATga                                                    <  1:1911739/93‑1 (MQ=255)
                                    agaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATga                                                    <  1:3333472/93‑1 (MQ=255)
                                     gaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgag                                                   <  2:1319946/93‑1 (MQ=255)
                                     gaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgag                                                   <  1:1173589/93‑1 (MQ=255)
                                     gaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgag                                                   >  2:523107/1‑93 (MQ=255)
                                        gaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGCAAACCTAAAAGGCACAATCCATGAGaaa                                                >  2:2637391/1‑93 (MQ=255)
                                        gaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaa                                                <  1:1289599/93‑1 (MQ=255)
                                         aaaaGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaaa                                               <  1:3663079/93‑1 (MQ=255)
                                             ggCACTGGCAACTCATTATGGTTTTGCCTGTGGTGTTCAACACGGCATCAGGTCACTGGAGAACCTAAAAGGCACAATCCATGAGAAAAATaa                                           <  2:1006385/93‑1 (MQ=255)
                                             ggCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaa                                           >  1:2915618/1‑93 (MQ=255)
                                              gCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaa                                          <  1:2712341/93‑1 (MQ=255)
                                              gCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaa                                          >  2:778204/1‑93 (MQ=255)
                                               cACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaaa                                         <  1:39940/93‑1 (MQ=255)
                                                aCTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAACAGGCAGAATCCATGAGAAAAATaaaaa                                        >  1:5012861/1‑93 (MQ=255)
                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                       >  1:4054795/1‑93 (MQ=255)
                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                       >  2:4201517/1‑93 (MQ=255)
                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                       >  1:5223576/1‑93 (MQ=255)
                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                       <  2:5091624/93‑1 (MQ=255)
                                                     cAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATgg                                   <  1:3117835/93‑1 (MQ=255)
                                                     cAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATgg                                   <  1:5845844/93‑1 (MQ=255)
                                                     cAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATgg                                   <  1:5861164/93‑1 (MQ=255)
                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                  <  2:1749272/93‑1 (MQ=255)
                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                  <  2:5016626/93‑1 (MQ=255)
                                                       aCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAt                                 >  2:4903249/1‑93 (MQ=255)
                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                <  1:5677531/93‑1 (MQ=255)
                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                <  1:2631129/93‑1 (MQ=255)
                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                >  2:2565612/1‑93 (MQ=255)
                                                         gCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTa                               >  1:3760179/1‑93 (MQ=255)
                                                         gCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTa                               >  2:3491438/1‑93 (MQ=255)
                                                         gCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTa                               >  1:1613218/1‑93 (MQ=255)
                                                          cATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTAc                              <  2:5712947/93‑1 (MQ=255)
                                                           aTTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACt                             <  2:4668126/93‑1 (MQ=255)
                                                            ttATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTc                            <  2:5786152/93‑1 (MQ=255)
                                                             tATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCa                           >  2:2403604/1‑93 (MQ=255)
                                                              aTGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCaa                          >  1:6246164/1‑93 (MQ=255)
                                                              aTGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCaa                          >  2:3972118/1‑93 (MQ=255)
                                                               tGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCaaa                         <  2:4088332/93‑1 (MQ=255)
                                                                ggTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAg                        <  1:2752796/93‑1 (MQ=255)
                                                                  ttttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGt                      >  1:2264055/1‑93 (MQ=255)
                                                                   tttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTg                     <  1:4458029/93‑1 (MQ=255)
                                                                    ttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGa                    <  1:738372/93‑1 (MQ=255)
                                                                          gtggtgTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGt              >  1:155145/1‑93 (MQ=255)
                                                                                 tCAACACTGCACCAGGTCATTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAcca       >  1:1694684/1‑93 (MQ=255)
                                                                                  cAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccat      >  1:6165059/1‑93 (MQ=255)
                                                                                   agcacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATCACTCAAAGGTGATGGCGTTCTAccatc     <  2:2079115/91‑1 (MQ=255)
                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatc     >  2:2246300/1‑93 (MQ=255)
                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatc     <  1:1504627/93‑1 (MQ=255)
                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatc     <  2:2669444/93‑1 (MQ=255)
                                                                                     cacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcca   <  1:3966741/93‑1 (MQ=255)
                                                                                     cacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcca   >  1:4409372/1‑93 (MQ=255)
                                                                                     cacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcca   >  2:1117732/1‑93 (MQ=255)
                                                                                     cacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcca   <  2:1868444/93‑1 (MQ=255)
                                                                                      acTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatccat  >  1:2116870/1‑93 (MQ=255)
                                                                                          |                                                                                        
CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAT  >  NC_000913/3873403‑3873581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAA  >  NC_000913/3873400‑3873586
                                                                                             |                                                                                             
TACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGC                                                                                               <  SRR4302533.4261765/94‑1 (MQ=60)
TACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGC                                                                                               >  SRR4302533.6072731/1‑94 (MQ=60)
   CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAATGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACC                                                                                            <  SRR4302533.3154002/94‑1 (MQ=60)
   CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACC                                                                                            >  SRR4302533.4129910/1‑94 (MQ=60)
   CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACC                                                                                            >  SRR4302533.1457503/1‑94 (MQ=60)
    GGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCA                                                                                           <  SRR4302533.1223981/94‑1 (MQ=60)
    GGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCA                                                                                           <  SRR4302533.1287222/94‑1 (MQ=60)
      CACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGG                                                                                         <  SRR4302533.1578834/94‑1 (MQ=60)
       ACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGT                                                                                        >  SRR4302533.5747456/1‑94 (MQ=60)
       ACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGT                                                                                        <  SRR4302533.1797343/94‑1 (MQ=60)
         AATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCA                                                                                      >  SRR4302533.591528/1‑94 (MQ=60)
           TATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACT                                                                                    >  SRR4302533.3764077/1‑94 (MQ=60)
           TATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACT                                                                                    <  SRR4302533.1628480/94‑1 (MQ=60)
            ATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTG                                                                                   <  SRR4302533.1121521/94‑1 (MQ=60)
             TCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGG                                                                                  >  SRR4302533.4076078/1‑94 (MQ=60)
              CGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGA                                                                                 <  SRR4302533.1372202/94‑1 (MQ=60)
               GCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAA                                                                                >  SRR4302533.3597714/1‑94 (MQ=60)
                CTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAA                                                                               >  SRR4302533.2476136/1‑94 (MQ=60)
                 TTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAA                                                                              >  SRR4302533.4533831/1‑94 (MQ=60)
                  TTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAAC                                                                             <  SRR4302533.623656/94‑1 (MQ=60)
                   TACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACC                                                                            >  SRR4302533.538715/1‑94 (MQ=60)
                   TACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACC                                                                            <  SRR4302533.5307031/94‑1 (MQ=60)
                   TACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACC                                                                            <  SRR4302533.5896170/94‑1 (MQ=60)
                      TACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAA                                                                         >  SRR4302533.4871330/1‑94 (MQ=60)
                       ACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAA                                                                        <  SRR4302533.3033704/94‑1 (MQ=60)
                       ACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAA                                                                        <  SRR4302533.3805767/94‑1 (MQ=60)
                       ACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAA                                                                        >  SRR4302533.449103/1‑94 (MQ=60)
                        CTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAA                                                                       <  SRR4302533.4177065/94‑1 (MQ=60)
                        CTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAA                                                                       <  SRR4302533.3357333/94‑1 (MQ=60)
                         TACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAG                                                                      <  SRR4302533.850641/94‑1 (MQ=60)
                         TACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAG                                                                      <  SRR4302533.3814092/94‑1 (MQ=60)
                         TACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAG                                                                      <  SRR4302533.4493646/94‑1 (MQ=60)
                          ACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGG                                                                     <  SRR4302533.198294/94‑1 (MQ=60)
                          ACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGG                                                                     <  SRR4302533.5002536/94‑1 (MQ=60)
                          ACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGG                                                                     >  SRR4302533.4221507/1‑94 (MQ=60)
                           CTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGC                                                                    >  SRR4302533.2541161/1‑94 (MQ=60)
                            TGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCA                                                                   <  SRR4302533.6000653/94‑1 (MQ=60)
                             GGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCAC                                                                  >  SRR4302533.1212034/1‑94 (MQ=60)
                             GGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCAC                                                                  <  SRR4302533.3113446/94‑1 (MQ=60)
                                ATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAAT                                                               <  SRR4302533.2589357/94‑1 (MQ=60)
                                ATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAAT                                                               <  SRR4302533.6268432/94‑1 (MQ=60)
                                  AACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATga                                                             >  SRR4302533.3247438/1‑92 (MQ=60)
                                   ACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCA                                                            <  SRR4302533.6100585/94‑1 (MQ=60)
                                   ACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCA                                                            >  SRR4302533.3270846/1‑94 (MQ=60)
                                    CGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAT                                                           <  SRR4302533.3852116/94‑1 (MQ=60)
                                    CGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAT                                                           >  SRR4302533.4362671/1‑94 (MQ=60)
                                    CGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAT                                                           >  SRR4302533.3679590/1‑94 (MQ=60)
                                    CGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAT                                                           >  SRR4302533.3639073/1‑94 (MQ=60)
                                    CGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAT                                                           <  SRR4302533.1375313/94‑1 (MQ=60)
                                     GAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATGCATG                                                          >  SRR4302533.4446694/1‑94 (MQ=60)
                                     GAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATG                                                          <  SRR4302533.1012799/94‑1 (MQ=60)
                                     GAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATG                                                          <  SRR4302533.4037908/94‑1 (MQ=60)
                                     GAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATG                                                          >  SRR4302533.3103972/1‑94 (MQ=60)
                                     GAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTAAACACTGAACCGGGTCACTGGAAAACCTAAAAGGAACATTCAATG                                                          >  SRR4302533.4954307/1‑94 (MQ=60)
                                      AAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGA                                                         >  SRR4302533.2866661/1‑94 (MQ=60)
                                      AAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGA                                                         <  SRR4302533.1400834/94‑1 (MQ=60)
                                       AGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCGCCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAG                                                        <  SRR4302533.1923402/94‑1 (MQ=60)
                                       AGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAG                                                        <  SRR4302533.4684661/94‑1 (MQ=60)
                                       AGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAG                                                        <  SRR4302533.3354070/94‑1 (MQ=60)
                                       AGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAG                                                        <  SRR4302533.2943968/94‑1 (MQ=60)
                                       AGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAG                                                        <  SRR4302533.2488681/94‑1 (MQ=60)
                                       AGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAG                                                        <  SRR4302533.2595969/94‑1 (MQ=60)
                                        GAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGA                                                       >  SRR4302533.5003881/1‑94 (MQ=60)
                                        GAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGA                                                       <  SRR4302533.1180632/94‑1 (MQ=60)
                                        GAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGA                                                       >  SRR4302533.1268324/1‑94 (MQ=60)
                                         AAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAA                                                      >  SRR4302533.4140602/1‑94 (MQ=60)
                                           GAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAA                                                    <  SRR4302533.1297361/94‑1 (MQ=60)
                                            AAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAA                                                   >  SRR4302533.6139816/1‑94 (MQ=60)
                                            AAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAA                                                   <  SRR4302533.3685597/94‑1 (MQ=60)
                                             AAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAAT                                                  >  SRR4302533.3784679/1‑94 (MQ=60)
                                             AAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAAT                                                  >  SRR4302533.2647793/1‑94 (MQ=60)
                                             AAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAAT                                                  >  SRR4302533.867067/1‑94 (MQ=60)
                                              AAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATA                                                 <  SRR4302533.2781508/94‑1 (MQ=60)
                                              AAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATA                                                 <  SRR4302533.737576/94‑1 (MQ=60)
                                              AAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATA                                                 <  SRR4302533.1546126/94‑1 (MQ=60)
                                               AGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAA                                                >  SRR4302533.2933498/1‑94 (MQ=60)
                                                 GCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAA                                              <  SRR4302533.2728995/94‑1 (MQ=60)
                                                  CACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAACAGGCAGAATCCATGAGAAAAATAAAAA                                             >  SRR4302533.5043234/1‑94 (MQ=60)
                                                  CACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAA                                             <  SRR4302533.40208/94‑1 (MQ=60)
                                                   ACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAC                                            >  SRR4302533.5255002/1‑94 (MQ=60)
                                                   ACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAC                                            >  SRR4302533.4079482/1‑94 (MQ=60)
                                                     TGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACAT                                          <  SRR4302533.5530586/94‑1 (MQ=60)
                                                       GCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGG                                        <  SRR4302533.234908/94‑1 (MQ=60)
                                                       GCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGG                                        <  SRR4302533.1760541/94‑1 (MQ=60)
                                                        CAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGA                                       <  SRR4302533.5896603/94‑1 (MQ=60)
                                                        CAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGA                                       <  SRR4302533.3137006/94‑1 (MQ=60)
                                                        CAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGA                                       <  SRR4302533.5881185/94‑1 (MQ=60)
                                                         AACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAT                                      <  SRR4302533.5343009/94‑1 (MQ=60)
                                                         AACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAT                                      >  SRR4302533.1544928/1‑94 (MQ=60)
                                                           CGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTA                                    >  SRR4302533.5022732/1‑94 (MQ=60)
                                                           CGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTA                                    >  SRR4302533.3312864/1‑94 (MQ=60)
                                                           CGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTA                                    >  SRR4302533.1622995/1‑94 (MQ=60)
                                                           CGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTA                                    <  SRR4302533.2647300/94‑1 (MQ=60)
                                                           CGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTA                                    <  SRR4302533.5711820/94‑1 (MQ=60)
                                                           CGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTA                                    >  SRR4302533.3783228/1‑94 (MQ=60)
                                                            GCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACNATCCATGAGANAAATANAAACNTNGNTNAC                                   >  SRR4302533.137837/1‑94 (MQ=60)
                                                            GCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTAC                                   >  SRR4302533.1331209/1‑94 (MQ=60)
                                                             CATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACT                                  <  SRR4302533.213437/94‑1 (MQ=60)
                                                              ATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTC                                 >  SRR4302533.5735935/1‑94 (MQ=60)
                                                                TATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAA                               >  SRR4302533.6283547/1‑94 (MQ=60)
                                                                 ATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAA                              >  SRR4302533.3001570/1‑94 (MQ=60)
                                                                  TGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAG                             <  SRR4302533.4448386/94‑1 (MQ=60)
                                                                   GGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGG                            >  SRR4302533.4769634/1‑94 (MQ=60)
                                                                   GGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGG                            <  SRR4302533.682800/94‑1 (MQ=60)
                                                                   GGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGG                            <  SRR4302533.2769676/94‑1 (MQ=60)
                                                                    GTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGT                           >  SRR4302533.2278065/1‑94 (MQ=60)
                                                                    GTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGT                           >  SRR4302533.5453696/1‑94 (MQ=60)
                                                                    GTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGT                           <  SRR4302533.5722008/94‑1 (MQ=60)
                                                                     TTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTG                          >  SRR4302533.5823806/1‑94 (MQ=60)
                                                                      TTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGA                         <  SRR4302533.4484993/94‑1 (MQ=60)
                                                                       TTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGAT                        <  SRR4302533.742720/94‑1 (MQ=60)
                                                                       TTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGAT                        <  SRR4302533.3446027/94‑1 (MQ=60)
                                                                        TGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATG                       >  SRR4302533.276180/1‑94 (MQ=60)
                                                                         GCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGG                      <  SRR4302533.3722719/94‑1 (MQ=60)
                                                                          CCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGC                     >  SRR4302533.1213331/1‑94 (MQ=60)
                                                                            TGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGT                   >  SRR4302533.156132/1‑94 (MQ=60)
                                                                             GTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTT                  >  SRR4302533.3945356/1‑94 (MQ=60)
                                                                             GTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTT                  >  SRR4302533.2654964/1‑94 (MQ=60)
                                                                              TGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTC                 >  SRR4302533.2242109/1‑94 (MQ=60)
                                                                               GGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCT                >  SRR4302533.5057821/1‑94 (MQ=60)
                                                                                 TGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAC              <  SRR4302533.164230/94‑1 (MQ=60)
                                                                                 TATTCAACACTGCACCAGATCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAC              <  SRR4302533.2435471/94‑1 (MQ=60)
                                                                                  GTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACC             <  SRR4302533.5093264/94‑1 (MQ=60)
                                                                                  GTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACC             <  SRR4302533.5963856/94‑1 (MQ=60)
                                                                                   TTCAACACTGCACCAGGTCATTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCA            >  SRR4302533.1704951/1‑94 (MQ=60)
                                                                                   TTCAACACTGCACCAGGTCACTGGAAAGCCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCA            <  SRR4302533.3996095/94‑1 (MQ=60)
                                                                                    TCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCAT           >  SRR4302533.6202120/1‑94 (MQ=60)
                                                                                    TCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCAT           >  SRR4302533.2074531/1‑94 (MQ=60)
                                                                                      AACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCC         <  SRR4302533.1513701/94‑1 (MQ=60)
                                                                                      AACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCC         <  SRR4302533.5786284/94‑1 (MQ=60)
                                                                                       ACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCA        >  SRR4302533.4435999/1‑94 (MQ=60)
                                                                                        CACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAT       >  SRR4302533.2755463/1‑94 (MQ=60)
                                                                                        CACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAT       <  SRR4302533.3990904/94‑1 (MQ=60)
                                                                                        CACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAT       >  SRR4302533.2129926/1‑94 (MQ=60)
                                                                                             CACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAA  <  SRR4302533.979118/94‑1 (MQ=60)
                                                                                             |                                                                                             
TACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAA  >  NC_000913/3873400‑3873586

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