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,1040AG97.6% 460.7 / ‑2.4 136pseudogene (11/426 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base G (56/76);  total (56/80)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGC  >  NC_000913/1803007‑1803200
                                                                                                 |                                                                                                
taACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCg                                                                                                >  2:1387241/1‑100 (MQ=255)
taACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCg                                                                                                <  1:1448697/100‑1 (MQ=255)
  aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                              >  2:1779150/1‑100 (MQ=255)
  aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                              >  2:3772796/1‑100 (MQ=255)
   cAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGata                                                                                             >  2:3432347/1‑100 (MQ=255)
    naTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGAGAACAGACTATAAGACAAAGCGat                                                                                              >  1:3565160/2‑98 (MQ=255)
    aaTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGatat                                                                                            >  2:1301116/1‑100 (MQ=255)
    aaTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGatat                                                                                            >  1:607280/1‑100 (MQ=255)
     aTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAAGCATCATCCTCAAAACCTCCAAATATATATAAGGCGAACAGATTATGAGTCAAAGCg                                                                                                <  1:3234116/95‑1 (MQ=255)
       ttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCa                                                                                          <  1:279936/99‑1 (MQ=255)
         tgcGTTCACATTTCGTTCTGCGGATAATAATTAATCATCATCCTCAAAACCTCCAGACATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTa                                                                                       <  1:3112085/98‑1 (MQ=255)
           cGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                     >  1:2846761/1‑100 (MQ=255)
           cGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                     >  2:2652363/1‑100 (MQ=255)
             tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                   <  1:3777217/100‑1 (MQ=255)
             tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                   <  2:1347089/100‑1 (MQ=255)
             tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAg                                                                                   <  1:1961563/100‑1 (MQ=255)
              gCACATTTCGTTCTGTCGATAATAATTAACAATCATACTCAACAAATACAGATATATATAAGGCGCAAAAATTATAAGTCAAAAAGATATCATTATCAg                                                                                   <  1:2211794/99‑1 (MQ=255)
                acaTTTCGTTCTGGCGATAATAAATAATCAGCATCCTCAACACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAAc                                                                                <  1:3045309/100‑1 (MQ=255)
                 caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt                                                                               <  2:1709290/100‑1 (MQ=255)
                   tttCGTTCTGGAGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCa                                                                             <  1:1887454/100‑1 (MQ=255)
                    ttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAt                                                                            <  2:2263079/100‑1 (MQ=255)
                    ttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAt                                                                            <  1:904657/100‑1 (MQ=255)
                     tCGTTCTGGCGATAATAATTAATCAACATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAtt                                                                           <  1:2852521/100‑1 (MQ=255)
                       gTTCTTGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATtat                                                                         <  1:405956/100‑1 (MQ=255)
                         tCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATaa                                                                       >  1:163794/1‑100 (MQ=255)
                          cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  2:2371767/100‑1 (MQ=255)
                          cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  2:929637/100‑1 (MQ=255)
                          cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                      <  2:90934/100‑1 (MQ=255)
                            ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                    <  2:471278/100‑1 (MQ=255)
                              cGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtctc                                                                  >  1:3153612/1‑100 (MQ=255)
                               gATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcc                                                                 <  2:1402504/100‑1 (MQ=255)
                                 taataaTTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcctc                                                               <  1:1444007/100‑1 (MQ=255)
                                 taataaTTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcctc                                                               >  2:3301703/1‑100 (MQ=255)
                                 taataaTTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTctcctc                                                               <  1:402783/100‑1 (MQ=255)
                                    taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                            >  1:2375056/1‑100 (MQ=255)
                                    taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                            <  2:1581434/100‑1 (MQ=255)
                                    taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                            <  2:71153/100‑1 (MQ=255)
                                     aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAg                                                           >  1:1990799/1‑100 (MQ=255)
                                      attaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGa                                                          <  1:2940526/100‑1 (MQ=255)
                                         aatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTg                                                       <  1:1379311/100‑1 (MQ=255)
                                         aatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTg                                                       <  1:3497463/100‑1 (MQ=255)
                                          atcatcatcCTCAACACCTCCAGCTATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTg                                                       <  2:706727/99‑1 (MQ=255)
                                          atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGGGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                      <  2:1830656/100‑1 (MQ=255)
                                            catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                    <  1:2775545/100‑1 (MQ=255)
                                             atcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcaca                                                   <  2:1003972/100‑1 (MQ=255)
                                             atcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcaca                                                   <  2:3380809/100‑1 (MQ=255)
                                                 tcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAATCACAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATATATTGCACATc                                                 >  1:1489218/1‑98 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                              >  2:2826024/1‑100 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCAAACAGATTATGAGTCAAAGCGATATTATCATCAGAACTCATTATAATTCTCTTCATAGATTGCACATCGAt                                              >  1:485554/1‑100 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                             >  1:2076816/1‑100 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTAAGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTACTCATAGATTGCACATc                                                 >  1:1326198/1‑96 (MQ=255)
                                                    tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg                                            <  1:2707593/100‑1 (MQ=255)
                                                     cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                           <  2:2541513/100‑1 (MQ=255)
                                                     cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                           <  1:837010/100‑1 (MQ=255)
                                                     cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                           <  2:2612705/100‑1 (MQ=255)
                                                      aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                           <  1:934879/99‑1 (MQ=255)
                                                      aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCg                                          >  2:2915254/1‑100 (MQ=255)
                                                       aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa                                         <  2:3128818/100‑1 (MQ=255)
                                                         aCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACa                                       <  2:2736662/100‑1 (MQ=255)
                                                          ccTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACAt                                      >  1:1519393/1‑100 (MQ=255)
                                                          ccTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACAt                                      >  2:2652203/1‑100 (MQ=255)
                                                          ccTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATAACATTATCAGAACTCATTATAAGTCTACTCATAGATTGAACATCGATATCGACAt                                      >  1:2571609/1‑100 (MQ=255)
                                                           cTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATa                                     >  2:3576634/1‑100 (MQ=255)
                                                           cTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTACACATCGATAGCGACATa                                     <  1:2719580/100‑1 (MQ=255)
                                                            tCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAc                                    >  2:1221957/1‑100 (MQ=255)
                                                            tCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAc                                    >  2:1221936/1‑100 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                >  1:1444082/1‑100 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                <  2:2267369/100‑1 (MQ=255)
                                                                    tatataAGGCGAACAGATTATGGGACAAAGCGATATCATTATCAGAACTCATTATAAGTTGCCTCATAGATTGCACATCGATATCGAGATACaaa                                 >  1:1001558/1‑95 (MQ=255)
                                                                    tatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                            <  1:3383570/100‑1 (MQ=255)
                                                                     atataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt                           <  1:1640194/100‑1 (MQ=255)
                                                                      tataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtc                          >  2:3006206/1‑100 (MQ=255)
                                                                       ataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                         <  1:107910/100‑1 (MQ=255)
                                                                        taAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGTACATCGATAGCGACATACCAACACCTtcat                        >  2:12793/1‑100 (MQ=255)
                                                                         aaGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                       >  2:3127765/1‑100 (MQ=255)
                                                                         aaGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                       <  1:1212710/100‑1 (MQ=255)
                                                                         aaGGCGAACAGATTATGAATCAAAGCTATATCATTATCAGAACTCATTACAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                         <  1:2094188/98‑1 (MQ=255)
                                                                          aGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatca                      >  1:938029/1‑100 (MQ=255)
                                                                           ggCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCGTCATAGATTGCACATCGATAGCGACATGTCAACACCTtca                         >  1:3135695/1‑96 (MQ=255)
                                                                             cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                   <  2:3668424/100‑1 (MQ=255)
                                                                              gAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGc                  >  1:3381478/1‑100 (MQ=255)
                                                                              gAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGc                  <  1:649460/100‑1 (MQ=255)
                                                                               aaCAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGcn                 <  1:2352867/100‑2 (MQ=255)
                                                                               aaCAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGcc                 <  1:1779143/100‑1 (MQ=255)
                                                                               aaCAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGcc                 <  1:2701694/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             >  2:1697316/1‑100 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             <  1:2907868/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt             <  2:552903/100‑1 (MQ=255)
                                                                                     ttATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTaa           >  2:2699179/1‑100 (MQ=255)
                                                                                        tGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa        <  2:3332960/100‑1 (MQ=255)
                                                                                         gAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCaa       <  2:3207774/100‑1 (MQ=255)
                                                                                          aGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAt      >  1:789711/1‑100 (MQ=255)
                                                                                          aGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAt      >  2:1428676/1‑100 (MQ=255)
                                                                                           gTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt     <  2:456555/100‑1 (MQ=255)
                                                                                           gTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAGTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATaaa       >  1:315744/1‑98 (MQ=255)
                                                                                            tCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAttt    <  1:2979580/100‑1 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAttt    >  1:1356813/1‑99 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   <  2:2045309/100‑1 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg   >  2:840638/1‑100 (MQ=255)
                                                                                              aaagCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgc  >  2:2972181/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGC  >  NC_000913/1803007‑1803200

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: