Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,594,256 (C)8→7 100% 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.91.4% 360.3 / 22.5 93pseudogene (3739/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base . (41/44);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

CACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACC  >  NC_000913/1594150‑1594342
                                                                                                   |                                                                                             
cACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTc                                                                                               <  2:2108641/100‑1 (MQ=255)
 aCCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcc                                                                                              <  2:1129790/100‑1 (MQ=255)
 aCCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATAGCAAAGTCGGAACCGAAAGCAATTTcc                                                                                              >  1:2000945/1‑100 (MQ=255)
  ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                                                             >  1:565425/1‑100 (MQ=255)
    tGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccc                                                                                           >  1:2008049/1‑100 (MQ=255)
    tGCAGTGGTGACATTTCCCGACATCACACCGCTGATATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccc                                                                                           <  1:2061824/100‑1 (MQ=255)
      cAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccc                                                                                         >  2:2353259/1‑100 (MQ=255)
      cAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccc                                                                                         <  2:1630869/100‑1 (MQ=255)
        gtggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccac                                                                                       >  1:1189825/1‑98 (MQ=255)
         tggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccacg                                                                                      <  1:148398/100‑4 (MQ=255)
          ggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATCTTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccacga                                                                                     >  1:1066665/1‑96 (MQ=255)
          ggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccacga                                                                                     >  1:1514848/1‑96 (MQ=255)
              aCATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATc                                                                                >  1:1677336/1‑100 (MQ=255)
               cATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATcc                                                                               >  2:1932081/1‑100 (MQ=255)
               cATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATcc                                                                               >  1:2142875/1‑100 (MQ=255)
                aTTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATcc                                                                               >  2:2375243/1‑100 (MQ=255)
                aTTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcg                                                                              <  1:1884977/100‑1 (MQ=255)
                 tttCCCGACATCACACCGCTGTTATGGATATTACTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcg                                                                              >  1:2364162/1‑99 (MQ=255)
                 tttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgt                                                                             <  2:1685532/100‑1 (MQ=255)
                 tttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgt                                                                             <  1:1806143/100‑1 (MQ=255)
                 tttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgt                                                                             <  2:1410758/100‑1 (MQ=255)
                  ttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtc                                                                            >  1:887809/1‑100 (MQ=255)
                    cccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtcgt                                                                          >  1:1811157/1‑100 (MQ=255)
                      cgtaaTCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTGCCCCCCCACGATATCCGTCGTAc                                                                        <  1:562566/96‑1 (MQ=255)
                      cGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTAc                                                                        <  2:2142889/100‑1 (MQ=255)
                          aTCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTg                                                                    <  2:1982594/100‑1 (MQ=255)
                          aTCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTg                                                                    >  1:1675133/1‑100 (MQ=255)
                            cacaCCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGCACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGAt                                                                  >  2:2439794/1‑100 (MQ=255)
                            cacaCCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGAt                                                                  <  1:748146/100‑1 (MQ=255)
                            cacaCCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGAt                                                                  >  2:396471/1‑100 (MQ=255)
                            cacaCCGCTGTTATGCATATTAATGTGTTGACTCGCCATATTAATGGCAGAGTCTGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGa                                                                   <  1:445285/99‑1 (MQ=255)
                             acacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATa                                                                 <  1:1132144/100‑1 (MQ=255)
                                ccGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAcgc                                                              >  2:616880/1‑100 (MQ=255)
                                    ttgtATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAcgccg                                                            <  1:276288/95‑1 (MQ=255)
                                     gTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCgg                                                         <  1:2439777/100‑1 (MQ=255)
                                     gTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCgg                                                         <  2:1061763/100‑1 (MQ=255)
                                     gTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCgg                                                         <  2:1316460/100‑1 (MQ=255)
                                     gTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCgg                                                         <  2:1932251/100‑1 (MQ=255)
                                     gTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCgg                                                         >  2:775974/1‑100 (MQ=255)
                                       tATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACTCCGCGGaa                                                       >  1:831181/1‑100 (MQ=255)
                                       tATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                       >  1:1824382/1‑100 (MQ=255)
                                       tATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGa                                                        >  2:2170933/1‑99 (MQ=255)
                                       tATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAACTTCCCCCCCACGATATCCGTCGTACCCTGATAcgccgc                                                           >  1:1782732/1‑96 (MQ=255)
                                         tGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATc                                                     >  1:130184/1‑100 (MQ=255)
                                         tGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACCCCGCGGCATc                                                     >  1:1877911/1‑100 (MQ=255)
                                          ggATATTCATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCg                                                    >  1:325182/1‑100 (MQ=255)
                                          ggATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCg                                                    <  1:2311608/100‑1 (MQ=255)
                                          ggATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCg                                                    <  2:1427440/100‑1 (MQ=255)
                                           gATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCg                                                    <  2:203548/100‑1 (MQ=255)
                                             tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAACGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgt                                                   <  1:1419297/98‑1 (MQ=255)
                                             tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                                 <  2:458863/100‑1 (MQ=255)
                                             tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                                 >  2:1812860/1‑100 (MQ=255)
                                             tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                                 <  2:153250/100‑1 (MQ=255)
                                             tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                                 <  1:2496463/100‑1 (MQ=255)
                                               ttAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                                 >  2:1531328/1‑98 (MQ=255)
                                                tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                              <  2:1811167/100‑1 (MQ=255)
                                                tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                              <  2:2198432/100‑1 (MQ=255)
                                                 aaTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                              <  2:1038866/100‑1 (MQ=255)
                                                  aTGTGTTGACTTGCCATATTAATGGCAGATTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAag                                            <  1:1896813/100‑1 (MQ=255)
                                                  aTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAag                                            <  2:130185/100‑1 (MQ=255)
                                                  aTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAag                                            <  1:2527115/100‑1 (MQ=255)
                                                   tgtgTTGACTTGCCATATGCATGGCAGAGTGGGGTTTGAACGCAATTCTCCTCCCCCGATATACGACGTACCCTGATACGCCCCCGAACGGTTCTTAaga                                           >  1:1700122/1‑100 (MQ=255)
                                                    gtgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAagag                                          <  2:1273258/100‑1 (MQ=255)
                                                      gtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTc                                        >  2:2209553/1‑100 (MQ=255)
                                                      gtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTc                                        >  2:693699/1‑100 (MQ=255)
                                                        tGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAg                                      <  2:1003494/100‑1 (MQ=255)
                                                         gACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAgg                                     <  1:498753/100‑1 (MQ=255)
                                                              gCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATccc                                >  2:368012/1‑100 (MQ=255)
                                                              gCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATccc                                >  2:2073551/1‑100 (MQ=255)
                                                              gCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTACCCCCCACGATATCAGTCGTACCCTGATACGCCGCGGACTCGTTGTTAAGAGTCAGGATccc                                >  1:758027/1‑100 (MQ=255)
                                                                cATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCgg                              <  2:236752/100‑1 (MQ=255)
                                                                cATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCgg                              <  1:1485393/100‑1 (MQ=255)
                                                                  tatTAGTGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCTTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTa                            <  1:1983946/100‑1 (MQ=255)
                                                                  tatTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTAAGGATCCCgg                              >  1:633680/1‑98 (MQ=255)
                                                                   atTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAc                           <  2:10133/100‑1 (MQ=255)
                                                                         ggCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtt                     >  2:1991690/1‑100 (MQ=255)
                                                                         ggCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGt                      >  2:1269606/1‑99 (MQ=255)
                                                                          gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttt                    <  2:1630596/100‑1 (MQ=255)
                                                                           cAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtttt                   >  1:955978/1‑100 (MQ=255)
                                                                            cgagTCAGTACCGAAAGCAATTTCCCCCACACGATATCCGTCGTACCCTGAGCCGCCGCGGAATCGTTGTTAAGAGTCAGGAGCCCGGTACCCAGttt                    <  1:2325342/97‑1 (MQ=255)
                                                                             gagTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttttn                  <  1:2375228/100‑2 (MQ=255)
                                                                              agTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTc                <  2:606180/100‑1 (MQ=255)
                                                                               gTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTca               <  1:880641/100‑1 (MQ=255)
                                                                               gTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTca               >  1:688787/1‑100 (MQ=255)
                                                                               gTCCGACCCGAAAGCAATTTCCCCCCCCCGATATCCGTCCTACCCTGATACGCCGCAGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGt                 <  1:2103074/98‑1 (MQ=255)
                                                                                tCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcac              <  2:698314/100‑1 (MQ=255)
                                                                                 cGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca             >  2:1494567/1‑100 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt            >  1:196990/1‑100 (MQ=255)
                                                                                   gAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAACAGTCAGGATCCCGGTACCCAGTTTTGTCACATc           >  1:294403/1‑100 (MQ=255)
                                                                                    aaCCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCa          <  1:970801/100‑1 (MQ=255)
                                                                                     aCCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCac         >  1:606920/1‑100 (MQ=255)
                                                                                     aCCGAAAGCAATTGCCCCCCCACGATATCCGTCGTACCGTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATTCCGGTACCCAGTTTTGTCACATCac         >  1:1595086/1‑100 (MQ=255)
                                                                                        gAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgc      <  2:327064/100‑1 (MQ=255)
                                                                                        gAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgc      >  2:1029408/1‑100 (MQ=255)
                                                                                        gAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccg       <  1:1168146/99‑1 (MQ=255)
                                                                                          aaGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgcta    >  2:1240603/1‑100 (MQ=255)
                                                                                           aGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctac   <  1:2209538/100‑1 (MQ=255)
                                                                                            gCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctacc  <  2:1211964/100‑1 (MQ=255)
                                                                                            gCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctacc  <  1:616877/100‑1 (MQ=255)
                                                                                            gCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctacc  >  2:1008652/1‑100 (MQ=255)
                                                                                                   |                                                                                             
CACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACC  >  NC_000913/1594150‑1594342

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: