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,1040AG99.2% 425.1 / ‑4.9 118pseudogene (11/426 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (57/60);  total (57/61)
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.20e-01

Only 100 of 119 total aligned reads displayed.

AACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGC  >  NC_000913/1803008‑1803200
                                                                                                |                                                                                                
aaCAATTTTTCGTGCACATTTCGTTGTGGCGATAATAATTAATCATCATCCTCAGAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGa                                                                                               <  1:1879873/100‑1 (MQ=255)
 aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                              >  1:441129/1‑100 (MQ=255)
 aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                              >  1:334410/1‑100 (MQ=255)
 aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                              >  2:2735473/1‑100 (MQ=255)
 aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGa                                                                                               <  1:278803/99‑1 (MQ=255)
  cAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGata                                                                                             >  2:779015/1‑100 (MQ=255)
    aTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATc                                                                                           >  2:1354517/1‑100 (MQ=255)
     tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCa                                                                                          <  1:1200616/100‑1 (MQ=255)
      tgttCGTGTACATTTCGTTCTGGCGATAATAATTAATCATAATCCTCAAAGCCACCAGACAGATATACGGCGAACAGATTATGAGTCAAAGCGATATc                                                                                           <  2:1989789/96‑1 (MQ=255)
       tttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCAtt                                                                                        <  2:2714554/100‑1 (MQ=255)
       tgtCGTGCACAGTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCa                                                                                          <  1:2575840/96‑1 (MQ=255)
          cGTGCACATTTCGTTCTGGCGATAATAATTTATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                     >  2:1894372/1‑100 (MQ=255)
           gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                    >  2:3079661/1‑100 (MQ=255)
            tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                   >  2:2146914/1‑100 (MQ=255)
            tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                   <  1:75546/100‑1 (MQ=255)
             gCACATTTCGTTCTGGCTATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAg                                                                                   <  2:2975656/99‑1 (MQ=255)
             gCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGa                                                                                  <  1:1008767/100‑1 (MQ=255)
             gCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGa                                                                                  >  1:1414467/1‑100 (MQ=255)
             gCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGCCGAACAGATTATGAGTCAAAGCGATATCATTATCAGa                                                                                  <  2:1648381/100‑1 (MQ=255)
                caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt                                                                               <  1:2746690/100‑1 (MQ=255)
                     cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATta                                                                          <  1:237257/100‑1 (MQ=255)
                      gTTCTGGCGATAATAATGAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATtat                                                                         <  2:2183117/100‑1 (MQ=255)
                        tCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATaa                                                                       >  1:1508311/1‑100 (MQ=255)
                         cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  1:1327424/100‑1 (MQ=255)
                         cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  1:2918232/100‑1 (MQ=255)
                         cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  2:1771708/100‑1 (MQ=255)
                         cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  1:270724/100‑1 (MQ=255)
                          tGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  2:1332573/99‑1 (MQ=255)
                           ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                    >  2:894686/1‑100 (MQ=255)
                             cGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtctc                                                                  >  2:965101/1‑100 (MQ=255)
                             cGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtctc                                                                  >  1:1669619/1‑100 (MQ=255)
                              gATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcc                                                                 <  1:2708474/100‑1 (MQ=255)
                                taataaTTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcctc                                                               <  2:646824/100‑1 (MQ=255)
                                taataaTTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcctc                                                               <  2:573594/100‑1 (MQ=255)
                                  ataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCAt                                                             >  2:673330/1‑100 (MQ=255)
                                   taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                            >  1:2013017/1‑100 (MQ=255)
                                   taattaagCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGAGATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                            >  2:1621335/1‑100 (MQ=255)
                                   taattaacCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCCCATa                                                            <  2:226574/100‑1 (MQ=255)
                                    aattaatCATCATCCGCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTATCCTCATa                                                            >  2:2079257/1‑99 (MQ=255)
                                       taatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCGAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAtt                                                        <  1:1784998/100‑1 (MQ=255)
                                        aatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTg                                                       >  1:1125619/1‑100 (MQ=255)
                                         atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                      <  1:2051364/100‑1 (MQ=255)
                                         atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                      <  1:87586/100‑1 (MQ=255)
                                          tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                     <  1:640582/100‑1 (MQ=255)
                                           catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                     <  1:2959319/99‑1 (MQ=255)
                                             tcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACAt                                                  >  2:89645/1‑100 (MQ=255)
                                             tcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACAt                                                  <  2:1028475/100‑1 (MQ=255)
                                                tcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGa                                               >  2:2370238/1‑100 (MQ=255)
                                                 ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                              <  2:778873/100‑1 (MQ=255)
                                                  cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                             >  2:2754639/1‑100 (MQ=255)
                                                  cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                             >  1:2442264/1‑100 (MQ=255)
                                                  cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                             >  2:579465/1‑100 (MQ=255)
                                                  cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATGATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                             >  2:1974349/1‑100 (MQ=255)
                                                    naaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                           >  2:2174026/2‑100 (MQ=255)
                                                    cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                           <  2:2286214/100‑1 (MQ=255)
                                                    cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                           >  1:1384873/1‑100 (MQ=255)
                                                     aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCg                                          <  2:1669601/100‑1 (MQ=255)
                                                        aCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACa                                       <  1:818646/100‑1 (MQ=255)
                                                        aCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACa                                       <  1:323691/100‑1 (MQ=255)
                                                          cTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATa                                     >  1:446087/1‑100 (MQ=255)
                                                          cTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATa                                     >  1:2658743/1‑100 (MQ=255)
                                                            ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                   <  1:2489499/100‑1 (MQ=255)
                                                              aGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCaa                                 >  2:797977/1‑100 (MQ=255)
                                                               gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                <  1:166211/100‑1 (MQ=255)
                                                                atatatatAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAaca                               <  1:673333/100‑1 (MQ=255)
                                                                  atatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAaca                               <  2:1366461/98‑1 (MQ=255)
                                                                   tatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                            <  1:66705/100‑1 (MQ=255)
                                                                   tatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                            >  1:2587830/1‑100 (MQ=255)
                                                                    atataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt                           >  1:868393/1‑100 (MQ=255)
                                                                     tataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtc                          >  1:442316/1‑100 (MQ=255)
                                                                      ataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                         <  1:1966297/100‑1 (MQ=255)
                                                                        aaGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                       >  1:850911/1‑100 (MQ=255)
                                                                         aGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatca                      <  1:972626/100‑1 (MQ=255)
                                                                            cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                   >  1:2451664/1‑100 (MQ=255)
                                                                            cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                   >  1:1892131/1‑100 (MQ=255)
                                                                            cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                   >  2:2436530/1‑100 (MQ=255)
                                                                             gAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGc                  >  1:2942481/1‑100 (MQ=255)
                                                                               aCAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCt                <  1:3079694/100‑1 (MQ=255)
                                                                               aCAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCt                >  2:2522700/1‑100 (MQ=255)
                                                                                  gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             <  1:971951/100‑1 (MQ=255)
                                                                                  gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             <  1:972441/100‑1 (MQ=255)
                                                                                  gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             >  1:1315914/1‑100 (MQ=255)
                                                                                  gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             <  2:2202831/100‑1 (MQ=255)
                                                                                   aTTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTa            <  1:1966659/100‑1 (MQ=255)
                                                                                    ttATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTaa           <  2:441127/100‑1 (MQ=255)
                                                                                       tGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa        >  2:123441/1‑100 (MQ=255)
                                                                                        gAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCaa       >  2:959676/1‑100 (MQ=255)
                                                                                         ngTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAt      >  2:2974140/2‑100 (MQ=255)
                                                                                          gTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt     <  1:2982658/100‑1 (MQ=255)
                                                                                           tCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAttt    <  1:649878/100‑1 (MQ=255)
                                                                                            cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCACttt    >  2:243190/1‑99 (MQ=255)
                                                                                            cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   >  1:2818619/1‑100 (MQ=255)
                                                                                            cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   <  1:2735500/100‑1 (MQ=255)
                                                                                            cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   >  1:2629791/1‑100 (MQ=255)
                                                                                            cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   <  1:2605265/100‑1 (MQ=255)
                                                                                            cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   <  2:3041202/100‑1 (MQ=255)
                                                                                            cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   <  1:581648/100‑1 (MQ=255)
                                                                                            cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   <  2:1113595/100‑1 (MQ=255)
                                                                                            cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   >  2:965356/1‑100 (MQ=255)
                                                                                             aaagCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgc  >  1:1687935/1‑100 (MQ=255)
                                                                                                |                                                                                                
AACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGC  >  NC_000913/1803008‑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: ‑