Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,810,324 G→A 17.5% intergenic (+20/+19) 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,3240GA17.5% 271.8 / 44.5 120intergenic (+20/+19)coaD/mutMpantetheine‑phosphate adenylyltransferase/formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (44/55);  new base A (10/11);  total (54/66)
Fisher's exact test for biased strand distribution p-value = 8.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

Only 100 of 124 total aligned reads displayed.

AGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACA  >  NC_000913/3810225‑3810423
                                                                                                   |                                                                                                   
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                                                                                             <  2:2175004/100‑1 (MQ=255)
        cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCgg                                                                                             <  1:1910008/100‑1 (MQ=255)
        cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCgg                                                                                             <  1:2497258/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)
           cATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgcgc                                                                                          <  2:2000117/100‑1 (MQ=255)
            gtCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgcgcg                                                                                         >  2:1500490/2‑100 (MQ=255)
            aTCAGGGCGATGTCACCCATTTCCTGCCGTAGCATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCAATcagac                                                                                            >  1:2349005/1‑93 (MQ=255)
            aTCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgcgcg                                                                                         <  2:12815/100‑1 (MQ=255)
            aTCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgcgcg                                                                                         <  1:2310439/100‑1 (MQ=255)
            aTCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCGgcat                                                                                          <  1:1927810/99‑3 (MQ=255)
                 ggCGATGTCTCCCATTGCCTGCCGGATAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCGgcataaacgc                                                                                    <  2:2320382/100‑9 (MQ=37)
                   cGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGGTGATGGCGAAGTTAGCGTAGCGTTGATGACGGATGGTATGCCATCCGGCGCGCATGaa                                                                                   <  1:1552258/99‑1 (MQ=255)
                   cGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAt                                                                                  >  1:684850/1‑100 (MQ=255)
                    gATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAtt                                                                                 <  1:838113/100‑1 (MQ=255)
                     aTGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGGTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAtt                                                                                 <  1:1785000/99‑1 (MQ=255)
                        tCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACtt                                                                             <  2:620504/100‑1 (MQ=255)
                         cACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTc                                                                            >  1:351292/1‑100 (MQ=255)
                           cccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTg                                                                          >  2:1595590/1‑100 (MQ=255)
                           cccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTg                                                                          >  1:1308050/1‑100 (MQ=255)
                           cccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTg                                                                          >  1:336955/1‑100 (MQ=255)
                             cATTTCCTGCCGGAGAATGTCCATCGGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGc                                                                        <  1:611742/100‑1 (MQ=255)
                             cATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGc                                                                        >  1:1800223/1‑100 (MQ=255)
                             cATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGc                                                                        <  2:628870/100‑1 (MQ=255)
                             cATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGc                                                                        <  2:1240015/100‑1 (MQ=255)
                                       cGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACa                                                              >  2:2579945/1‑100 (MQ=255)
                                       cGGAGAACGTCCATCAGGCGCAGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACa                                                              <  1:616643/100‑1 (MQ=255)
                                          agaATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATa                                                           <  2:1612158/100‑1 (MQ=255)
                                             aTGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATaaaa                                                        >  2:704463/1‑100 (MQ=255)
                                              tGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATaaaa                                                        >  1:2092552/1‑99 (MQ=255)
                                              tGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAAc                                                       <  2:1249613/100‑1 (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)
                                                       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                                        >  1:1337116/1‑100 (MQ=255)
                                                             gATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCa                                        <  2:556938/100‑1 (MQ=255)
                                                              aTGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCa                                        >  1:776970/1‑99 (MQ=255)
                                                               tGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATg                                      >  2:465758/1‑100 (MQ=255)
                                                                 gCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGtt                                    >  1:1643933/1‑100 (MQ=255)
                                                                 gCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGtt                                    >  2:484839/1‑100 (MQ=255)
                                                                 gCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGtt                                    >  2:2375398/1‑100 (MQ=255)
                                                                 gCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGtt                                    <  1:1967436/100‑1 (MQ=255)
                                                                  cGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGttt                                   <  1:2480756/100‑1 (MQ=255)
                                                                    aaGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAg                                 <  1:760037/100‑1 (MQ=255)
                                                                         agcgtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGcc                            >  2:2494622/1‑100 (MQ=255)
                                                                          gcgtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGcc                            <  1:1786842/99‑1 (MQ=255)
                                                                          gcgtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCca                           >  2:2525748/1‑100 (MQ=255)
                                                                          gcgtagcgtTTATGCCGGATGGAATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGcc                            <  1:540678/99‑1 (MQ=255)
                                                                           cgtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCcac                          <  1:644502/100‑1 (MQ=255)
                                                                           cgtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCcac                          <  1:2294977/100‑1 (MQ=255)
                                                                                cgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCg                     >  2:505005/1‑100 (MQ=255)
                                                                                 atgcgcgCCGGATGGCATACCATCCGGCGCGCATGAATGACTACTGGCACTGCCGACAATAAAACGTTGCTCGCTGCGCCTGTTTAGACGCCACAATc                      >  1:1189499/7‑98 (MQ=25)
                                                                                  tttATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGc                   <  1:1471773/100‑1 (MQ=255)
                                                                                  tttATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGc                   <  2:2497363/100‑1 (MQ=255)
                                                                                   ttATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCg                  <  2:1435395/100‑1 (MQ=255)
                                                                                   ttATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCg                  <  1:988419/100‑1 (MQ=255)
                                                                                    tATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATTTTTAGTCGCCACAATCGGCg                  <  1:422893/99‑1 (MQ=255)
                                                                                    cgcgCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGt                 >  2:837550/4‑100 (MQ=255)
                                                                                      tGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAc               >  2:953734/1‑100 (MQ=255)
                                                                                       gCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAcc              <  1:2125763/100‑1 (MQ=255)
                                                                                       gCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAcc              >  2:338211/1‑100 (MQ=255)
                                                                                       gCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAcc              >  1:166225/1‑100 (MQ=255)
                                                                                       gCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAcc              >  2:381916/1‑100 (MQ=255)
                                                                                        ccGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCg             <  2:284889/100‑1 (MQ=255)
                                                                                        ccGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCg             >  2:2226392/1‑100 (MQ=255)
                                                                                        ccGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCg             >  1:2443558/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)
                                                                                          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:1203544/1‑100 (MQ=255)
                                                                                                tATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCg     <  1:1557259/100‑1 (MQ=255)
                                                                                                tATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCg     <  2:2438585/100‑1 (MQ=255)
                                                                                                 aTGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGa    <  1:659665/100‑1 (MQ=255)
                                                                                                  tGCCATCCGGCGCGCATGACTTACTTCTGGCACTGCCGCCAATAAAACGTTGCCCGCTGCGCCTGTTTAGTCTCCCCAATCGGCGTACCGCACACCCGac   >  1:1165341/1‑100 (MQ=255)
                                                                                                  tGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGac   >  1:1746748/1‑100 (MQ=255)
                                                                                                   gCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGaca  <  2:2387251/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACA  >  NC_000913/3810225‑3810423

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: