Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,420,057 C→G 6.3% pseudogene (1016/1025 nt) yhdW → putative ABC transporter periplasmic binding protein YhdW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,420,0570CG 6.3% 201.6 / 2.4 61pseudogene (1016/1025 nt)yhdWputative ABC transporter periplasmic binding protein YhdW
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base G (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TAATAAATGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTCGCGCCTGGTTATTCC  >  NC_000913/3419913‑3420206
                                                                                                                                                |                                                                                                                                                     
taataaATGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACGGGCGGTGCg                                                                                                                                                  >  1:643742/1‑150 (MQ=255)
taataaATGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCg                                                                                                                                                  >  1:243734/1‑150 (MQ=255)
taataaATGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCg                                                                                                                                                  <  1:208477/150‑1 (MQ=255)
  ttaaTTGGGGCTATAAGATGATTTAAGAGGTGGGTAAGTAGTCGGAATTTTTTGAGCGTAAGGTTGGTTCAGAAAGGCCGCTGAAAATTAAACGTGGGCAAAATAATCTGTGGAATAACGGCGGTATTCAGTACGCACCGCCGGTGCGtt                                                                                                                                                <  2:665382/149‑1 (MQ=255)
  ataaATGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGTAATAACGTCGGTATTCAGTACGCACCGCCCGTGCGtt                                                                                                                                                <  2:1360166/150‑1 (MQ=255)
   taaATGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACGGGCGGTGCGtac                                                                                                                                               <  2:643742/150‑3 (MQ=255)
   taaATGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTa                                                                                                                                               <  2:116139/150‑1 (MQ=255)
   taaATGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTa                                                                                                                                               <  1:896266/150‑1 (MQ=255)
   taaATGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTa                                                                                                                                               <  1:1265156/150‑1 (MQ=255)
     aaTGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAg                                                                                                                                             <  1:737331/150‑1 (MQ=255)
         ggCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCgg                                                                                                                                         >  2:947660/1‑150 (MQ=255)
         ggCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCgg                                                                                                                                         >  1:1251368/1‑150 (MQ=255)
           ccTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGgtg                                                                                                                                       <  2:1413456/150‑1 (MQ=255)
             tataACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGtgtg                                                                                                                                     >  1:897122/1‑150 (MQ=255)
              ataACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGa                                                                                                                                    >  2:874151/1‑150 (MQ=255)
               taACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAg                                                                                                                                   <  2:790335/150‑1 (MQ=255)
                 aCATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGt                                                                                                                                 >  1:1014882/1‑150 (MQ=255)
                 aCATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGt                                                                                                                                 >  1:236092/1‑150 (MQ=255)
                  catcatCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTa                                                                                                                                >  1:553381/1‑150 (MQ=255)
                  catcatCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACCCACCGCCCGTGCGTTAAGTCGGTGTGAGGTa                                                                                                                                >  2:1109564/1‑150 (MQ=255)
                     catcaAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACg                                                                                                                             >  1:616678/1‑150 (MQ=255)
                                tGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTAc                                                                                                                  >  2:1337309/1‑150 (MQ=255)
                                 gggTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTAcg                                                                                                                 >  2:1011937/1‑150 (MQ=255)
                                  ggTAACTTCTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGtgtg                                                                                                                                     <  2:1065261/129‑1 (MQ=255)
                                  ggTAACTTCTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGtgtg                                                                                                                                     >  1:1065261/1‑129 (MQ=255)
                                   gTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTAcggc                                                                                                               <  1:1109564/150‑1 (MQ=255)
                                        tactCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGgcgc                                                                                                          <  2:553381/150‑1 (MQ=255)
                                           tCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGctct                                                                                                       <  1:947660/150‑1 (MQ=255)
                                            cGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGctctc                                                                                                      >  2:654100/1‑150 (MQ=255)
                                               aaaTTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTc                                                                                                   >  1:1351660/1‑150 (MQ=255)
                                                  tttttGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGc                                                                                                >  2:40606/1‑150 (MQ=255)
                                                   ttttGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCg                                                                                               >  1:803107/1‑150 (MQ=255)
                                                    tttGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGt                                                                                              >  1:1251298/1‑150 (MQ=255)
                                                       gAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTAt                                                                                           <  1:1356961/150‑1 (MQ=255)
                                                           gTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCtt                                                                                       <  2:1014882/150‑1 (MQ=255)
                                                                    gTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGtt                                                                              >  1:1258172/1‑150 (MQ=255)
                                                                      tCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTc                                                                            <  2:616678/150‑1 (MQ=255)
                                                                       cAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCt                                                                           >  2:778939/1‑150 (MQ=255)
                                                                          aaaGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATg                                                                        <  1:188437/150‑1 (MQ=255)
                                                                            aGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGtc                                                                      <  2:66387/150‑1 (MQ=255)
                                                                                   tGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATcgc                                                               >  1:511825/1‑150 (MQ=255)
                                                                                   tGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATcgc                                                               >  1:417836/1‑150 (MQ=255)
                                                                                           aaaCGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAAc                                                       >  2:468252/1‑150 (MQ=255)
                                                                                            aaCGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAAcc                                                      >  2:1163648/1‑150 (MQ=255)
                                                                                                     aaaaTAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGc                                             <  2:897122/150‑1 (MQ=255)
                                                                                                     aaaaTAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGc                                             <  2:511825/150‑1 (MQ=255)
                                                                                                     aaaaTAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGc                                             <  1:778939/150‑1 (MQ=255)
                                                                                                             ctctGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCActctc                                     >  1:1257149/1‑150 (MQ=255)
                                                                                                                 ggAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCtttt                                 <  1:1011937/150‑1 (MQ=255)
                                                                                                                         ggcggTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGGTTCTATGGCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAAccc                         <  1:916115/150‑1 (MQ=255)
                                                                                                                           cggTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTa                       <  1:654100/150‑1 (MQ=255)
                                                                                                                            ggTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTAc                      <  1:1037159/150‑1 (MQ=255)
                                                                                                                             gTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACg                     >  1:1269757/1‑150 (MQ=255)
                                                                                                                                ttCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGtt                  >  1:1161660/1‑150 (MQ=255)
                                                                                                                                    taaCGCACGGGCGGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCACCCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTcgcg              >  1:925209/3‑150 (MQ=255)
                                                                                                                                     tACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTCgcgc             >  1:886844/1‑150 (MQ=255)
                                                                                                                                       cGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTCGCGCCt           <  2:1136040/150‑1 (MQ=255)
                                                                                                                                       cGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTCGCGCCt           >  1:976033/1‑150 (MQ=255)
                                                                                                                                       cGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTCGCGCCt           <  1:40606/150‑1 (MQ=255)
                                                                                                                                         cACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTCGCGCCTgg         >  1:391136/1‑150 (MQ=255)
                                                                                                                                          aCCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTCGCGCCTGGt        >  1:1139061/1‑150 (MQ=255)
                                                                                                                                                cGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTCGCGCCTGGTTATTcc  <  2:167226/150‑1 (MQ=255)
                                                                                                                                                cGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTCGCGCCTGGTTATTcc  <  2:1251368/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TAATAAATGGGCCTATAACATCATCAAACAGGTGGGTAACTACTCGGAAATTTTTGAGCGTAACGTAGGTTCAGAAAGCCCGCTGAAAATTAAACGTGGGCAAAATAATCTCTGGAATAACGGCGGTATTCAGTACGCACCGCCCGTGCGTTAAGTCGGTGTGAGGTAACGAGCGTCGCTACGGCGGCGCTCTCTTCAGCGTTATGCTTTCCGAGGTTTCTATGTCTCATCGCCGCTCAACCGTTAAAGGCTCACTCTCTTTTGCCAACCCTACGGTTCGCGCCTGGTTATTCC  >  NC_000913/3419913‑3420206

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