Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 556,641 G→T 6.6% intergenic (‑88/+20) ybcI ← / ← ybcJ DUF457 family inner membrane protein/ribosome‑associated protein; putative RNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913556,6410GT 6.6% 329.0 / 9.7 116intergenic (‑88/+20)ybcI/ybcJDUF457 family inner membrane protein/ribosome‑associated protein; putative RNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (61/47);  new base T (3/5);  total (64/52)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 120 total aligned reads displayed.

TTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTGC  >  NC_000913/556492‑556790
                                                                                                                                                     |                                                                                                                                                        
ttCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTg                                                                                                                                                          >  1:621069/1‑150 (MQ=255)
ttCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTACGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTg                                                                                                                                                          <  1:1196908/150‑1 (MQ=255)
 tCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGa                                                                                                                                                         >  2:1637238/1‑150 (MQ=255)
 tCGAACCCAACACTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTAGa                                                                                                                                                         >  1:1945096/1‑150 (MQ=255)
   gAACCCAACCCTTAACCAATGCATTGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGTTGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATa                                                                                                                                                       <  1:425420/150‑1 (MQ=255)
   gAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTGGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATa                                                                                                                                                       >  1:763292/1‑150 (MQ=255)
   gAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGGCGCTGCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATa                                                                                                                                                       <  1:1836477/150‑1 (MQ=255)
    aacccaacccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATAc                                                                                                                                                      >  1:1597194/1‑150 (MQ=255)
    aacccaacccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGAGGAACGGGGCGTAGAGGGCAAAAGTCTGATAAGAGAACCGGCGTGTTGATAc                                                                                                                                                      >  2:1962453/1‑150 (MQ=255)
           cccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCgg                                                                                                                                               <  1:1990312/150‑1 (MQ=255)
            ccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCggt                                                                                                                                              <  1:573190/150‑2 (MQ=255)
            ccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCggt                                                                                                                                              <  1:964500/150‑2 (MQ=255)
             cTAAACCAATGCCTAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGa                                                                                                                                             <  2:1772859/150‑1 (MQ=255)
             cTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCCATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGa                                                                                                                                             >  1:922979/1‑150 (MQ=255)
              tAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCggttc                                                                                                                                            <  2:1905588/150‑4 (MQ=255)
              aaaaCCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCggttc                                                                                                                                            <  2:511909/149‑4 (MQ=255)
                aaCCACTGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCAGCCTGTTGATACAGGCCGGGAAA‑‑‑g                                                                                                                                       <  2:1687471/150‑1 (MQ=255)
                aaCCAATGCATAGTGTAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCTTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑g                                                                                                                                       <  1:391323/150‑1 (MQ=255)
                      tGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCa                                                                                                                                 >  1:254363/1‑150 (MQ=255)
                       gCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAg                                                                                                                                >  2:300248/1‑150 (MQ=255)
                         aTAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGTGAAA‑‑‑GGGATCAGGc                                                                                                                              >  1:1455838/1‑150 (MQ=255)
                         aTAGGGGAACGGCGGCGTGCGTAATAACGCTTGTCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGGCTGAAAAGAGAACCGGCCTGATGATACAGGCCGGGAAA‑‑‑GGGATCAGGc                                                                                                                              <  1:62575/150‑1 (MQ=255)
                              ggAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGcaacaa                                                                                                                         >  1:1297784/1‑150 (MQ=255)
                                      ggcgggcgTAAAAACGGTTGGCATGATGGTCGCTTCGGCAAAAAGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCGGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACg                                                                                                                 <  1:1400792/150‑1 (MQ=255)
                                      ggcgTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACg                                                                                                                 >  2:1199611/1‑150 (MQ=255)
                                      ggcgTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACg                                                                                                                 >  1:255684/1‑150 (MQ=255)
                                       gcgtgcgtAATAACGGTTGGCATGATGGGCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGc                                                                                                                <  1:435750/150‑1 (MQ=255)
                                        cgtgcgtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCt                                                                                                               >  2:165915/1‑150 (MQ=255)
                                           gcgtAATAACGGGGGGCATGATGGTCGCTGCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCtgtg                                                                                                            <  2:1536276/150‑1 (MQ=255)
                                              taataaCGGTTGGCATGAAGGTCGCTTCGGCAAAATGTCGAGCCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGACTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGAcc                                                                                                         <  1:1349509/150‑1 (MQ=255)
                                                 taaCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGGTGATAAAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACCCTGTGACCTGc                                                                                                      >  2:1841559/1‑150 (MQ=255)
                                                        tGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTc                                                                                               <  1:555624/150‑1 (MQ=255)
                                                          gCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGTCGTAGAGTGCAAAAGTCTGAAACGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCAc                                                                                             <  2:1039431/150‑1 (MQ=255)
                                                           cATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACt                                                                                            >  1:1108906/1‑150 (MQ=255)
                                                             tgatgGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGCCCTGCAAAACTCACTGt                                                                                          >  1:1124116/1‑150 (MQ=255)
                                                             tgatgGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGt                                                                                          >  2:962045/1‑150 (MQ=255)
                                                              gatgGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTc                                                                                         <  1:1264091/150‑1 (MQ=255)
                                                              gatgGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTc                                                                                         <  1:1051176/150‑1 (MQ=255)
                                                                ggggcGCTGCGTCAAAATGTCGATGCTATCAGCATGGATGATCGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGAGCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTg                                                                                       <  1:2028721/146‑1 (MQ=255)
                                                                  gTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGCCCTGCAAAACTCACTGTCTGAc                                                                                     >  2:1435457/1‑150 (MQ=255)
                                                                  gTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGAc                                                                                     <  2:1761587/150‑1 (MQ=255)
                                                                          ggCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGcgacg                                                                             >  1:1709337/1‑150 (MQ=255)
                                                                          ggCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGcgacg                                                                             >  1:1209005/1‑150 (MQ=255)
                                                                           gCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGcgacga                                                                            >  2:878542/1‑150 (MQ=255)
                                                                                 tGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGc                                                                      <  2:877434/150‑1 (MQ=255)
                                                                                 tGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGc                                                                      >  1:616906/1‑150 (MQ=255)
                                                                                 tGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGc                                                                      <  1:1101412/150‑1 (MQ=255)
                                                                                        gCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGttt                                                               >  1:1202051/1‑150 (MQ=255)
                                                                                         cTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGTAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGCACGGCGCCGATTTTGCAGCGTTTg                                                              >  1:363813/1‑150 (MQ=255)
                                                                                         cTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTg                                                              >  2:1621188/1‑150 (MQ=255)
                                                                                          tATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTTATACAAGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTgc                                                             >  2:120238/1‑150 (MQ=255)
                                                                                          tATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAACACTCACTGTCTGACCGGCGACGAGTTTGCAGCGTTTgc                                                             >  1:342099/1‑150 (MQ=255)
                                                                                           aTCAGCATGTATGAACGGTGCTTAGAGTGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTgcg                                                            <  2:643782/150‑1 (MQ=255)
                                                                                            tCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGAAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTgcgc                                                           >  2:445131/1‑150 (MQ=255)
                                                                                             cAGCATGGATGAACGGGGCGTAGAGGGCATAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGTCCGGCGACGATTTTGCAGCGTTTGcgcg                                                          <  1:262878/150‑1 (MQ=255)
                                                                                             cAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGcgcg                                                          >  2:465852/1‑150 (MQ=255)
                                                                                             cAGAATGGATGAAAGGGGCGTAGAGGGCAAAAGTGTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGcgcg                                                          <  1:303080/150‑1 (MQ=255)
                                                                                              aGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGAAGCGTTTGCTCGt                                                         >  2:1095247/1‑150 (MQ=255)
                                                                                              aGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCAGCAAAACTCACTGTCTGACCGGCGAAGATTTTGCAGCGTTTTCGCGt                                                         >  2:1557725/1‑150 (MQ=255)
                                                                                                cATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGttt                                                       <  1:958691/150‑1 (MQ=255)
                                                                                                cATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCGGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCAGCAAAACTCACTGTCTGACCCGCGACGATTTTGCCGGGTTTGCGCGttt                                                       >  2:1570764/1‑150 (MQ=255)
                                                                                                   ggATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCaa                                                    <  2:252111/150‑1 (MQ=255)
                                                                                                      tGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCg                                                 <  1:438568/150‑1 (MQ=255)
                                                                                                       gAACGGGGCGTATAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAT‑‑‑GGGATAAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGCCGATTTTGCAGCGTTTGCGCGTTTCAACCGc                                                >  2:1282326/1‑150 (MQ=255)
                                                                                                        aaCGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACGCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCa                                               >  1:247128/1‑150 (MQ=255)
                                                                                                             ggCGTAGAGGGCAATAGGCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTc                                          <  2:1687038/150‑1 (MQ=255)
                                                                                                             ggCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTc                                          >  1:427087/1‑150 (MQ=255)
                                                                                                                 tAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACt                                      <  1:1887468/150‑1 (MQ=255)
                                                                                                                 tAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCGGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGAATTTGCAGCGTTTGCGCGTTTCAACCGCCCCGTCGagg                                      >  1:1968342/1‑148 (MQ=255)
                                                                                                                  agagGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACtt                                     <  1:772442/150‑1 (MQ=255)
                                                                                                                   gagGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACttt                                    >  1:561877/1‑150 (MQ=255)
                                                                                                                     gggCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCa                                  >  2:1705451/1‑150 (MQ=255)
                                                                                                                      ggCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCAc                                 >  1:1821174/1‑150 (MQ=255)
                                                                                                                       gCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCAcc                                <  2:1297783/150‑1 (MQ=255)
                                                                                                                       gCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCAcc                                <  1:1199612/150‑1 (MQ=255)
                                                                                                                        cAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAAACGGAATGAGGCTACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGAGGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCa                                  >  2:314454/1‑150 (MQ=255)
                                                                                                                           aactCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGc                            <  1:1638657/147‑1 (MQ=255)
                                                                                                                             gTCTGAAAAGAGAACCGGCCGGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCt                          <  1:1557726/150‑1 (MQ=255)
                                                                                                                                tGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCg                       >  1:905370/1‑150 (MQ=255)
                                                                                                                                  aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGc                     <  1:221151/150‑1 (MQ=255)
                                                                                                                                  aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGc                     >  2:1713145/1‑150 (MQ=255)
                                                                                                                                  aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGc                     <  1:176135/150‑1 (MQ=255)
                                                                                                                                  aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGACTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGc                     <  2:255684/150‑1 (MQ=255)
                                                                                                                                      gagaACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCCGCGTTTGCGCGTTTCAACCGCACCGTCGGCGTTCACCTGGCCGTCGGCAATc                 >  2:1301428/1‑150 (MQ=255)
                                                                                                                                      gagaACCGGCCTGTTGAAACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCCCCGTCGACTTTCACCTGGCCTTCGGCAATc                 >  2:858701/1‑150 (MQ=255)
                                                                                                                                        gaACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATcgc               >  2:249138/1‑150 (MQ=255)
                                                                                                                                        gaACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATcgc               >  1:383223/1‑150 (MQ=255)
                                                                                                                                         aaCCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCATCAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGAGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACGGTCGACTTTCACCTGGCCTTCGGCAATcgcg              >  1:1291287/1‑150 (MQ=255)
                                                                                                                                           ccGGCCTGTATCAACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAt            >  2:1188240/1‑150 (MQ=255)
                                                                                                                                           ccGGCCTGTATCAACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAt            >  1:538446/1‑150 (MQ=255)
                                                                                                                                             ggCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAttt          >  1:1091021/1‑150 (MQ=255)
                                                                                                                                             ggCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAttt          >  1:15137/1‑150 (MQ=255)
                                                                                                                                             ggCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAttt          <  1:1825180/150‑1 (MQ=255)
                                                                                                                                             ggCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAttt          >  2:1132865/1‑150 (MQ=255)
                                                                                                                                              gCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAtttt         >  1:1402250/1‑150 (MQ=255)
                                                                                                                                                 tGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGc      >  1:46725/1‑150 (MQ=255)
                                                                                                                                                  gTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGcc     <  2:453814/150‑1 (MQ=255)
                                                                                                                                                  gTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGcc     <  2:1311332/150‑1 (MQ=255)
                                                                                                                                                     gATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTgc  >  2:78646/1‑150 (MQ=255)
                                                                                                                                                     gATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTgc  >  1:451300/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
TTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAA‑‑‑GGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTGC  >  NC_000913/556492‑556790

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