Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,609,203 A→G 6.9% intergenic (‑181/+26) yneE ← / ← uxaB bestrophin family putative inner membrane protein/altronate oxidoreductase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,609,2030AG 6.9% 293.0 / 10.0 101intergenic (‑181/+26)yneE/uxaBbestrophin family putative inner membrane protein/altronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (38/56);  new base G (2/5);  total (40/61)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

Only 100 of 102 total aligned reads displayed.

AAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCTGCC  >  NC_000913/1609056‑1609352
                                                                                                                                                    |                                                                                                                                                     
aaaCGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAAc                                                                                                                                                     >  1:35354/1‑150 (MQ=255)
   cGTGATTACGATCACATTCTCAAGACGCTGCTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGTCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCgg                                                                                                                                                  >  1:198722/1‑150 (MQ=255)
   cGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCgg                                                                                                                                                  <  1:397000/150‑1 (MQ=255)
    gTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCggg                                                                                                                                                 <  2:1417503/150‑1 (MQ=255)
    gTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCggg                                                                                                                                                 >  1:398270/1‑150 (MQ=255)
    gTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCggg                                                                                                                                                 <  2:654166/150‑1 (MQ=255)
    gTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCggg                                                                                                                                                 <  1:258104/150‑1 (MQ=255)
    gTGATTAAGTTCACATTCTCTAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCggg                                                                                                                                                 <  2:238620/150‑1 (MQ=255)
       aTTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTc                                                                                                                                              <  1:402486/150‑1 (MQ=255)
        ttACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTTTCGACCGGGtttc                                                                                                                                             <  1:1597186/150‑3 (MQ=255)
         tACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTTTCGACCGGGtttct                                                                                                                                            <  1:1110971/150‑4 (MQ=255)
         tACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTTTCGACCGGGtttct                                                                                                                                            <  1:1612987/150‑4 (MQ=255)
            gATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAg                                                                                                                                         <  2:1130516/150‑1 (MQ=255)
            gATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAg                                                                                                                                         <  1:775233/150‑1 (MQ=255)
            gATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAg                                                                                                                                         >  1:706988/1‑150 (MQ=255)
            gATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAg                                                                                                                                         <  2:1501487/150‑1 (MQ=255)
             aTCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGt                                                                                                                                        >  2:1111823/1‑150 (MQ=255)
               cacaTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTaa                                                                                                                                      >  1:984154/1‑150 (MQ=255)
                    tctcAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGgtg                                                                                                                                 <  2:491822/150‑1 (MQ=255)
                     ctcAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGGTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGgtgt                                                                                                                                >  1:797912/1‑150 (MQ=255)
                         aGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGttatt                                                                                                                            <  2:484420/150‑1 (MQ=255)
                           aCGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAg                                                                                                                          <  1:1416938/150‑1 (MQ=255)
                              cTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGcac                                                                                                                       <  1:804057/150‑1 (MQ=255)
                                 cTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAAc                                                                                                                    >  1:604938/1‑150 (MQ=255)
                                 attACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATATGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAAc                                                                                                                    <  2:1378586/149‑1 (MQ=255)
                                  ttACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGTCATATAGGATGTTGTAACTAATATGGTGAAATAAGATACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACg                                                                                                                   <  1:1127144/150‑1 (MQ=255)
                                  ttACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACg                                                                                                                   <  2:173834/150‑1 (MQ=255)
                                    aCAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTTTCGACCGGGtttct                                                                                                                                            <  2:867290/123‑4 (MQ=255)
                                    aCAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTTTCGACCGGGtttct                                                                                                                                            >  1:867290/1‑120 (MQ=255)
                                     cAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGAc                                                                                                                <  2:1392371/150‑1 (MQ=255)
                                       aaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGt                                                                                                              >  1:587112/1‑150 (MQ=255)
                                        aatACCTACAAGCTTACAACATTAGCAGGCGATGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACAAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTa                                                                                                             <  1:69303/148‑1 (MQ=255)
                                           aCCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACa                                                                                                           <  1:399845/150‑1 (MQ=255)
                                             cTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTTTCGACCGGGtttctt                                                                                                                                           >  1:519808/1‑111 (MQ=255)
                                             cTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTTTCGACCGGGtttctt                                                                                                                                           <  2:519808/115‑5 (MQ=255)
                                               aCAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGCTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTc                                                                                                      >  2:23703/1‑150 (MQ=255)
                                                  aGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTcgcg                                                                                                   >  2:1461580/1‑150 (MQ=255)
                                                    cctttatcaTTAACAGGTGCCACATGCAGCGCCTGCTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCTCAACGGACGTACAGCTTCGCGCa                                                                                                 >  1:376156/8‑150 (MQ=37)
                                                    cTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCa                                                                                                 <  1:1039004/150‑1 (MQ=255)
                                                     ttACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCAt                                                                                                <  1:1111823/150‑1 (MQ=255)
                                                        tatcaTTAACAGGTGCCACATGCAGCGCCTGCTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACAGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATccc                                                                                             >  1:143191/4‑150 (MQ=38)
                                                         aacaTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCt                                                                                            <  2:5013/150‑1 (MQ=255)
                                                              tAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTc                                                                                       <  2:913034/150‑1 (MQ=255)
                                                               aGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCg                                                                                      <  2:1081092/150‑1 (MQ=255)
                                                                ccAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCga                                                                                     <  2:1274981/149‑1 (MQ=255)
                                                                  tggCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTCAAATAAGAAACCCGGTCGCATCCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCgaga                                                                                   <  2:1053645/149‑1 (MQ=255)
                                                                   ggCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGaa                                                                                  >  1:441555/1‑150 (MQ=255)
                                                                   ggCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGaa                                                                                  <  1:853720/150‑1 (MQ=255)
                                                                     cgcTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATc                                                                                <  1:1078477/150‑1 (MQ=255)
                                                                       cTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGc                                                                              >  1:619952/1‑150 (MQ=255)
                                                                          cATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATc                                                                           >  2:23093/1‑150 (MQ=255)
                                                                          cATGTGGCACCTGTTAAAGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGTTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTAAGAGAATCGCATc                                                                           >  1:276131/1‑150 (MQ=255)
                                                                           aTGTGGCACCTGTTAATGATAAATGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCCAAACCGGGTTCAGAAGTAACGGTGTAATTAGCCCAACGGACGTACCGCTTCGCGCACCCCTTTTTCGAGAATCGCATcc                                                                          <  2:82101/150‑1 (MQ=255)
                                                                            tgtgGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAg                                                                        >  1:996669/1‑150 (MQ=255)
                                                                              tgGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAgg                                                                       >  1:1158130/1‑150 (MQ=255)
                                                                               ggCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGt                                                                      >  2:154548/1‑150 (MQ=255)
                                                                                    cTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTa                                                                 >  1:1485105/1‑150 (MQ=255)
                                                                                     tGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACTGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCAATAg                                                                >  1:98318/1‑150 (MQ=255)
                                                                                     tGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAg                                                                <  1:749732/150‑1 (MQ=255)
                                                                                         aaTGTTATAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAAGCGGTTTCAGAAGTAACGGTGTTATTAGCACATCGGACGTACAGCTACGCGCATCCAGTTTTCGAGAATCGCATCCAGGTCATTAGCacc                                                            >  2:1131008/1‑148 (MQ=255)
                                                                                          aTGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAAcc                                                           >  1:1430679/1‑150 (MQ=255)
                                                                                          aTGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAAcc                                                           <  1:1284211/150‑1 (MQ=255)
                                                                                          aTGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAAcc                                                           <  2:1465839/150‑1 (MQ=255)
                                                                                            gATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCCGGTCATTAGCAACCTg                                                         >  2:1395218/1‑150 (MQ=255)
                                                                                               aaaGGCATATAGTATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTc                                                      <  1:373690/150‑1 (MQ=255)
                                                                                               aaaGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCta                                                      >  2:896044/1‑149 (MQ=255)
                                                                                               aaaGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTc                                                      <  2:198722/150‑1 (MQ=255)
                                                                                                aaGGAATATAGGATGTTGTAACTAATATGGTGATATATGAATCCCGTTCGAAACAGGGTTCAGAAGTAACGGTGTTATTAGCACTACGTACGTACAGCTTCGTGCTTCGCTTTTTCGAGATTCGCATCCAGGTCTTTAGCTACGTGCTCg                                                     <  1:1107333/150‑1 (MQ=255)
                                                                                                  ggCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgac                                                   <  1:521681/150‑1 (MQ=255)
                                                                                                      tataGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccaga                                               <  2:1232136/150‑1 (MQ=255)
                                                                                                        taGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccagacc                                             >  1:1515711/1‑150 (MQ=255)
                                                                                                          ggATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGaccagaccag                                           <  2:1636637/150‑1 (MQ=255)
                                                                                                           gATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAgg                                          <  1:1021624/150‑1 (MQ=255)
                                                                                                           gATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGGGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAgg                                          <  1:193795/150‑1 (MQ=255)
                                                                                                               ttgtAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACt                                      <  2:638047/150‑1 (MQ=255)
                                                                                                               ttgtAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACt                                      >  2:465885/1‑150 (MQ=255)
                                                                                                               ttgtAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACt                                      <  1:1718828/150‑1 (MQ=255)
                                                                                                                 gtAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTg                                    >  2:548076/1‑150 (MQ=255)
                                                                                                                  tAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGc                                   <  2:1376704/150‑1 (MQ=255)
                                                                                                                     cTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTTCTCGACCAGACCAGGCTCTTGCGTc                                >  1:377140/1‑150 (MQ=255)
                                                                                                                     cTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTc                                <  2:732990/150‑1 (MQ=255)
                                                                                                                       aaTATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAAACGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCCCGTGGGTCAg                              >  2:1129821/1‑150 (MQ=255)
                                                                                                                          atGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTc                           >  2:988696/1‑150 (MQ=255)
                                                                                                                           tGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTcc                          <  1:1396834/150‑1 (MQ=255)
                                                                                                                              tGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGc                       >  2:563627/1‑150 (MQ=255)
                                                                                                                               gAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCt                      >  1:429936/1‑150 (MQ=255)
                                                                                                                                     aGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGt                <  2:629171/150‑1 (MQ=255)
                                                                                                                                     aGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGt                <  2:398270/150‑1 (MQ=255)
                                                                                                                                     aGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGt                <  2:706988/150‑1 (MQ=255)
                                                                                                                                        aaaCCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTgg              <  2:984154/150‑1 (MQ=255)
                                                                                                                                         aaCCCGGTCGAAACCGGGTTCAGAAGTAACGGCGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGt             >  1:373418/1‑150 (MQ=255)
                                                                                                                                          aCCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTc            <  2:488454/150‑1 (MQ=255)
                                                                                                                                            ccGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCtt          >  1:621436/1‑150 (MQ=255)
                                                                                                                                              ggTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCtttt        <  1:392377/150‑1 (MQ=255)
                                                                                                                                              ggTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCtttt        <  1:472808/150‑1 (MQ=255)
                                                                                                                                                tCGAATCCGGGTTCTGAAGTAAGGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGATTCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTTGTCTTTTTc      <  1:909775/150‑1 (MQ=255)
                                                                                                                                                 cGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCACCCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCt     >  2:881985/1‑150 (MQ=255)
                                                                                                                                                  gAAACCGTGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCTg    <  2:824026/150‑1 (MQ=255)
                                                                                                                                                  gAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCTg    <  2:1419119/150‑1 (MQ=255)
                                                                                                                                                    aaCCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTATCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCTGcc  <  2:98318/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCTGCC  >  NC_000913/1609056‑1609352

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