Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,891,179 G→C S240R (AGC→AGG cysJ ← sulfite reductase, alpha subunit, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,891,1790GC99.2% 462.6 / ‑4.0 130S240R (AGC→AGGcysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (56/73);  total (56/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 130 total aligned reads displayed.

GGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCG  >  NC_000913/2891084‑2891273
                                                                                               |                                                                                              
ggCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAg                                                                                            <  2:2535225/100‑1 (MQ=255)
 gCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGc                                                                                           <  2:1994811/100‑1 (MQ=255)
  cTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGcc                                                                                          <  1:1491941/100‑1 (MQ=255)
    ggTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGccac                                                                                        >  2:643093/1‑100 (MQ=255)
     gTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGccacc                                                                                       >  1:1461217/1‑100 (MQ=255)
       aaCGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAg                                                                                     >  1:58676/1‑100 (MQ=255)
        aCGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAgc                                                                                    <  2:599238/100‑1 (MQ=255)
         cGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAgcg                                                                                   >  1:324855/1‑100 (MQ=255)
          gCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAgcgg                                                                                  <  1:1851967/100‑1 (MQ=255)
          gCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAgcgg                                                                                  <  2:2420146/100‑1 (MQ=255)
           cATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAgcggc                                                                                 >  2:1821206/1‑100 (MQ=255)
           cATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAgcggc                                                                                 >  2:2361295/1‑100 (MQ=255)
            aTGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAgcggcg                                                                                >  2:881727/1‑100 (MQ=255)
             tGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGgcgc                                                                               >  1:2366503/1‑100 (MQ=255)
              gCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTTATTTTTTCGTCAGGTAGTTTGCTGGTTAACAGATAGCCTAGCCACCAgcgtttcg                                                                              >  2:2615023/1‑95 (MQ=255)
              gCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAACCACCAGCGGcgcg                                                                              >  2:2103404/1‑100 (MQ=255)
               cccGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGt                                                                             <  1:2194923/100‑1 (MQ=255)
                ccGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGcgcgcc                                                                            >  1:731106/1‑98 (MQ=255)
                  gAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCtt                                                                          <  1:1683125/100‑1 (MQ=255)
                      cACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCt                                                                      >  2:554265/1‑100 (MQ=255)
                      cACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCt                                                                      >  1:2410917/1‑100 (MQ=255)
                       aCCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGc                                                                       <  2:1867608/98‑1 (MQ=255)
                       aCCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTg                                                                     >  1:1830689/1‑100 (MQ=255)
                        ccTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGt                                                                    >  2:2189774/1‑100 (MQ=255)
                        ccTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGt                                                                    >  1:1174307/1‑100 (MQ=255)
                           aaGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACg                                                                 <  2:2834864/100‑1 (MQ=255)
                            aGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACgg                                                                >  1:1004086/1‑100 (MQ=255)
                             gTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACggg                                                               >  1:1002126/1‑100 (MQ=255)
                             gTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACggg                                                               >  1:1297987/1‑100 (MQ=255)
                             gTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACggg                                                               >  1:1039945/1‑100 (MQ=255)
                              tCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGc                                                              <  1:3057992/100‑1 (MQ=255)
                               cAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCt                                                             >  1:3053342/1‑100 (MQ=255)
                               cAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCt                                                             >  2:2412465/1‑100 (MQ=255)
                                aagttCGATATAGCGAACGTCTTTTTCAGAGTTAAGCCCGGTAATTTTCTGGTTAAGAGAGAGCATAGCCACCAGCGGCGCGAGTTTGCTATACGGACTg                                                            >  2:2015443/4‑100 (MQ=255)
                                 aTTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTgg                                                           <  2:1461203/100‑1 (MQ=255)
                                 aTTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTgg                                                           >  1:1166662/1‑100 (MQ=255)
                                  tttCGATATGGCGAACGTCTTTTTTTTTTTTTTTCCCGGTAATTTTCTGGTTAACAGAGAGCTTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGgt                                                          >  2:1317820/1‑100 (MQ=255)
                                  tttCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGgt                                                          <  2:720471/100‑1 (MQ=255)
                                   ttCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGgtg                                                         <  1:1584509/100‑1 (MQ=255)
                                   ttCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGgtg                                                         <  1:1127946/100‑1 (MQ=255)
                                     ccgatctgCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGtgtg                                                       <  2:2240487/93‑1 (MQ=255)
                                     cGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGtgtg                                                       <  2:2148160/100‑1 (MQ=255)
                                       agatGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGa                                                     <  2:113364/98‑1 (MQ=255)
                                        tatGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGAt                                                    <  1:1739638/100‑1 (MQ=255)
                                         atGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGAtt                                                   <  2:1119067/100‑1 (MQ=255)
                                         atGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGAtt                                                   >  2:842649/1‑100 (MQ=255)
                                          tGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGAttt                                                  >  2:2105642/1‑100 (MQ=255)
                                          tGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGAttt                                                  >  2:2661176/1‑100 (MQ=255)
                                           ggCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTc                                                 >  1:481458/1‑100 (MQ=255)
                                           ggCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTc                                                 >  1:366795/1‑100 (MQ=255)
                                            gCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGtgtgagatcgg                                                >  1:2240510/1‑93 (MQ=255)
                                            gCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCa                                                <  2:2636272/100‑1 (MQ=255)
                                            gCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGAGGATTTCa                                                >  1:270465/1‑100 (MQ=255)
                                             cGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCAt                                               <  1:2054287/100‑1 (MQ=255)
                                             cGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCAt                                               <  2:324853/100‑1 (MQ=255)
                                              gAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCAtt                                              <  2:652431/100‑1 (MQ=255)
                                               aaCGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTc                                                 <  2:1930201/96‑1 (MQ=255)
                                               aaCGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCAttt                                             >  2:2184523/1‑100 (MQ=255)
                                                 cGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTAc                                           <  1:2274547/100‑1 (MQ=255)
                                                  gTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTAcc                                          >  2:1101507/1‑100 (MQ=255)
                                                   tCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACcg                                         <  1:66191/100‑1 (MQ=255)
                                                    cTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACcgc                                        >  2:2362566/1‑100 (MQ=255)
                                                       tttCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGcc                                     <  2:58676/100‑1 (MQ=255)
                                                       tttCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGcc                                     <  1:2799821/100‑1 (MQ=255)
                                                       tttCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGcc                                     <  1:3055153/100‑1 (MQ=255)
                                                         tCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAg                                   <  1:2445644/100‑1 (MQ=255)
                                                          cAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGt                                  <  2:1933818/100‑1 (MQ=255)
                                                          cAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGt                                  <  1:1327750/100‑1 (MQ=255)
                                                             agTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGc                               <  2:2766372/100‑1 (MQ=255)
                                                              gTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCg                              <  2:1756404/100‑1 (MQ=255)
                                                              gTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCg                              <  2:3036747/100‑1 (MQ=255)
                                                              gTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCg                              >  1:986319/1‑100 (MQ=255)
                                                                ccggatgggTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGAc                            <  1:2502962/93‑1 (MQ=255)
                                                                 aCGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACg                           <  1:2302252/100‑1 (MQ=255)
                                                                 aCGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACg                           <  1:2361318/100‑1 (MQ=255)
                                                                  cGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACgg                          >  1:1073674/1‑100 (MQ=255)
                                                                    cccGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGAt                        <  1:1642033/100‑1 (MQ=255)
                                                                     ccGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGAtt                       >  2:1549927/1‑100 (MQ=255)
                                                                      cGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTg                      <  1:1354942/100‑1 (MQ=255)
                                                                       ggTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGc                     <  1:1409261/100‑1 (MQ=255)
                                                                       ggTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGc                     >  1:397154/1‑100 (MQ=255)
                                                                       ggTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGc                     >  1:395820/1‑100 (MQ=255)
                                                                        gTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCg                    >  2:231787/1‑100 (MQ=255)
                                                                        gTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCg                    <  2:3012288/100‑1 (MQ=255)
                                                                           aTTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAg                 <  2:1617075/100‑1 (MQ=255)
                                                                           aTTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAg                 >  2:1700122/1‑100 (MQ=255)
                                                                           aTTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAg                 <  1:753275/100‑1 (MQ=255)
                                                                           aTTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCAACGGATTGCGaa                  >  2:802259/1‑99 (MQ=255)
                                                                            ttttCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAgg                <  2:2923115/100‑1 (MQ=255)
                                                                            ttttCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAgg                >  2:3052475/1‑100 (MQ=255)
                                                                             tttCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGt               >  1:2150480/1‑100 (MQ=255)
                                                                              ttCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTg              <  1:1682434/100‑1 (MQ=255)
                                                                                  ngTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgc          >  2:2974159/2‑100 (MQ=255)
                                                                                  ggTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgc          >  1:592505/1‑100 (MQ=255)
                                                                                  ggTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgc          <  1:1490362/100‑1 (MQ=255)
                                                                                   gTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTCGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgcg         <  2:1711920/100‑1 (MQ=255)
                                                                                     tAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGAc       <  1:2189797/100‑1 (MQ=255)
                                                                                      aaCAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACa      <  1:1837241/100‑1 (MQ=255)
                                                                                         agagagCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGgc   >  2:1786636/1‑100 (MQ=255)
                                                                                          gagagCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGgcg  >  1:1565695/1‑100 (MQ=255)
                                                                                               |                                                                                              
GGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCG  >  NC_000913/2891084‑2891273

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