Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,580 G→A 6.8% intergenic (+25/+216) 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,5800GA 6.8% 363.2 / 9.5 120intergenic (+25/+216)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base G (69/41);  new base A (6/2);  total (76/44)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

Only 100 of 121 total aligned reads displayed.

TTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATT  >  NC_000913/857431‑857729
                                                                                                                                                     |                                                                                                                                                     
ttGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCGGAACGCGTGATCGACTTTAGCGGTAATTACTAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGt                                                                                                                                                      >  2:888387/1‑150 (MQ=255)
  gTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTc                                                                                                                                                     >  1:1469172/1‑150 (MQ=255)
    ttCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGTAATGACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                                                                                   >  1:993994/1‑150 (MQ=255)
    ttCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATTACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                                                                                   >  1:1702021/1‑150 (MQ=255)
    ttCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                                                                                   >  2:626221/1‑150 (MQ=255)
     tCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGAGCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTCTCggg                                                                                                                                                  >  1:265180/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                                                                                                                                            <  2:1775469/150‑1 (MQ=255)
              cGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcggcg                                                                                                                                         <  2:938105/150‑1 (MQ=255)
                tGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACggcggc                                                                                                                                       >  1:1800144/1‑150 (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)
                        tAAGCTCCCGGCCCACCCGATTTCTGTAATACACCCCGCACCGGGTTTTCGACTTTAGGGTTAATTCCCAATATAACCTGCTTTTTAACGGTTTCGTTTAAGTTATTACCCGCCACCCCCAAACTGTCGGGGTGACGGCGGCTGATGACa                                                                                                                               <  2:606308/150‑1 (MQ=255)
                              cccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCa                                                                                                                         >  1:1106950/1‑150 (MQ=255)
                              cccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCa                                                                                                                         >  2:690690/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:1737912/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)
                                                          cccGGACCGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAg                                                                                             <  2:1774617/150‑1 (MQ=255)
                                                          cccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAg                                                                                             >  1:448047/1‑150 (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)
                                                                         cGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCt                                                                              >  2:1620776/1‑150 (MQ=255)
                                                                          gACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTg                                                                             >  1:586464/1‑150 (MQ=255)
                                                                            cTTTAGCGTTAATTACGAAGATTACCTGCGTAGTAATGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                                                                           <  2:12184/150‑1 (MQ=255)
                                                                                  cGGTAATTACAAAGATTACCAGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTAGGGGTAACGGCGGCTGATGACAAATTCAAAAAATCCTGATTCCATACGCTTAAAAGACCTACGCAGTACCAGAAaaatat                                                                     >  2:1152485/1‑150 (MQ=255)
                                                                                         tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCTTCACCCCGATAATGTCGGGGTGCCGGCGGCTGATGACAAATGCAAAATAGCATGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt                                                              >  2:14040/1‑150 (MQ=255)
                                                                                         tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt                                                              >  1:1325118/1‑150 (MQ=255)
                                                                                         tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACACCGATACTGTCGGGGTGACGGCTGCTGATGACAAAAGAAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt                                                              >  2:496301/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)
                                                                                            gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCa                                                           >  2:1051453/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)
                                                                                                    ccTGCGTAGAAAAGGGAGCTCGTAAGTTATAAGCCGTCCCCCCCATACTGGGGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTTCCTGCAATATATTGAATCTGCAGGACTTTg                                                   <  1:626221/150‑1 (MQ=25)
                                                                                                     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: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)
                                                                                                                   gATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGAGAAGGc                                    >  2:1327927/1‑150 (MQ=255)
                                                                                                                    aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCg                                   >  1:875517/1‑150 (MQ=255)
                                                                                                                       gAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTa                                >  1:1350301/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)
                                                                                                                               gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCa                        >  2:77476/1‑150 (MQ=255)
                                                                                                                               gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCa                        >  2:1588649/1‑150 (MQ=255)
                                                                                                                                aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAata                                                                        <  2:315364/101‑1 (MQ=255)
                                                                                                                                aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTACAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAAGCCGCAt                       >  1:1281750/1‑150 (MQ=255)
                                                                                                                                  tATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCa                     >  1:698423/1‑150 (MQ=255)
                                                                                                                                  tATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCa                     >  2:618638/1‑150 (MQ=255)
                                                                                                                                  tATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCa                     >  2:1511221/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)
                                                                                                                                    gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAgg                   >  1:112116/2‑150 (MQ=255)
                                                                                                                                     ccGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGc                  <  1:1135333/150‑1 (MQ=255)
                                                                                                                                     ccGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGc                  >  2:581377/1‑150 (MQ=255)
                                                                                                                                     ccGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGc                  >  1:481320/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)
                                                                                                                                          aCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTaa             >  1:1792357/1‑150 (MQ=255)
                                                                                                                                          aCCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTaa             >  2:1648176/1‑150 (MQ=255)
                                                                                                                                           ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTATGCCGTATAAGGCGTTAACGCCGCATCAGGCATTAAc            >  2:965198/1‑150 (MQ=255)
                                                                                                                                           ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc            >  1:813899/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)
                                                                                                                                                     |                                                                                                                                                     
TTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATT  >  NC_000913/857431‑857729

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: