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,1040AG98.2% 621.0 / ‑4.9 168pseudogene (11/426 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (70/95);  total (71/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTG  >  NC_000913/1803007‑1803199
                                                                                                 |                                                                                               
taACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCg                                                                                               <  2:2132832/100‑1 (MQ=255)
   cAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGata                                                                                            <  2:2637557/100‑1 (MQ=255)
     aTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATc                                                                                          >  1:2314846/1‑100 (MQ=255)
      tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCa                                                                                         >  1:59049/1‑100 (MQ=255)
       ttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCAt                                                                                        >  2:2012059/1‑100 (MQ=255)
           cGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                    <  1:483771/100‑1 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCGAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                   <  2:2098658/100‑1 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                   >  1:275204/1‑100 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                   <  2:732711/100‑1 (MQ=255)
             tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                  >  2:1004725/1‑100 (MQ=255)
               cacaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGaa                                                                                >  1:758802/1‑100 (MQ=255)
                acaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAAc                                                                               <  2:1662933/100‑1 (MQ=255)
                 caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt                                                                              >  1:2773049/1‑100 (MQ=255)
                  aTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTc                                                                             >  2:939735/1‑100 (MQ=255)
                    ttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAt                                                                           <  2:2900095/100‑1 (MQ=255)
                      cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATta                                                                         <  2:232427/100‑1 (MQ=255)
                       gTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATtat                                                                        >  2:2337789/1‑100 (MQ=255)
                        ttCTTGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATATGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCAtt                                                                          <  2:1444199/97‑1 (MQ=255)
                        ttCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATtata                                                                       <  1:1872902/100‑1 (MQ=255)
                            ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                   >  2:1235699/1‑100 (MQ=255)
                            ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                   <  1:1352531/100‑1 (MQ=255)
                            ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                   >  1:44323/1‑100 (MQ=255)
                              cGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGtct                                                                  <  2:2846941/99‑1 (MQ=255)
                               gATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcc                                                                <  1:797047/100‑1 (MQ=255)
                                 taataaTTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcct                                                               >  2:1559094/1‑99 (MQ=255)
                                   ataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCAt                                                            <  1:2625976/100‑1 (MQ=255)
                                   ataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCAt                                                            >  1:2397407/1‑100 (MQ=255)
                                    taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                           >  2:1704514/1‑100 (MQ=255)
                                    taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                           <  2:758799/100‑1 (MQ=255)
                                    taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                           >  1:1047544/1‑100 (MQ=255)
                                    taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                           <  1:1061694/100‑1 (MQ=255)
                                     aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                           <  1:659082/99‑1 (MQ=255)
                                       ttaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAt                                                        <  1:2118706/100‑1 (MQ=255)
                                       ttaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAt                                                        <  2:499447/100‑1 (MQ=255)
                                        taatCATCATCCTCAAAACCTCCAGAAATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAtt                                                       >  1:1815989/1‑100 (MQ=255)
                                         aatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTg                                                      >  1:2238529/1‑100 (MQ=255)
                                          atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     <  1:1757979/100‑1 (MQ=255)
                                          atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     <  2:906977/100‑1 (MQ=255)
                                          atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     >  1:403916/1‑100 (MQ=255)
                                           tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                    <  2:1082999/100‑1 (MQ=255)
                                           tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                    <  1:1000436/100‑1 (MQ=255)
                                           tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                    >  2:1335159/1‑100 (MQ=255)
                                           tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                    >  1:2058825/1‑100 (MQ=255)
                                            catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                   <  2:326958/100‑1 (MQ=255)
                                            catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                   >  1:2181132/1‑100 (MQ=255)
                                             atcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcaca                                                  >  2:2553898/1‑100 (MQ=255)
                                              tcatcCTCAAAACCTCCAGATCTATATAAGGCGAACAGATTATGAGTCAAAGCGCTATCATTATCAGAACTCATTATAACTCTCCTCATACATTGccca                                                  >  2:2706007/1‑99 (MQ=255)
                                               catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                <  1:350180/100‑1 (MQ=255)
                                               catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                <  2:1572099/100‑1 (MQ=255)
                                                atcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCg                                               >  1:2284474/1‑100 (MQ=255)
                                                 tcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGa                                              >  2:998357/1‑100 (MQ=255)
                                                 tcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGa                                              >  2:2964872/1‑100 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                            <  1:1406955/100‑1 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                            >  1:2378023/1‑100 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                            <  1:2197316/100‑1 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                            >  2:1051177/1‑100 (MQ=255)
                                                    tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg                                           >  1:1831104/1‑100 (MQ=255)
                                                      aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCg                                         <  1:1998015/100‑1 (MQ=255)
                                                        aaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGAc                                       <  1:1461066/100‑1 (MQ=255)
                                                         aCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACa                                      <  2:1009591/100‑1 (MQ=255)
                                                             ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                  <  2:2688871/100‑1 (MQ=255)
                                                              cAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCa                                 >  1:1218867/1‑100 (MQ=255)
                                                              cAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCa                                 <  2:1047539/100‑1 (MQ=255)
                                                              cAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCa                                 >  2:2063677/1‑100 (MQ=255)
                                                               aGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCaa                                >  1:2153741/1‑100 (MQ=255)
                                                               aGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCaa                                >  2:1780600/1‑100 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTTTACTCATAGATTACACATCGATAGCGACATACCa                                 >  2:2537286/1‑98 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                               <  1:2986160/100‑1 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                               <  1:227155/100‑1 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                               <  2:1973207/100‑1 (MQ=255)
                                                                 atatatatAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAaca                              <  1:1581852/100‑1 (MQ=255)
                                                                   atatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACAcc                            >  2:592279/1‑100 (MQ=255)
                                                                   atatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACAcc                            <  1:980913/100‑1 (MQ=255)
                                                                    tatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                           >  2:1362486/1‑100 (MQ=255)
                                                                     atataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt                          <  1:1993475/100‑1 (MQ=255)
                                                                      tataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtc                         >  2:489581/1‑100 (MQ=255)
                                                                        taAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCagatcggaaga                       >  1:2740854/1‑90 (MQ=255)
                                                                        taAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcat                       >  2:2117926/1‑100 (MQ=255)
                                                                        taAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                        >  1:895588/1‑99 (MQ=255)
                                                                         aaGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                      <  1:2886629/100‑1 (MQ=255)
                                                                          aGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatca                     >  2:974562/1‑100 (MQ=255)
                                                                            gCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga                   <  2:2397392/100‑1 (MQ=255)
                                                                            gCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga                   <  2:22655/100‑1 (MQ=255)
                                                                               aaCAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGcc                <  1:527816/100‑1 (MQ=255)
                                                                                 cAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTa              >  1:2349315/1‑100 (MQ=255)
                                                                                  aGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAt             <  1:1381805/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAGAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt            >  1:1961195/1‑100 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt            <  2:2378009/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt            <  1:1559103/100‑1 (MQ=255)
                                                                                    aTTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTa           <  2:1772573/100‑1 (MQ=255)
                                                                                      tATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAAt         >  2:2223742/1‑100 (MQ=255)
                                                                                       aTGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATc        <  1:2509123/100‑1 (MQ=255)
                                                                                       aTGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATc        <  2:2482329/100‑1 (MQ=255)
                                                                                        tGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa       >  1:1312140/1‑100 (MQ=255)
                                                                                         gAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCaa      <  1:1157576/100‑1 (MQ=255)
                                                                                          aGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAt     >  1:130201/1‑100 (MQ=255)
                                                                                           gTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt    >  1:1910358/1‑100 (MQ=255)
                                                                                           gTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt    <  1:2964888/100‑1 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg  <  1:1144704/100‑1 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg  >  2:803104/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTG  >  NC_000913/1803007‑1803199

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