Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,803,104 A→G pseudogene (11/426 nt) arpB → pseudogene, ankyrin repeats

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,803,1040AG97.0% 458.9 / ‑0.1 135pseudogene (11/426 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (58/73);  total (60/75)
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.98e-01

Only 100 of 137 total aligned reads displayed.

TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGC  >  NC_000913/1803007‑1803200
                                                                                                 |                                                                                                
 aaCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGa                                                                                               >  2:775402/1‑100 (MQ=255)
  aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                              >  1:2508790/1‑100 (MQ=255)
   cAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGata                                                                                             >  1:2359374/1‑100 (MQ=255)
     aTTTTTCTTTCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATc                                                                                           >  2:1034084/1‑100 (MQ=255)
     aTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATc                                                                                           <  2:830597/100‑1 (MQ=255)
       ttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCAt                                                                                         <  2:93511/100‑1 (MQ=255)
       ttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCAt                                                                                         <  1:2265295/100‑1 (MQ=255)
         ttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTa                                                                                       <  1:2462750/100‑1 (MQ=255)
         ttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTa                                                                                       <  1:2102200/100‑1 (MQ=255)
           cGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                     <  1:1921209/100‑1 (MQ=255)
           cGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCCCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                     >  2:1558807/1‑100 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                    >  1:2360160/1‑100 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                    <  2:2870034/100‑1 (MQ=255)
             tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAGGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                   >  2:1773710/1‑100 (MQ=255)
             tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                   <  2:1321297/100‑1 (MQ=255)
             tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                   <  2:820289/100‑1 (MQ=255)
             tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                   <  1:3100750/100‑1 (MQ=255)
              gCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGa                                                                                  <  1:2560250/100‑1 (MQ=255)
                acaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAAc                                                                                >  1:1572021/1‑100 (MQ=255)
                acaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAAc                                                                                <  2:2333133/100‑1 (MQ=255)
                 caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCACAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt                                                                               <  2:1775650/100‑1 (MQ=255)
                 caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt                                                                               <  1:2063350/100‑1 (MQ=255)
                 caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt                                                                               <  1:1067032/100‑1 (MQ=255)
                     tCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAtt                                                                           >  2:2280716/1‑100 (MQ=255)
                      cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATta                                                                          <  1:296624/100‑1 (MQ=255)
                      cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATta                                                                          >  2:142015/1‑100 (MQ=255)
                       gTTCTGGCGATAAGAATTAATCATCATCCTCAAAACCTACAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATtat                                                                         >  2:2011030/1‑100 (MQ=255)
                         tCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATaa                                                                       >  1:1271352/1‑100 (MQ=255)
                         gtgtgttacaatgatTAATCATCATCCTAAAACCCTCCAGGTATATATAAGGCGAACAGATTATGAGTCAAAGTGATATCATTAGCAGAACTCATTATaa                                                                       <  2:2937128/91‑1 (MQ=255)
                          cTGGCAATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  2:347989/100‑1 (MQ=255)
                           tGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  2:1890166/99‑1 (MQ=255)
                            ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                    <  1:1504092/100‑1 (MQ=255)
                             gCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  2:1890135/97‑1 (MQ=255)
                             gCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtct                                                                   >  2:2258626/1‑100 (MQ=255)
                              cGATAATAATTAATCATTATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCTATATCAGTATCAGAACTCATTATAAg                                                                      >  2:1948069/1‑96 (MQ=255)
                              cGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtct                                                                   >  2:716252/1‑99 (MQ=255)
                                ataataATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcct                                                                <  1:622784/100‑1 (MQ=255)
                                 taataaTTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcctc                                                               <  2:2343779/100‑1 (MQ=255)
                                  aataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCa                                                              <  1:116374/100‑1 (MQ=255)
                                     aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAg                                                           <  1:3089532/100‑1 (MQ=255)
                                     aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAg                                                           >  2:57573/1‑100 (MQ=255)
                                       ttaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAt                                                         <  2:956599/100‑1 (MQ=255)
                                         aatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTg                                                       <  2:1682782/100‑1 (MQ=255)
                                          atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                      <  2:657908/100‑1 (MQ=255)
                                           tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                     >  1:2777182/1‑100 (MQ=255)
                                           tcatcaccCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                     <  2:1059568/100‑1 (MQ=255)
                                            catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                    <  1:2583710/100‑1 (MQ=255)
                                             atcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcaca                                                   <  1:2261168/100‑1 (MQ=255)
                                               catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                 <  2:2626504/100‑1 (MQ=255)
                                                atcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCg                                                <  1:956800/100‑1 (MQ=255)
                                                atcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCg                                                <  2:1397395/100‑1 (MQ=255)
                                                 tcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGa                                               <  2:2037620/100‑1 (MQ=255)
                                                 tcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGa                                               <  2:2867733/100‑1 (MQ=255)
                                                 tcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGa                                               <  2:1049215/100‑1 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                              >  2:702100/1‑100 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                             >  2:2632920/1‑100 (MQ=255)
                                                     cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGTTATCATTATCAGAACTTATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                           <  2:1309247/100‑1 (MQ=255)
                                                     cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                           <  1:664798/100‑1 (MQ=255)
                                                     cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                           <  2:1943757/100‑1 (MQ=255)
                                                       aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa                                         <  1:1928319/100‑1 (MQ=255)
                                                       aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa                                         <  1:2534127/100‑1 (MQ=255)
                                                        aaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGAc                                        >  1:1430434/1‑100 (MQ=255)
                                                        aaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGAc                                        <  2:1650923/100‑1 (MQ=255)
                                                          ccTCCAGAGATATAGAAGTCGAACAGATTATTAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACAt                                      <  2:2908740/100‑1 (MQ=255)
                                                           cTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATa                                     >  2:594366/1‑100 (MQ=255)
                                                           cTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATa                                     >  1:1442165/1‑100 (MQ=255)
                                                           cTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATa                                     <  1:932525/100‑1 (MQ=255)
                                                             ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                   >  2:1344096/1‑100 (MQ=255)
                                                               aGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCaa                                 >  1:2007636/1‑100 (MQ=255)
                                                                    tatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                            >  1:1946329/1‑100 (MQ=255)
                                                                       ataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                         <  1:3242808/100‑1 (MQ=255)
                                                                       ataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                         >  2:1641232/1‑100 (MQ=255)
                                                                           ggCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCag                     >  2:3147998/1‑100 (MQ=255)
                                                                           ggCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCag                     >  2:649959/1‑100 (MQ=255)
                                                                            gCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTGCATCaga                    <  1:1711709/100‑1 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                   >  1:936006/1‑100 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                   >  1:760702/1‑100 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                   <  1:1571573/100‑1 (MQ=255)
                                                                              gAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGc                  <  2:1271338/100‑1 (MQ=255)
                                                                              gAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGc                  <  1:3230071/100‑1 (MQ=255)
                                                                              gAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGc                  <  1:2589492/100‑1 (MQ=255)
                                                                                  aGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAt              >  2:524007/1‑100 (MQ=255)
                                                                                  aGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAt              <  2:1018388/100‑1 (MQ=255)
                                                                                  aGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAt              <  1:2315313/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             >  1:1392679/1‑100 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             >  2:3247185/1‑100 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             <  2:2360133/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             >  1:2441605/1‑100 (MQ=255)
                                                                                    aacgaccgccAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTa            <  2:3046005/91‑1 (MQ=255)
                                                                                     ttATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTaa           <  2:1554783/100‑1 (MQ=255)
                                                                                       aTGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATc         <  1:3262270/100‑1 (MQ=255)
                                                                                       aTGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATc         <  1:488616/100‑1 (MQ=255)
                                                                                          aGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAt      >  2:2624915/1‑100 (MQ=255)
                                                                                           gTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt     >  2:1065442/1‑100 (MQ=255)
                                                                                           gTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt     >  1:490073/1‑100 (MQ=255)
                                                                                             cAAAGCGGTATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   >  1:1722480/1‑100 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   <  2:2219742/100‑1 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   >  2:852489/1‑100 (MQ=255)
                                                                                              aaagCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAt      >  2:2180101/1‑96 (MQ=255)
                                                                                              aaagCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgc  >  2:1732479/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGC  >  NC_000913/1803007‑1803200

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: