Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,952 T→C 16.1% intergenic (‑39/+34) ridA ← / ← pyrI enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,470,9520TC16.1% 285.7 / 45.5 124intergenic (‑39/+34)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base T (48/56);  new base C (14/6);  total (62/62)
Fisher's exact test for biased strand distribution p-value = 8.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 126 total aligned reads displayed.

GCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATC  >  NC_000913/4470854‑4471051
                                                                                                  |                                                                                                   
gCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAtt                                                                                                    >  2:743973/1‑100 (MQ=255)
   cTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATtata                                                                                                 <  2:1227404/100‑1 (MQ=255)
        cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAg                                                                                            >  2:1486903/1‑100 (MQ=255)
        cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAg                                                                                            >  2:1444023/1‑100 (MQ=255)
        cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAg                                                                                            <  2:1045883/100‑1 (MQ=255)
          tAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGcc                                                                                          >  2:1938527/1‑100 (MQ=255)
           aaGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCa                                                                                         <  2:2040975/100‑1 (MQ=255)
           aaGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCa                                                                                         <  1:2253221/100‑1 (MQ=255)
            aGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAg                                                                                        <  2:2770155/100‑1 (MQ=255)
             ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGa                                                                                       >  2:1689936/1‑100 (MQ=255)
             ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGACAGa                                                                                       >  2:1876115/1‑100 (MQ=255)
              gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                      <  1:2891310/100‑1 (MQ=38)
              gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                      <  1:450243/100‑1 (MQ=38)
              gACCGATAGCTGCCGGTGCATTTTCCATCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                      >  2:1774518/1‑100 (MQ=255)
               aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                     >  2:1102803/1‑100 (MQ=255)
               aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTAGCTCATGATTTCTCCTTTAATACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                     <  1:1583713/100‑1 (MQ=37)
                ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                    <  1:1382138/100‑1 (MQ=38)
                   atggCTGCCGGTGCATTTTACGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                 <  1:34720/97‑1 (MQ=255)
                    tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                                <  2:2181913/100‑1 (MQ=255)
                       cTGCCGGTGCATTTTCCGTCGCGATTGTTTTGCGCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                              <  1:41803/99‑1 (MQ=255)
                       cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAc                                                                             <  1:2064992/100‑1 (MQ=255)
                        tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                            >  1:308636/1‑100 (MQ=255)
                        tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                            <  1:2647731/100‑1 (MQ=255)
                           cGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAc                                                                         <  2:1202648/100‑1 (MQ=255)
                            ggTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                        >  2:347923/1‑100 (MQ=255)
                             gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                       <  2:1510927/100‑1 (MQ=255)
                             gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                       <  2:365411/100‑1 (MQ=255)
                              tGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAccac                                                                      <  1:1119862/100‑3 (MQ=255)
                              tGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGc                                                                      >  1:592993/1‑100 (MQ=255)
                               gCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCa                                                                     >  2:2297395/1‑100 (MQ=255)
                                cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaa                                                                    >  1:1757041/1‑100 (MQ=255)
                                  ttttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACTGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatt                                                                  <  1:1862475/100‑1 (MQ=255)
                                  ttttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatt                                                                  >  2:645221/1‑100 (MQ=255)
                                   tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatta                                                                 <  2:821539/100‑1 (MQ=255)
                                       cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGGCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                               >  2:2760685/1‑98 (MQ=255)
                                           gcgATTGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTg                                                             <  1:743971/96‑1 (MQ=255)
                                           gcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCa                                                         <  2:2402554/100‑1 (MQ=255)
                                            cgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAACTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTCATTTGc                                                           >  1:2432506/1‑97 (MQ=255)
                                               tAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                     <  2:1031187/100‑1 (MQ=255)
                                                aGTTTTGCTCATGCTTCCTCCTTTACTACAGCGGTAAAAAGTCTGCCTCATTATAGGGAGCCAGACTTATATTACCACCCGCAATTAATCCGCCAGcccc                                                    <  1:2318959/100‑1 (MQ=255)
                                                aGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacc                                                    <  2:745342/100‑1 (MQ=255)
                                                 gTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacca                                                   <  2:2085864/100‑1 (MQ=255)
                                                 gTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacca                                                   >  2:833695/1‑100 (MQ=255)
                                                   tttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACcaca                                                 <  2:1350349/100‑1 (MQ=255)
                                                   tttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACcaca                                                 <  2:2842560/100‑1 (MQ=255)
                                                    ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                                <  1:2435550/100‑1 (MQ=255)
                                                      gCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTa                                              >  2:1407626/1‑100 (MQ=255)
                                                       cTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                             >  2:1889980/1‑100 (MQ=255)
                                                        tCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                            >  2:1635901/1‑100 (MQ=255)
                                                         cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTACTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                           <  1:1552685/100‑1 (MQ=255)
                                                          aTGATTTCTCCTTGATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATggg                                          <  1:566425/100‑1 (MQ=255)
                                                             aTTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                       <  2:58306/100‑1 (MQ=255)
                                                              tttCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaaa                                      >  2:60952/1‑100 (MQ=255)
                                                              tttCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaaa                                      >  2:801694/1‑100 (MQ=255)
                                                                tctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAct                                    <  2:2242792/100‑1 (MQ=255)
                                                                 ctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                   >  2:760890/1‑100 (MQ=255)
                                                                 ctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTCGCCAGCACCACATTATGGGAAAActc                                   >  1:336176/1‑100 (MQ=255)
                                                                  tcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActct                                  >  1:2902235/1‑100 (MQ=255)
                                                                   ccTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                 >  2:464797/1‑100 (MQ=255)
                                                                   ccTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                 >  2:230508/1‑100 (MQ=255)
                                                                    cTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                >  1:1309391/1‑100 (MQ=255)
                                                                    cTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                >  1:2240026/1‑100 (MQ=255)
                                                                    cTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                <  2:497570/100‑1 (MQ=255)
                                                                      ttattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttttt                              <  2:1983072/100‑1 (MQ=255)
                                                                       tattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                             >  1:2058452/1‑100 (MQ=255)
                                                                        attaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                            >  2:2527228/1‑100 (MQ=255)
                                                                         ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                           >  1:2863885/1‑100 (MQ=255)
                                                                         ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                           <  1:1641421/100‑1 (MQ=255)
                                                                           gttggtaataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                         >  1:2850264/10‑100 (MQ=21)
                                                                           aCAGCGGTAAAAAGTCTGGCTCATTATAGGTAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                           <  1:2666521/98‑1 (MQ=255)
                                                                           aCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAACTGGCCAGCACCACATTATGGGAAATCTCTTTTTCACAg                         <  1:1798008/100‑1 (MQ=255)
                                                                            cAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGt                        >  1:2293029/1‑100 (MQ=255)
                                                                                ggTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACAACATTATGGGAAAACTCTTTTTCACAGTAttt                    >  1:1705096/1‑100 (MQ=255)
                                                                                 gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                   <  2:1309393/100‑1 (MQ=255)
                                                                                 gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                   >  1:2202323/1‑100 (MQ=255)
                                                                                 gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                   <  2:2240036/100‑1 (MQ=255)
                                                                                  tAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                  <  2:2810714/100‑1 (MQ=255)
                                                                                  ataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                  >  2:1566626/3‑100 (MQ=37)
                                                                                   taaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCa                 >  2:1884069/2‑100 (MQ=37)
                                                                                   taaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCa                 >  2:1075931/2‑100 (MQ=37)
                                                                                     aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt               >  1:2018534/1‑100 (MQ=37)
                                                                                     aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGACCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt               >  1:2172910/1‑100 (MQ=37)
                                                                                      aaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt              >  2:2179948/1‑100 (MQ=38)
                                                                                       aGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg             >  1:182233/1‑100 (MQ=38)
                                                                                       aGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg             <  2:1723345/100‑1 (MQ=255)
                                                                                       aGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg             <  2:1480584/100‑1 (MQ=255)
                                                                                        gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg             >  1:83751/1‑99 (MQ=38)
                                                                                        gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGa            >  2:2560029/1‑100 (MQ=38)
                                                                                         tCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg           >  1:899361/1‑100 (MQ=38)
                                                                                         tCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg           <  2:2430189/100‑1 (MQ=255)
                                                                                         tCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg           <  2:2364879/100‑1 (MQ=255)
                                                                                          cTGGCTCCCTATAATGAGCCAGACTCTTATTACCAACCGCAAATAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc          >  1:1163676/1‑100 (MQ=255)
                                                                                            ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc        >  1:671131/1‑100 (MQ=255)
                                                                                             gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg       >  1:615718/1‑100 (MQ=255)
                                                                                               tCATTATAGGGAGCCAGACTTTTATTACCAATCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGat     <  1:1958564/100‑1 (MQ=255)
                                                                                               tCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGat     >  2:45733/1‑100 (MQ=255)
                                                                                               tCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGat     >  2:1842410/1‑100 (MQ=255)
                                                                                                cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGata    <  1:618755/100‑1 (MQ=255)
                                                                                                 aTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGatat   <  2:592994/100‑1 (MQ=255)
                                                                                                 aTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGatat   >  2:1085808/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATC  >  NC_000913/4470854‑4471051

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