Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,810,321 T→C 17.9% intergenic (+17/+22) coaD → / ← mutM pantetheine‑phosphate adenylyltransferase/formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,810,3210TC17.9% 240.9 / 45.3 117intergenic (+17/+22)coaD/mutMpantetheine‑phosphate adenylyltransferase/formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (41/54);  new base C (10/11);  total (51/66)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-02

Only 100 of 121 total aligned reads displayed.

AAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCG  >  NC_000913/3810222‑3810420
                                                                                                   |                                                                                                   
   agGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATg                                                                                                  <  2:2585627/100‑1 (MQ=255)
   agGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATg                                                                                                  <  2:711828/100‑1 (MQ=255)
       ggCGCGCTATCAGGGTGGTGTCTCCCATTTCCTGCCGGAGAATGTCCATCAGCCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCAt                                                                                              <  1:164586/100‑1 (MQ=255)
       ggCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCAt                                                                                              >  2:2225/1‑100 (MQ=255)
           cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGTCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCgg                                                                                          <  1:2144417/100‑1 (MQ=255)
           cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCgg                                                                                          <  2:593124/100‑1 (MQ=255)
           cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCgg                                                                                          <  1:2497258/100‑1 (MQ=255)
           cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCgg                                                                                          <  1:649277/100‑1 (MQ=255)
           cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCgg                                                                                          <  1:1910008/100‑1 (MQ=255)
           cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCgg                                                                                          <  2:2175004/100‑1 (MQ=255)
           cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCgg                                                                                          <  2:1848611/100‑1 (MQ=255)
           cgcCATCAGGGCGATGTCACCCATTGCCTGCCGGTGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCgg                                                                                          <  2:769844/100‑1 (MQ=255)
            gcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCGgc                                                                                         <  1:684805/100‑1 (MQ=255)
             ccATTAGGGCGATGGCACCCATTTCCTGCCAGAGAATGGCCATCATGTGCTGATGGCGAAGTTAGCTTACCGTTTATTCCGAATGGCAAACCATCCgg                                                                                          <  1:640709/98‑1 (MQ=255)
             ccATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgcg                                                                                        <  2:1308380/100‑1 (MQ=255)
             ccATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCGgca                                                                                        <  2:645824/100‑2 (MQ=255)
              cATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgcgc                                                                                       >  2:857542/1‑100 (MQ=255)
               gtCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgcgcg                                                                                      >  2:1500490/2‑100 (MQ=255)
               aTCAGGGCGATGTCACCCATTTCCTGCCGTAGCATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCAATcagac                                                                                         >  1:2349005/1‑93 (MQ=255)
               aTCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgcgcg                                                                                      <  1:2310439/100‑1 (MQ=255)
               aTCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgcgcg                                                                                      <  2:12815/100‑1 (MQ=255)
               aTCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCGgcat                                                                                       <  1:1927810/99‑3 (MQ=255)
                  aGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCAt                                                                                   <  1:155494/100‑1 (MQ=255)
                    ggCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGa                                                                                 >  2:1166864/1‑100 (MQ=255)
                      cGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGGTGATGGCGAAGTTAGCGTAGCGTTGATGACGGATGGTATGCCATCCGGCGCGCATGaa                                                                                <  1:1552258/99‑1 (MQ=255)
                       gATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAtt                                                                              <  1:838113/100‑1 (MQ=255)
                        aTGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGGTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAtt                                                                              <  1:1785000/99‑1 (MQ=255)
                           tCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACtt                                                                          <  2:620504/100‑1 (MQ=255)
                            catccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTc                                                                         >  2:1282191/4‑100 (MQ=255)
                            cACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTc                                                                         >  1:351292/1‑100 (MQ=255)
                              cccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTg                                                                       >  2:1595590/1‑100 (MQ=255)
                              cccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTg                                                                       >  1:336955/1‑100 (MQ=255)
                              cccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTg                                                                       >  1:1308050/1‑100 (MQ=255)
                              cccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTAATTCTg                                                                       >  2:2435761/1‑100 (MQ=255)
                               ccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTgg                                                                      <  2:443819/100‑1 (MQ=255)
                               ccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTgg                                                                      <  1:1422454/100‑1 (MQ=255)
                                cATTTCCTGCCGGAGAATGTCCATCGGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGc                                                                     <  1:611742/100‑1 (MQ=255)
                                cATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGc                                                                     >  2:1543815/1‑100 (MQ=255)
                                cATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGc                                                                     >  1:1800223/1‑100 (MQ=255)
                                cATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGc                                                                     <  2:1240015/100‑1 (MQ=255)
                                cATTTCCTGCCGGAGAATGTCCATCAGCCGGTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCt                                                                        <  1:1787596/97‑1 (MQ=255)
                                         ccGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGAc                                                            <  1:753559/100‑1 (MQ=255)
                                          cGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACa                                                           >  2:2579945/1‑100 (MQ=255)
                                             agaATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATa                                                        <  2:1612158/100‑1 (MQ=255)
                                                aTGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATaaaa                                                     >  2:704463/1‑100 (MQ=255)
                                                 tGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATaaaa                                                     >  1:2092552/1‑99 (MQ=255)
                                                  gTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACg                                                   >  1:628189/1‑100 (MQ=255)
                                                   tCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGACATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACCATAAAAc                                                    >  1:2416819/1‑98 (MQ=255)
                                                    ccATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGtt                                                 >  1:305355/1‑100 (MQ=255)
                                                    ccATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGtt                                                 >  1:1734979/1‑100 (MQ=255)
                                                    ccATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGtt                                                 >  2:2366310/1‑100 (MQ=255)
                                                     gatctggCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTg                                                <  2:587475/95‑1 (MQ=255)
                                                     cATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCCTGAATTACTTCTGGCACTGCCGACAATAAAACg                                                   >  1:2441954/1‑97 (MQ=255)
                                                     cATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTg                                                <  2:1392788/100‑1 (MQ=255)
                                                     cATCAGGCGCTGATCGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATGCGGCGCGCATGTCTTCCTTCTCGCACTGCCGACAAAAAAACg                                                   >  1:515947/1‑97 (MQ=255)
                                                      aTCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCGGACAATAAAACGTTGc                                               <  2:637455/100‑1 (MQ=255)
                                                          ggCGCTGATGTCGCAGTTAGAGTAGCGTTTATGCCGGATGGTACGCCATCCATCTCGCATGTATTACTTCTGGCCCTGCCGACAATAAAAc                                                    <  1:2097790/91‑1 (MQ=255)
                                                          ggCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTgagatc                                           >  1:587473/1‑95 (MQ=255)
                                                           gcgcTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTACGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTAccc                                             >  1:575879/1‑97 (MQ=255)
                                                               tGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTgcgc                                      >  2:34513/1‑100 (MQ=255)
                                                               tGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTgcgc                                      <  2:2238602/100‑1 (MQ=255)
                                                                gATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCa                                     <  2:556938/100‑1 (MQ=255)
                                                                gATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCa                                     >  1:1337116/1‑100 (MQ=255)
                                                                 aTGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCa                                     >  1:776970/1‑99 (MQ=255)
                                                                  tGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATg                                   >  2:465758/1‑100 (MQ=255)
                                                                    gCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGtt                                 <  1:1967436/100‑1 (MQ=255)
                                                                    gCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGtt                                 >  2:484839/1‑100 (MQ=255)
                                                                    gCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGtt                                 >  2:2375398/1‑100 (MQ=255)
                                                                     cGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGttt                                <  1:2579926/100‑1 (MQ=255)
                                                                     cGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGttt                                <  1:614693/100‑1 (MQ=255)
                                                                     cGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGttt                                <  1:2480756/100‑1 (MQ=255)
                                                                       aaGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAg                              <  1:760037/100‑1 (MQ=255)
                                                                           tagcgtagcgTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGc                          <  2:57384/100‑1 (MQ=255)
                                                                            agcgtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGcc                         >  2:2494622/1‑100 (MQ=255)
                                                                             gcgtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCca                        >  2:2525748/1‑100 (MQ=255)
                                                                              cgtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCcac                       <  1:644502/100‑1 (MQ=255)
                                                                              cgtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCcac                       <  1:2294977/100‑1 (MQ=255)
                                                                                    atgcgcgCCGGATGGCATACCATCCGGCGCGCATGAATGACTACTGGCACTGCCGACAATAAAACGTTGCTCGCTGCGCCTGTTTAGACGCCACAATc                   >  1:1189499/7‑98 (MQ=25)
                                                                                      ttATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCg               <  1:988419/100‑1 (MQ=255)
                                                                                       tATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATTTTTAGTCGCCACAATCGGCg               <  1:422893/99‑1 (MQ=255)
                                                                                       cgcgCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGTCCGCTGCGCATGTTTAGTCGCCACAATCGGCGt              >  1:1686566/4‑100 (MQ=255)
                                                                                       cgcgCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGt              >  2:837550/4‑100 (MQ=255)
                                                                                         tGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAc            >  2:953734/1‑100 (MQ=255)
                                                                                          gCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAcc           <  1:2125763/100‑1 (MQ=255)
                                                                                          gCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAcc           >  1:166225/1‑100 (MQ=255)
                                                                                          gCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAcc           >  2:381916/1‑100 (MQ=255)
                                                                                           ccGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCg          <  2:284889/100‑1 (MQ=255)
                                                                                           ccGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCg          >  1:2443558/1‑100 (MQ=255)
                                                                                           ccGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCg          >  2:2226392/1‑100 (MQ=255)
                                                                                           ccGAATGGTATGCCATCCGGCGCGCATGAATTACTTCTGACACTGCCGACAAAAAAACGTTACCCGCTGCCCACGTTTCGTCCCcaca                      >  1:1542912/1‑88 (MQ=255)
                                                                                            cGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGc         <  2:657898/100‑1 (MQ=255)
                                                                                             ggATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGca        <  1:992873/100‑1 (MQ=255)
                                                                                             ggATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGca        <  2:2162385/100‑1 (MQ=255)
                                                                                             ggATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGca        >  1:241861/1‑100 (MQ=255)
                                                                                              gATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGcac       <  1:7202/100‑1 (MQ=255)
                                                                                                  gTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACAccc   <  1:1547234/100‑1 (MQ=255)
                                                                                                  gTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACAccc   <  1:1345639/100‑1 (MQ=255)
                                                                                                   tATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCg  <  2:2438585/100‑1 (MQ=255)
                                                                                                   tATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCg  <  1:1557259/100‑1 (MQ=255)
                                                                                                   tATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCg  >  2:1203544/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCG  >  NC_000913/3810222‑3810420

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