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.5% 799.6 / ‑5.9 217pseudogene (11/426 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (98/117);  total (99/118)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 217 total aligned reads displayed.

TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTG  >  NC_000913/1803007‑1803199
                                                                                                 |                                                                                               
taACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCg                                                                                               >  2:40111/1‑100 (MQ=255)
 aaCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGa                                                                                              <  2:4042478/100‑1 (MQ=255)
  aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                             >  2:3829987/1‑100 (MQ=255)
  aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                             <  2:3337414/100‑1 (MQ=255)
   cAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGata                                                                                            >  2:1722472/1‑100 (MQ=255)
    aaTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGatat                                                                                           >  2:1641361/1‑100 (MQ=255)
    aaTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGata                                                                                            <  2:519241/99‑1 (MQ=255)
      tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCa                                                                                         <  2:1140827/100‑1 (MQ=255)
       ttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCAt                                                                                        <  1:2375592/100‑1 (MQ=255)
         ttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTa                                                                                      <  2:1767616/100‑1 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                   <  1:3449532/100‑1 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                   >  1:763428/1‑100 (MQ=255)
               cacaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGTCGAACAGATTATGAGTCAAAGCGATATCATTATCAGaa                                                                                >  2:578925/1‑100 (MQ=255)
                 caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt                                                                              >  2:3498534/1‑100 (MQ=255)
                  aTTTCGTTCTGGGGATAATAATTAATCATCATGCTCAAATCGTCCAGATATATATGATGCGAACAGATTATGAGTAAGAGCGATATAAGTATGAgta                                                                                >  2:2580646/1‑95 (MQ=255)
                  aTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTc                                                                             <  1:2954498/100‑1 (MQ=255)
                   tttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCa                                                                            <  1:2724034/100‑1 (MQ=255)
                    ttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAt                                                                           <  2:1400582/100‑1 (MQ=255)
                    ttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAt                                                                           <  2:1845387/100‑1 (MQ=255)
                      cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATta                                                                         >  1:1168686/1‑100 (MQ=255)
                         tCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATaa                                                                      <  2:3710133/100‑1 (MQ=255)
                          caggCGATAATAAATAATCATCATCCTCAAAACCTCCAGATAGATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                     >  1:2820221/3‑100 (MQ=255)
                          cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                     >  2:1409921/1‑100 (MQ=255)
                          cTGGAGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                     <  2:761484/100‑1 (MQ=255)
                           tGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                     <  2:2976347/99‑1 (MQ=255)
                            ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                   >  2:1377883/1‑100 (MQ=255)
                            ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                   <  1:208220/100‑1 (MQ=255)
                            ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                   <  1:923046/100‑1 (MQ=255)
                            ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                   <  1:2315758/100‑1 (MQ=255)
                             gCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtct                                                                  >  1:228803/1‑100 (MQ=255)
                               gATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcc                                                                >  1:1656352/1‑100 (MQ=255)
                                    taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                           >  2:2033526/1‑100 (MQ=255)
                                     aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAg                                                          <  1:3351855/100‑1 (MQ=255)
                                      attaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGa                                                         >  1:2933482/1‑100 (MQ=255)
                                          atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     >  2:1986026/1‑100 (MQ=255)
                                          atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     <  1:1307661/100‑1 (MQ=255)
                                           tcatcattCTCAAAACATTCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAAGTAATTATAAGTCTCCTCATAGATTGca                                                    >  2:308559/1‑100 (MQ=255)
                                           tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                    <  2:586401/100‑1 (MQ=255)
                                            catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                   <  1:252582/100‑1 (MQ=255)
                                            catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                   <  1:665367/100‑1 (MQ=255)
                                            catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                   <  2:3184989/100‑1 (MQ=255)
                                             atcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcaca                                                  >  1:4032847/1‑100 (MQ=255)
                                             atcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                   <  1:1236063/99‑1 (MQ=255)
                                              tcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATCTCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACAt                                                 >  2:2673931/1‑100 (MQ=255)
                                               catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                <  1:316573/100‑1 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                             >  2:3662337/1‑100 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                             <  2:2890828/100‑1 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                             >  1:983618/1‑100 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                             >  1:2130113/1‑100 (MQ=255)
                                                    tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg                                           >  1:1700292/1‑100 (MQ=255)
                                                    tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg                                           <  1:2865213/100‑1 (MQ=255)
                                                    ncAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg                                           >  1:2589092/2‑100 (MQ=255)
                                                     cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                <  2:3580291/94‑1 (MQ=255)
                                                       aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa                                        >  1:2169119/1‑100 (MQ=255)
                                                       aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa                                        <  1:3830030/100‑1 (MQ=255)
                                                       aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa                                        <  1:688776/100‑1 (MQ=255)
                                                          ccTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACAt                                     <  2:2284830/100‑1 (MQ=255)
                                                            tCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAc                                   >  2:2863062/1‑100 (MQ=255)
                                                            tCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAc                                   <  2:2097282/100‑1 (MQ=255)
                                                            tCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAc                                   >  1:3073420/1‑100 (MQ=255)
                                                             ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                  <  1:2303441/100‑1 (MQ=255)
                                                             caaGATATATATAAGGCGAACAGATGATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                  <  2:2648884/98‑1 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                               <  2:1137642/100‑1 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                               <  1:2893650/100‑1 (MQ=255)
                                                                 atatatatAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAaca                              <  2:410968/100‑1 (MQ=255)
                                                                 atatatatAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAaca                              <  1:607860/100‑1 (MQ=255)
                                                                  tatatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAacac                             >  2:1538911/1‑100 (MQ=255)
                                                                  tatatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAacac                             >  2:2514977/1‑100 (MQ=255)
                                                                   atatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACAcc                            <  2:220278/100‑1 (MQ=255)
                                                                   atatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACAcc                            <  1:388841/100‑1 (MQ=255)
                                                                       ataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                        >  1:3499525/1‑100 (MQ=255)
                                                                       ataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                        <  2:1080638/100‑1 (MQ=255)
                                                                       ataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                        >  1:2549513/1‑100 (MQ=255)
                                                                        taAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcat                       >  2:3196478/1‑100 (MQ=255)
                                                                        taAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcat                       >  1:2606019/1‑100 (MQ=255)
                                                                        taAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                        <  2:3787464/99‑1 (MQ=255)
                                                                          aGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatca                     <  2:1491528/100‑1 (MQ=255)
                                                                           ggCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCag                    >  2:2325124/1‑100 (MQ=255)
                                                                           ggCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCag                    >  2:1651885/1‑100 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                  >  2:4058013/1‑100 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                  >  1:2484157/1‑100 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                  >  1:331770/1‑100 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                  >  1:723302/1‑100 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                  >  1:1874506/1‑100 (MQ=255)
                                                                              gAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGc                 <  2:3030846/100‑1 (MQ=255)
                                                                                    aTTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTa           >  2:775050/1‑100 (MQ=255)
                                                                                    aTTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTa           <  1:4149090/100‑1 (MQ=255)
                                                                                    aTTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTa           <  1:3860117/100‑1 (MQ=255)
                                                                                    aTTATGAGTCAAAGCGATATCATTATCAGAACTCATATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTaa          >  1:3700527/1‑100 (MQ=255)
                                                                                     ttATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTaa          <  2:2750726/100‑1 (MQ=255)
                                                                                     ttATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTaa          <  1:2012273/100‑1 (MQ=255)
                                                                                     ttATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTaa          >  1:1069973/1‑100 (MQ=255)
                                                                                     ttATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTaa          >  1:1026158/1‑100 (MQ=255)
                                                                                      tATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAAt         <  1:1395204/100‑1 (MQ=255)
                                                                                       aTGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATc        <  1:1641373/100‑1 (MQ=255)
                                                                                        tGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATTAGAGCCTATTAATCa       <  1:3932637/100‑1 (MQ=255)
                                                                                          aGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAt     >  2:1197002/1‑100 (MQ=255)
                                                                                           gTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt    >  2:161087/1‑100 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg  >  2:3465725/1‑100 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg  <  1:4058056/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTG  >  NC_000913/1803007‑1803199

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: