Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,803,104 A→G 100% 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% 440.1 / ‑4.8 129pseudogene (11/426 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (62/65);  total (62/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 129 total aligned reads displayed.

TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTG  >  NC_000913/1803007‑1803199
                                                                                                 |                                                                                               
taACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCg                                                                                               <  2:1862279/100‑1 (MQ=255)
 aaCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGa                                                                                              >  1:147411/1‑100 (MQ=255)
 aaCAATGTTTCGTGCACATTTCGTTCTGGCGATAATAAGTAATCATCATCCTCAGAACCGCCAGATATATATAAGGCGAACAGATTAGGAGTCAAAGCGa                                                                                              <  2:1157380/100‑1 (MQ=255)
  aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                             >  1:3886113/1‑100 (MQ=255)
  aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                             >  1:3012594/1‑100 (MQ=255)
  aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                             >  2:3657605/1‑100 (MQ=255)
       ttttCGTGCACATTTCGTTCTGGCGATAATAATTAATGATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCAt                                                                                        >  2:2075606/1‑100 (MQ=255)
         ttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGTCGAACAGATTATGAGTCAAAGCGATATCATTa                                                                                      <  2:3169034/100‑1 (MQ=255)
         ttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTa                                                                                      <  1:1757380/100‑1 (MQ=255)
          tCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTa                                                                                      <  1:1304821/99‑1 (MQ=255)
           cGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                    <  1:1090391/100‑1 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATTCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                    <  1:1374064/99‑1 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                   <  1:2978334/100‑1 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                   >  2:1909727/1‑100 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                   <  2:21337/100‑1 (MQ=255)
             tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                  >  2:17363/1‑100 (MQ=255)
              gCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATAACATTATCAGa                                                                                 >  2:706742/1‑100 (MQ=255)
               cacaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGaa                                                                                <  1:381986/100‑1 (MQ=255)
                acaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAAc                                                                               <  2:1170112/100‑1 (MQ=255)
                 caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt                                                                              <  1:1721129/100‑1 (MQ=255)
                   tttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCa                                                                            >  2:512515/1‑100 (MQ=255)
                   tttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCa                                                                            >  2:518332/1‑100 (MQ=255)
                    ttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAt                                                                           <  2:1364392/100‑1 (MQ=255)
                    ttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAt                                                                           <  2:3167249/100‑1 (MQ=255)
                     tCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAtt                                                                          >  2:2958926/1‑100 (MQ=255)
                      cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATta                                                                         <  2:1152101/100‑1 (MQ=255)
                       gTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATtat                                                                        <  1:3535869/100‑1 (MQ=255)
                         tCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATaa                                                                      <  1:1980161/100‑1 (MQ=255)
                          cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                     >  1:1306657/1‑100 (MQ=255)
                          cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                     <  2:3703311/100‑1 (MQ=255)
                            ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                     >  2:1807591/1‑98 (MQ=255)
                                ataataATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcct                                                               <  2:3860408/100‑1 (MQ=255)
                                   ataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCAt                                                            >  2:1023212/1‑100 (MQ=255)
                                   ataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCAt                                                            <  1:3287957/100‑1 (MQ=255)
                                    taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                           >  1:3640199/1‑100 (MQ=255)
                                     aattaatCATCATCCTCAAATCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAg                                                          <  1:1716876/100‑1 (MQ=255)
                                     aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACAAATTATAAATATCCTCATa                                                           >  1:3531862/1‑99 (MQ=255)
                                      attaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGa                                                         <  1:3569130/100‑1 (MQ=255)
                                        taatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTCTCAGAACTCATTATAAGTCTCCTCATAGaat                                                       <  1:2637201/100‑3 (MQ=255)
                                        taatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAtt                                                       >  2:1980124/1‑100 (MQ=255)
                                        taatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAtt                                                       <  1:1354652/100‑1 (MQ=255)
                                        taatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAtt                                                       <  1:2119861/100‑1 (MQ=255)
                                          tctatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     <  1:366132/97‑1 (MQ=255)
                                          gtcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     <  2:505602/99‑1 (MQ=255)
                                          atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     <  1:3708756/100‑1 (MQ=255)
                                          atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     <  1:2996560/100‑1 (MQ=255)
                                          atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     <  2:54971/100‑1 (MQ=255)
                                           tcatcatcCTCATATCCTCCTGATATATATAAGGCGAACAGATTGTGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                    <  1:2477331/100‑1 (MQ=255)
                                           tcatcatcCTCAAAACCTCCAGATATATATAAGTCTAACATATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                    <  1:1133633/100‑1 (MQ=255)
                                            catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                   <  1:1529067/100‑1 (MQ=255)
                                             atcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcaga                                                  >  2:366126/1‑98 (MQ=255)
                                               catcCTCAAATCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                <  1:1454355/100‑1 (MQ=255)
                                               catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                <  2:3408527/100‑1 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                             >  2:2589699/1‑100 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                             >  1:1010882/1‑100 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTAAAAGCGATATCATTATCAGAACTCATTATAAGTATCCTCATAGATTGCACATCGAt                                             >  2:759667/1‑100 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                            <  1:3199457/100‑1 (MQ=255)
                                                    tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg                                           >  1:788008/1‑100 (MQ=255)
                                                    tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg                                           >  2:3745771/1‑100 (MQ=255)
                                                     cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                          >  1:3297413/1‑100 (MQ=255)
                                                      aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTTTAAGTCTCCTCATAGATTGCACATCGATAGCg                                         <  1:3726026/100‑1 (MQ=255)
                                                       aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGcaa                                        >  2:2377560/1‑98 (MQ=255)
                                                        aaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGAc                                       <  1:465025/100‑1 (MQ=255)
                                                          ccTCCAGATATATATAAGGCGAATAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGAc                                       <  1:530986/98‑1 (MQ=255)
                                                          ccTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACAt                                     <  2:536934/100‑1 (MQ=255)
                                                           cTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATa                                    >  2:1047476/1‑100 (MQ=255)
                                                             ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                  <  2:1663628/100‑1 (MQ=255)
                                                             ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                  <  1:1963306/100‑1 (MQ=255)
                                                              cAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCa                                 <  2:76526/100‑1 (MQ=255)
                                                              cAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCa                                 <  1:3695562/100‑1 (MQ=255)
                                                              cAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCa                                 >  1:3697864/1‑100 (MQ=255)
                                                              cAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCa                                 >  1:2516343/1‑100 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                               <  2:3739441/100‑1 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                               <  2:472197/100‑1 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                               >  1:3771266/1‑100 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                               >  1:2257266/1‑100 (MQ=255)
                                                                   atatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACAcc                            >  2:1509178/1‑100 (MQ=255)
                                                                    tatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                           >  2:2201037/1‑100 (MQ=255)
                                                                     atataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt                          >  2:642657/1‑100 (MQ=255)
                                                                     atataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt                          >  1:1990724/1‑100 (MQ=255)
                                                                      tataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt                          >  1:2766033/1‑99 (MQ=255)
                                                                        taAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcat                       >  2:3556461/1‑100 (MQ=255)
                                                                         aaGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                      >  1:2258712/1‑100 (MQ=255)
                                                                         aaGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                      <  1:479900/100‑1 (MQ=255)
                                                                            gCGAACAGATCATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga                   >  2:2758267/1‑100 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                  >  2:605356/1‑100 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                  <  1:2032596/100‑1 (MQ=255)
                                                                                 cAGATTATGAGTCAAATCTATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCt               <  1:874444/99‑1 (MQ=255)
                                                                                  aGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAt             <  2:3627990/100‑1 (MQ=255)
                                                                                  aGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAt             <  1:2589723/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt            <  2:3525573/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt            <  1:2608956/100‑1 (MQ=255)
                                                                                    ccgatctgTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTa           <  2:1091644/93‑1 (MQ=255)
                                                                                      tATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAAt         >  1:2462809/1‑100 (MQ=255)
                                                                                       aTGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGATCCTATTAATc        >  2:3418574/1‑100 (MQ=255)
                                                                                           gTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAACGAAATACCAACACCTTCATCAGAGCCTATTaagatcgg    >  1:1091655/1‑94 (MQ=255)
                                                                                           gTCAAAGCAATATCATTATCAGAACTCATTATAAGTATCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa       >  2:1046612/1‑97 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg  >  2:484593/1‑100 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg  >  2:604458/1‑100 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACAAACACCATCATAAGAGCCTATTAATCaaat    >  1:2529418/1‑96 (MQ=255)
                                                                                                 |                                                                                               
TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTG  >  NC_000913/1803007‑1803199

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