Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,256 (C)8→7 pseudogene (3732/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,2490C.96.6% 374.6 / 5.4 88pseudogene (3739/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (48/37);  total (49/39)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 106 total aligned reads displayed.

CACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGC  >  NC_000913/1594150‑1594344
                                                                                                   |                                                                                               
cACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTc                                                                                                 >  1:964913/1‑100 (MQ=255)
 aCCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcc                                                                                                >  2:1455273/1‑100 (MQ=255)
 aCCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcc                                                                                                <  1:2534220/100‑1 (MQ=255)
  ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                                                               >  2:556550/1‑100 (MQ=255)
  ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                                                               >  2:2211017/1‑100 (MQ=255)
  ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                                                               >  2:352864/1‑100 (MQ=255)
    tGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccc                                                                                             <  2:1389604/100‑1 (MQ=255)
    tGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccc                                                                                             >  1:2594924/1‑100 (MQ=255)
    tGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccc                                                                                             >  2:2261456/1‑100 (MQ=255)
    tGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccc                                                                                             >  1:1967218/1‑100 (MQ=255)
     gCAGTGGTGACATTTCCCTACATCACACAGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccc                                                                                             <  2:2450310/99‑1 (MQ=255)
     gCAGTGGTGACATTTCCCGACATCACACCGCTGTTATAGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccc                                                                                            >  2:1550907/1‑100 (MQ=255)
       aGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccca                                                                                          <  1:352865/100‑2 (MQ=255)
        gtggtgACGTTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccac                                                                                         <  1:2311649/100‑3 (MQ=255)
        gtggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccac                                                                                         <  1:2383408/100‑3 (MQ=255)
         tggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccacg                                                                                        <  1:154393/100‑4 (MQ=255)
         tggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCAc                                                                                        <  1:220741/100‑1 (MQ=255)
          ggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccacga                                                                                       >  1:1886703/1‑96 (MQ=255)
            tgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGata                                                                                    >  1:1273912/1‑100 (MQ=255)
             gACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGatat                                                                                   <  1:657717/100‑1 (MQ=255)
              aCATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATc                                                                                  <  2:1815543/100‑1 (MQ=255)
               cATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATcc                                                                                 >  2:217869/1‑100 (MQ=255)
                aTTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcg                                                                                >  2:2189664/1‑100 (MQ=255)
                    cccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtcgt                                                                            >  2:2383618/1‑100 (MQ=255)
                    cccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtcgt                                                                            <  1:1303520/100‑1 (MQ=255)
                     ccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTa                                                                           >  2:2367221/1‑100 (MQ=255)
                     ccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTa                                                                           >  2:1378546/1‑100 (MQ=255)
                      cGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATCAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTa                                                                           >  1:2576458/1‑99 (MQ=255)
                              cacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAc                                                                  <  2:2268865/100‑1 (MQ=255)
                              cacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAc                                                                  >  1:415730/1‑100 (MQ=255)
                                ccGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAcgc                                                                >  1:287418/1‑100 (MQ=255)
                                    tGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCcgcg                                                            >  2:1220640/1‑100 (MQ=255)
                                    tGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCcgcg                                                            >  2:2041810/1‑100 (MQ=255)
                                     gTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCgg                                                           >  2:2192996/1‑100 (MQ=255)
                                     gTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCgg                                                           <  2:2793097/100‑1 (MQ=255)
                                     gTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCgg                                                           >  1:42386/1‑100 (MQ=255)
                                       tATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                         >  2:1933877/1‑100 (MQ=255)
                                       tATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                         >  2:2394981/1‑100 (MQ=255)
                                       tATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                         >  2:2492483/1‑100 (MQ=255)
                                       naTGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                         >  1:1068449/2‑100 (MQ=255)
                                          ggATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCg                                                      >  1:2107673/1‑100 (MQ=255)
                                           gATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgt                                                     <  1:918996/100‑1 (MQ=255)
                                             tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                                   >  2:2863267/1‑100 (MQ=255)
                                               ttAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                 <  2:229305/100‑1 (MQ=255)
                                               ttAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                 >  2:1096218/1‑100 (MQ=255)
                                                  aTGTGCTGAGTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCACCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAag                                              <  2:2500391/100‑1 (MQ=255)
                                                    gtgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAagag                                            >  1:2570759/1‑100 (MQ=255)
                                                    gtgtTGACTTGCCACATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAagag                                            <  2:1967204/100‑1 (MQ=255)
                                                      gtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTc                                          >  1:2582422/1‑100 (MQ=255)
                                                       ttGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCa                                         <  2:848433/100‑1 (MQ=255)
                                                       ttGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCa                                         <  1:1487647/100‑1 (MQ=255)
                                                        tGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAg                                        >  2:121068/1‑100 (MQ=255)
                                                        tGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAg                                        <  2:2023234/100‑1 (MQ=255)
                                                        tGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAg                                        <  2:2298274/100‑1 (MQ=255)
                                                          ctctGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGa                                      >  1:2074914/4‑100 (MQ=255)
                                                            ttGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATc                                    >  1:1798963/1‑100 (MQ=255)
                                                              gCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATccc                                  <  1:731182/100‑1 (MQ=255)
                                                              gCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATccc                                  >  2:2797910/1‑100 (MQ=255)
                                                               ccATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCg                                 >  2:1567025/1‑100 (MQ=255)
                                                                cATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACCATATCCGTCGTACCCTGATACGCCGCGCCATCGCTGTTAAGAGTCAGCATccc                                  >  2:951772/1‑98 (MQ=255)
                                                                  tatTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGt                               >  1:285180/1‑99 (MQ=255)
                                                                  tatTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTAACCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTa                              >  1:625992/1‑100 (MQ=255)
                                                                   atTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAc                             <  1:2406990/100‑1 (MQ=255)
                                                                   atTAATGGCAGAGGCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAc                             >  1:1965148/1‑100 (MQ=255)
                                                                     tAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAccc                           <  2:1693394/100‑1 (MQ=255)
                                                                      aaTGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCa                          <  1:2363148/100‑1 (MQ=255)
                                                                       aTGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAg                         <  2:283035/100‑1 (MQ=255)
                                                                          gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttt                      >  1:2823362/1‑100 (MQ=255)
                                                                          gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttt                      <  2:2191840/100‑1 (MQ=255)
                                                                           cAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtttt                     <  2:1273907/100‑1 (MQ=255)
                                                                            agagTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                    >  1:2576300/1‑100 (MQ=255)
                                                                                tgagcaCCAAAAGGAATTTCCCCCCCACTATATCCGTCGTGCCCTGATACGCGGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcac                <  2:2594905/95‑1 (MQ=255)
                                                                                tCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcac                <  1:2372146/100‑1 (MQ=255)
                                                                                tCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcac                <  1:496041/100‑1 (MQ=255)
                                                                                 cGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca               <  1:686724/100‑1 (MQ=255)
                                                                                 cGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca               <  1:2526288/100‑1 (MQ=255)
                                                                                 cGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca               >  1:2896/1‑100 (MQ=255)
                                                                                 cGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca               <  1:556552/100‑1 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCCCGATATCCGTCGCACCCTCATACGCCGCGGAATCGTTGTTCCGCCCCAGGATCCCGGCACCCCCTTTTCTCACAt              >  2:461682/1‑100 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt              >  1:467133/1‑100 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt              >  2:925266/1‑100 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt              >  2:1050370/1‑100 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt              >  1:86655/1‑100 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt              >  1:2665327/1‑100 (MQ=255)
                                                                                   gAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATc             <  2:285179/100‑1 (MQ=255)
                                                                                   gAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATc             <  1:1056774/100‑1 (MQ=255)
                                                                                     aCCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCa            <  1:2709715/100‑1 (MQ=255)
                                                                                      ccGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCacc          <  1:889422/100‑1 (MQ=255)
                                                                                      ccGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCacc          <  2:2693274/100‑1 (MQ=255)
                                                                                       cGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccg         >  2:963149/1‑100 (MQ=255)
                                                                                       cGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccg         >  1:1581538/1‑100 (MQ=255)
                                                                                        gAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgc        >  2:343393/1‑100 (MQ=255)
                                                                                        gAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgc        >  2:407199/1‑100 (MQ=255)
                                                                                        gAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgc        >  1:916106/1‑100 (MQ=255)
                                                                                          aaGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgcta      <  1:1860770/100‑1 (MQ=255)
                                                                                           aGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctac     <  1:1601493/100‑1 (MQ=255)
                                                                                            gCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctacc    <  2:287417/100‑1 (MQ=255)
                                                                                            gCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctacc    >  2:1846763/1‑100 (MQ=255)
                                                                                            gCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgct       >  2:1797066/1‑97 (MQ=255)
                                                                                              aaTTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctaccgc  <  2:2823341/100‑1 (MQ=255)
                                                                                                   |                                                                                               
CACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGC  >  NC_000913/1594150‑1594344

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