Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 740,004 G→A 84.3% M166I (ATG→ATA phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913740,0040GA84.3% 246.8 / 36.0 107M166I (ATG→ATAphrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (10/7);  new base A (46/44);  total (56/51)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

Only 100 of 107 total aligned reads displayed.

CCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATAC  >  NC_000913/739908‑740099
                                                                                                |                                                                                               
ccTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCg                                                                                              >  2:1258815/1‑100 (MQ=255)
 cTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCgg                                                                                             >  1:1230536/1‑100 (MQ=255)
 cTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCgg                                                                                             <  1:2399331/100‑1 (MQ=255)
    ccgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGt                                                                                          <  1:2800611/100‑1 (MQ=255)
    ccgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGt                                                                                          <  2:1756231/100‑1 (MQ=255)
    ccgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGt                                                                                          <  2:1756221/100‑1 (MQ=255)
     cgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTg                                                                                         <  2:1117117/100‑1 (MQ=255)
      gccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGc                                                                                        >  2:1818697/1‑100 (MQ=255)
      gccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGc                                                                                        >  2:1538886/1‑100 (MQ=255)
        cTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGt                                                                                      <  2:1699750/100‑1 (MQ=255)
          ggCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCg                                                                                    >  1:1264286/1‑100 (MQ=255)
          ggCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCg                                                                                    >  1:204113/1‑100 (MQ=255)
           gcgcGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGc                                                                                   <  1:1593685/100‑1 (MQ=255)
           gcgcGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGc                                                                                   <  2:434499/100‑1 (MQ=255)
            cgcgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCt                                                                                  >  1:493453/1‑100 (MQ=255)
               ggTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTgcg                                                                               <  1:2548275/100‑1 (MQ=255)
               ggTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTgcg                                                                               <  2:2646364/100‑1 (MQ=255)
                 tgatgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGcc                                                                             <  2:110315/100‑1 (MQ=255)
                  gatgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGcc                                                                             >  1:1343510/1‑99 (MQ=255)
                    tgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCaaa                                                                          <  2:2474808/100‑1 (MQ=255)
                     gaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTg                                                                                 >  1:2754386/1‑92 (MQ=255)
                     gaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCaaaa                                                                         <  1:1952889/100‑1 (MQ=255)
                     gaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCGTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCaaaa                                                                         <  2:973935/100‑1 (MQ=255)
                      aCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAg                                                                        <  1:2213526/100‑1 (MQ=255)
                      aCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAg                                                                        >  2:2839389/1‑100 (MQ=255)
                         ggTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTc                                                                     >  2:764392/1‑100 (MQ=255)
                         ggTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTc                                                                     >  2:2686206/1‑100 (MQ=255)
                         ggTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTc                                                                     >  2:2686670/1‑100 (MQ=255)
                          gTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCg                                                                    >  2:1800387/1‑100 (MQ=255)
                             aTCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCgtag                                                                 >  2:1854858/1‑100 (MQ=255)
                              tCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgtagt                                                                >  1:1563332/1‑100 (MQ=255)
                               cACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgtagta                                                               >  2:1229367/1‑100 (MQ=255)
                               cACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgtagta                                                               <  2:614943/100‑1 (MQ=255)
                               cACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGACTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgtagta                                                               >  1:1560001/1‑100 (MQ=255)
                                aCGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGtagtag                                                              >  2:13552/1‑100 (MQ=255)
                                aCGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGtagtag                                                              <  1:47535/100‑1 (MQ=255)
                                 cGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATAACGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGc                                                             >  1:2265645/1‑100 (MQ=255)
                                   agaTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCgg                                                           >  1:479247/1‑100 (MQ=255)
                                    gaTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGa                                                          >  2:351075/1‑100 (MQ=255)
                                    gaTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGa                                                          >  2:749400/1‑100 (MQ=255)
                                     aTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGAt                                                         <  2:1989955/100‑1 (MQ=255)
                                      tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                        <  2:260364/100‑1 (MQ=255)
                                      tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                        <  2:232738/100‑1 (MQ=255)
                                      tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                        <  2:2037889/100‑1 (MQ=255)
                                       gTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCg                                                       <  1:642904/100‑1 (MQ=255)
                                        tACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTTGCTGCGCCAAAAGTTCGTAGTAGCgg                                                           >  1:782949/1‑95 (MQ=255)
                                         aCAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGCAGGGATGCCGGAGTGCGTCGCTGCGCCCAAAGCTAGTAGTAGCGGATCGAt                                                     >  1:2808933/1‑100 (MQ=255)
                                         aCAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGAt                                                     >  1:1474845/1‑100 (MQ=255)
                                         aCAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGAt                                                     >  2:2027549/1‑100 (MQ=255)
                                          cAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATa                                                    <  1:48792/100‑1 (MQ=255)
                                          cAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATa                                                    >  2:23614/1‑100 (MQ=255)
                                          cAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATa                                                    <  2:2286018/100‑1 (MQ=255)
                                             aGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAATGCGTCGCTGCGCCAAAAGTTCGTAGTACCGGATCGATagtg                                                 <  2:1968840/100‑3 (MQ=255)
                                              gTCTTGACGCCTTTTAAGAATGCCGGGCTGAAACGGCTGCGGGAAGGGAGGCCGGAGTGCGTCGCTGCGCCAACAGTTCGTAGGAGCGGATCGATag                                                   <  1:1988859/97‑1 (MQ=255)
                                               tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGcc                                               >  1:333100/1‑100 (MQ=255)
                                                cTTTACGCCTTTTAAGAATGCCTTGCTGAAACGGCTGCGGGAATGGATACCGGAGTGCGTCGCTGCGCCACAAGTTAGTAGTAGCGGATCGATaga                                                  >  1:542454/1‑96 (MQ=255)
                                                cTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGccc                                              >  2:1752931/1‑100 (MQ=255)
                                                 tttACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCt                                             <  1:13552/100‑1 (MQ=255)
                                                  ttACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTc                                            <  1:2775685/100‑1 (MQ=255)
                                                    aCGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGc                                          >  1:1190391/1‑100 (MQ=255)
                                                    aCGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGc                                          >  1:2218438/1‑100 (MQ=255)
                                                     cGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGcc                                         <  1:509353/100‑1 (MQ=255)
                                                     cGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGcc                                         >  2:2215752/1‑100 (MQ=255)
                                                     cGCCCTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGcc                                         <  2:1408034/100‑1 (MQ=255)
                                                       ccTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccat                                       >  2:1656347/1‑100 (MQ=255)
                                                       ccTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccat                                       >  1:49559/1‑100 (MQ=255)
                                                        cTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatc                                      >  1:2318270/1‑100 (MQ=255)
                                                        cTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatc                                      >  2:2443977/1‑100 (MQ=255)
                                                        cTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatc                                      <  2:1177471/100‑1 (MQ=255)
                                                         ttttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcc                                     <  2:824122/100‑1 (MQ=255)
                                                         ttttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcc                                     >  2:36450/1‑100 (MQ=255)
                                                         ttttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcc                                     >  1:268619/1‑100 (MQ=255)
                                                          tttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcca                                    <  1:784126/100‑1 (MQ=255)
                                                          tttAAGAATGCCTGGCTGAAAAGGCTGCGGGAAGGGATGACGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcca                                    >  1:1344700/1‑100 (MQ=255)
                                                           ttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatacat                                   >  2:675487/1‑100 (MQ=255)
                                                            tAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCAtt                                  <  2:2861836/100‑1 (MQ=255)
                                                             aaGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGCCGCTGCGCCAAAAGTTAGTAGTAGCGGCTCGACAGAGCCCTCGCCATCCACTa                                 >  1:790637/1‑100 (MQ=255)
                                                              aGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTAc                                >  2:27355/1‑100 (MQ=255)
                                                              aGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTAc                                >  1:257311/1‑100 (MQ=255)
                                                                  tGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCACTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                            <  1:1867625/100‑1 (MQ=255)
                                                                    ccTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGaa                          >  2:748716/1‑100 (MQ=255)
                                                                      tGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAAtt                        >  2:342425/1‑100 (MQ=255)
                                                                         cTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATc                     >  2:922280/1‑100 (MQ=255)
                                                                           gAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCt                   >  2:2597962/1‑100 (MQ=255)
                                                                             aaCGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGAGAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCg                 <  2:2690027/100‑1 (MQ=255)
                                                                               cGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTc               <  2:1282889/100‑1 (MQ=255)
                                                                                  cTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt            >  1:981776/1‑100 (MQ=255)
                                                                                   tGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTc           <  1:1650416/100‑1 (MQ=255)
                                                                                     cGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTc           >  2:2019661/1‑98 (MQ=255)
                                                                                     cGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCtt         >  2:2304528/1‑100 (MQ=255)
                                                                                      gggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCttt        <  2:2754404/100‑1 (MQ=255)
                                                                                       ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc       <  1:2894961/100‑1 (MQ=255)
                                                                                       ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc       <  1:2858233/100‑1 (MQ=255)
                                                                                       ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc       <  2:1461615/100‑1 (MQ=255)
                                                                                       ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc       <  1:158843/100‑1 (MQ=255)
                                                                                         aaGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCgn     <  2:1560922/100‑2 (MQ=255)
                                                                                         aaGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGa     >  2:1912209/1‑100 (MQ=255)
                                                                                         aaGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGa     <  1:1450398/100‑1 (MQ=255)
                                                                                           gggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATa   >  2:1570151/1‑100 (MQ=255)
                                                                                            ggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTATCGATAc  <  2:1533097/100‑1 (MQ=255)
                                                                                                |                                                                                               
CCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATAC  >  NC_000913/739908‑740099

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: