Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,578 C→G 6.7% intergenic (+23/+218) ybiT → / ← ybiU ABC‑F family putative regulatory ATPase/DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,5780CG 6.7% 335.8 / 15.9 119intergenic (+23/+218)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base C (68/41);  new base G (6/2);  total (75/44)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

Only 100 of 122 total aligned reads displayed.

CTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCA  >  NC_000913/857429‑857727
                                                                                                                                                     |                                                                                                                                                     
cTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTTATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAc                                                                                                                                                       >  2:116409/1‑150 (MQ=255)
  ttGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCGGAACGCGTGATCGACTTTAGCGGTAATTACTAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGt                                                                                                                                                    >  2:888387/1‑150 (MQ=255)
    gTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTc                                                                                                                                                   >  1:1469172/1‑150 (MQ=255)
      ttCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGTAATGACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                                                                                 >  1:993994/1‑150 (MQ=255)
      ttCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                                                                                 >  2:626221/1‑150 (MQ=255)
        cacaCGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgggg                                                                                                                                               >  1:189749/1‑150 (MQ=255)
        cacaCGACCGTGACTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACGTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgggg                                                                                                                                               <  2:562715/150‑1 (MQ=255)
         acacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                              >  2:1217749/1‑150 (MQ=255)
         acacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                                                                              <  2:1878986/150‑1 (MQ=255)
          cacGACCGTGAGTTCGTAAGCTCCCTGGAGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGCCTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACACCGATACTGTCGGGGtt                                                                                                                                             >  1:98547/1‑149 (MQ=255)
           acGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGa                                                                                                                                            >  1:1184310/1‑150 (MQ=255)
             gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcg                                                                                                                                          <  1:1381417/150‑1 (MQ=255)
             gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcg                                                                                                                                          <  2:1775469/150‑1 (MQ=255)
                cGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcggcg                                                                                                                                       <  2:938105/150‑1 (MQ=255)
                   gAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCt                                                                                                                                    <  2:579683/150‑1 (MQ=255)
                    aGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtg                                                                                                                                   >  2:35711/1‑150 (MQ=255)
                     gTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtga                                                                                                                                  >  1:1827369/1‑150 (MQ=255)
                      ttCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgat                                                                                                                                 <  2:23632/150‑1 (MQ=255)
                         gTAAGCTCCCTGGCGACCCGCATTCTGGATATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGAc                                                                                                                              <  2:96896/150‑1 (MQ=255)
                          tAAGCTCCCTTGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACa                                                                                                                             <  1:1121720/150‑1 (MQ=255)
                                cccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCa                                                                                                                       >  2:690690/1‑150 (MQ=255)
                                cccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCa                                                                                                                       >  1:1106950/1‑150 (MQ=255)
                                cccTGGCGACCCGCATTCTGGAAAACACCCCGGAACGCGTGATCGACTTAAGCGGTAATTACGAATATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGGCGGGGTGACGGCGGCTGATGACAAATGCa                                                                                                                       <  1:774988/150‑1 (MQ=255)
                                 ccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                                                                      >  1:212646/1‑150 (MQ=255)
                                    ggCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACTAAGATTACCTGCGTAGTAAAGGGATCGAGAATGTTATTATCCGACACCCCGATACTGGCGGGGTGACGGCGGCTGATGACAAATGCAAAAt                                                                                                                   <  1:1896801/150‑1 (MQ=255)
                                         cccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTTGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGCCGGCGGCTGATGACAAATGAATAATAGCCt                                                                                                              >  2:80691/1‑150 (MQ=255)
                                             cATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATg                                                                                                          <  2:35492/150‑1 (MQ=255)
                                                 cTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCt                                                                                                      >  1:505309/1‑150 (MQ=255)
                                                 cTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCt                                                                                                      >  1:1737912/1‑150 (MQ=255)
                                                 cTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAAGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCt                                                                                                      >  1:709369/1‑150 (MQ=255)
                                                   ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTAc                                                                                                    >  1:1623917/1‑150 (MQ=255)
                                                      aaTCCCCCGGTAACGGGTTATCGACTTTAGGGTTAATTACCAATATTACTTGGTTTTTAATGGTATCGTTTATGTTTTTTTCGGTCACCCTTATACTGTGGGGGTGAGGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCt                                                                                                 <  2:1696535/150‑1 (MQ=255)
                                                      aaTCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGAGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCt                                                                                                 >  2:1671098/1‑150 (MQ=255)
                                                         cACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                              >  1:552756/1‑150 (MQ=255)
                                                         cACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                              >  1:341004/1‑150 (MQ=255)
                                                           ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCTAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCa                                                                                            <  2:545684/150‑1 (MQ=255)
                                                             ccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATACCCTGATGCCCAACCATTATCAgg                                                                                          >  1:1159920/1‑150 (MQ=255)
                                                               gcaaCGCGTTATCGACTTTCGCGTTAATTACTATTATTACCTGCGTAGTAACGGTATCGTTTAAGTTATTATCCGTCCCCCCTATACTGGCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGcc                                                                                        <  2:1838440/148‑1 (MQ=25)
                                                                        gATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTa                                                                               >  2:1463804/1‑150 (MQ=255)
                                                                          tCGACTTTAGCGGTATTTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTagc                                                                             <  1:1217749/150‑3 (MQ=255)
                                                                          tCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAcc                                                                             <  2:1532710/150‑1 (MQ=255)
                                                                            gACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTg                                                                           >  1:586464/1‑150 (MQ=255)
                                                                              cTTTAGCGTTAATTACGAAGATTACCTGCGTAGTAATGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                                                                         <  2:12184/150‑1 (MQ=255)
                                                                                ttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                                                                       <  1:974427/150‑1 (MQ=255)
                                                                                           tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCTTCACCCCGATAATGTCGGGGTGCCGGCGGCTGATGACAAATGCAAAATAGCATGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt                                                            >  2:14040/1‑150 (MQ=255)
                                                                                           tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACACCGATACTGTCGGGGTGACGGCTGCTGATGACAAAAGAAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt                                                            >  2:496301/1‑150 (MQ=255)
                                                                                            aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTg                                                           >  1:18506/1‑150 (MQ=255)
                                                                                             cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCAGCCCGTTACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc                                                          <  2:311954/150‑1 (MQ=255)
                                                                                             cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAATTGCAAAATAGCCTGCTGCGCTACGCTTATCAGGCCTACGCCGTACCTGCAATATATTGAATCTGc                                                          >  1:126245/1‑150 (MQ=255)
                                                                                             cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc                                                          <  2:1247080/150‑1 (MQ=255)
                                                                                             cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc                                                          >  1:1067907/1‑150 (MQ=255)
                                                                                               aaGATTACCTGCGTAGTAAAGGGATCGAGTATTTGATTATACGTCACCCCGATACTGTCGGTGTGACGTAGGCTGATGCGATATGCAAAATAGCATGATGCGCTAAGCTTATCAGGCCTACGAAGTACATGCAATTTATTGAATCTGCAg                                                        >  1:107097/1‑150 (MQ=255)
                                                                                                aGATTACGTGCGTAGTAAAGGGATCGAGTAAGGGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAgg                                                       <  2:285198/150‑1 (MQ=255)
                                                                                                      ccTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTg                                                 >  1:1663145/1‑150 (MQ=255)
                                                                                                       cTGCGTAGTAAACGGATCGAGTAAGTGATTATACGTCACCCCCATACTGTCGGGGTGACGGCGGCTGATGACAAAAGCAAAATACCCTGATGAGATACGCTTATCAGGCCTACGCAGTAACTCCAATATATTTAATCTCCAGGACTTTga                                                >  2:619127/1‑149 (MQ=38)
                                                                                                        tGCGTAGTAAAGGTTGCGAGTATGTGATTATCCGTCACCCCGATACTGGCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTa                                               <  1:1186737/150‑1 (MQ=38)
                                                                                                        tGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTa                                               >  1:45422/1‑150 (MQ=255)
                                                                                                          cGTTGTAAAGGGATCGAGTAAGTGTTTATCCGTCACCCAGATTCTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCGTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGATTCTGCAGGACTTTGTAgg                                             <  1:34195/150‑1 (MQ=38)
                                                                                                             agtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCg                                          >  2:1103234/1‑150 (MQ=255)
                                                                                                             agtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCg                                          >  2:772920/1‑150 (MQ=255)
                                                                                                               taAAGGGATCTATTAAGTGATTATCCTTCACCCCTATACTGTCGGGGTGAGGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGa                                        <  2:189749/150‑1 (MQ=37)
                                                                                                                 aaGTGATCTAGAAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATa                                      <  1:911827/150‑1 (MQ=255)
                                                                                                                    gtttcttgTAAGTGATTATCCGTCACCCGGTTAATGTCGGGTTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGGTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAgg                                   <  2:1699448/143‑1 (MQ=38)
                                                                                                                      aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCg                                 >  1:875517/1‑150 (MQ=255)
                                                                                                                         gAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTa                              >  1:1350301/1‑150 (MQ=255)
                                                                                                                           gTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGAGTTTGTAGGCCGGATAAGGCGTTAAc                            >  2:1610759/1‑150 (MQ=255)
                                                                                                                               ttGAATATCCGTCACCCCGACACTGTCGGGGTGACGGCGGCTGAGGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATAGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccg                        <  2:107872/149‑1 (MQ=255)
                                                                                                                               gTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccg                        <  1:690690/150‑1 (MQ=255)
                                                                                                                                tGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccgc                       <  2:1903513/150‑1 (MQ=255)
                                                                                                                                 gATTATCCGTCACCCCGATACTGTCGGGGTGCCGGCGGCTGATGACAAATGCAAAATAGCCTGATGCCCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGAGTTTGTAGGCCGGATAAGGCGTTAACCCCGCa                      >  2:1749618/1‑150 (MQ=255)
                                                                                                                                 gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCa                      >  2:77476/1‑150 (MQ=255)
                                                                                                                                 gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCa                      >  2:1588649/1‑150 (MQ=255)
                                                                                                                                  aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAata                                                                      <  2:315364/101‑1 (MQ=255)
                                                                                                                                  aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAata                                                                      >  1:315364/1‑101 (MQ=255)
                                                                                                                                  aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTACAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAAGCCGCAt                     >  1:1281750/1‑150 (MQ=255)
                                                                                                                                    tATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCa                   >  2:618638/1‑150 (MQ=255)
                                                                                                                                    tATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCa                   >  1:698423/1‑150 (MQ=255)
                                                                                                                                    tATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCa                   >  2:1518900/1‑150 (MQ=255)
                                                                                                                                     aTCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGAGAAGGCGTTAACGCCGCATCAg                  >  1:292186/1‑150 (MQ=255)
                                                                                                                                     aTCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACATGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGAAAAGGCGTTATCGCCtcatcat                  >  2:1797652/1‑149 (MQ=255)
                                                                                                                                      tCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAgg                 >  2:1248108/1‑150 (MQ=255)
                                                                                                                                      gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAgg                 >  1:112116/2‑150 (MQ=255)
                                                                                                                                       ccGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGc                <  1:1135333/150‑1 (MQ=255)
                                                                                                                                       ccGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGc                >  1:481320/1‑150 (MQ=255)
                                                                                                                                       ccGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGc                >  2:581377/1‑150 (MQ=255)
                                                                                                                                        cGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCa               >  2:925833/1‑150 (MQ=255)
                                                                                                                                        cGTCACCCCGATACTGTCGCGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCCCATCAGGCa               >  2:1810425/1‑150 (MQ=255)
                                                                                                                                         gTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAt              >  2:1362313/1‑150 (MQ=255)
                                                                                                                                          tCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAtt             <  2:1469172/150‑1 (MQ=255)
                                                                                                                                          ccACCCCCATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAtt             <  1:1758916/149‑1 (MQ=255)
                                                                                                                                           cACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAAGCTGCAGGACTTTGTAGGCCGGCTAAGGCGTTAACGCCGCATCAGGCATTa            >  1:1202183/1‑150 (MQ=255)
                                                                                                                                            aCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTaa           >  1:1792357/1‑150 (MQ=255)
                                                                                                                                            aCCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTaa           >  2:1648176/1‑150 (MQ=255)
                                                                                                                                             ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTATGCCGTATAAGGCGTTAACGCCGCATCAGGCATTAAc          >  2:965198/1‑150 (MQ=255)
                                                                                                                                             ccccGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc          >  2:1101884/1‑150 (MQ=255)
                                                                                                                                             ccccGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc          >  2:1109223/1‑150 (MQ=255)
                                                                                                                                              cccGATACTTTCGGTGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACa         <  2:993994/150‑1 (MQ=255)
                                                                                                                                               ccGATACTGTCGGTGTGACGGCGGCTGATGACAAATGCAAAATAGCATGATGCGCTACGCTTATCAGACCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGCTAACGCCGCATCAGGCATTAACaa        <  1:116409/150‑1 (MQ=255)
                                                                                                                                                 ttaaCTGTGGGGGTTACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAg      <  2:1414773/147‑1 (MQ=255)
                                                                                                                                                     cTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCa  >  1:1765522/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCA  >  NC_000913/857429‑857727

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: