Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,578 C→G 11.3% 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,5780CG11.3% 332.5 / 26.1 114intergenic (+23/+218)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base C (61/39);  new base G (5/8);  total (66/48)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 114 total aligned reads displayed.

TTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAA  >  NC_000913/857431‑857722
                                                                                                                                                   |                                                                                                                                                 
ttGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑Tg                                                                                                                                                <  1:1595021/150‑1 (MQ=255)
     tCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGCGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCggg                                                                                                                                           >  2:180958/1‑150 (MQ=255)
      cacaCGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCgggg                                                                                                                                          >  2:716473/1‑150 (MQ=255)
       acacGAGCGTGAGTGTGTAAGCTCCCTGGCGACCCTCATTCTGGAAAGCACCCCGTATCGCGTGATCGACTTTAGCGGTAATTAAGAAGATTAAATGCGTAGTAAAGTTATCGAGTAAGTGATTATCCGTCTCCCCGACAG‑TATCGGGGt                                                                                                                                         <  1:138280/150‑1 (MQ=255)
       acacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAGGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAG‑TATCGGGGt                                                                                                                                         <  1:415745/150‑1 (MQ=255)
        cacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAG‑TATCGGGGTg                                                                                                                                        <  1:435371/150‑1 (MQ=255)
        cacGACCGTGAGTTCGTAAGCTCCATGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGAAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAG‑TATCTGGGTg                                                                                                                                        <  2:1268360/150‑1 (MQ=255)
         acGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGa                                                                                                                                       >  1:1309290/1‑150 (MQ=255)
          cGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGAc                                                                                                                                      >  2:478478/1‑150 (MQ=255)
          cGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAG‑TATCGGGGTGAc                                                                                                                                      <  2:1253340/150‑1 (MQ=255)
          cGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGGCACCCCGATAC‑TGTCGGGGTGAc                                                                                                                                      <  1:1634604/150‑1 (MQ=255)
           gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAG‑TATCGGGGTGAcg                                                                                                                                     <  2:1690568/150‑1 (MQ=255)
            aCCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAG‑TATCGGGGTGAcgg                                                                                                                                    <  1:1228652/150‑1 (MQ=255)
             ccGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGAcggc                                                                                                                                   <  2:1663573/150‑1 (MQ=255)
             ccGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAG‑TATCGGGGTGAcgga                                                                                                                                   <  1:1273429/150‑2 (MQ=255)
              cGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGGAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGCTAC‑TGTCGGCGAGAcggcg                                                                                                                                  >  2:356203/1‑150 (MQ=255)
                tGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTATAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACggcggc                                                                                                                                <  2:890597/150‑1 (MQ=255)
                    ttCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGACCGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCtgat                                                                                                                            <  1:967816/150‑1 (MQ=255)
                    ttCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCtgat                                                                                                                            >  2:8554/1‑150 (MQ=255)
                      cGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCtgatga                                                                                                                          >  2:1029177/1‑150 (MQ=255)
                       gTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGAc                                                                                                                         >  2:422419/1‑150 (MQ=255)
                       gTAAGCTCCCTGGCGACCCGCAATCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGAc                                                                                                                         <  1:647102/150‑1 (MQ=255)
                           gCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAAt                                                                                                                     >  2:1209930/1‑150 (MQ=255)
                           gCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAAt                                                                                                                     >  2:300600/1‑150 (MQ=255)
                            cTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTTATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATg                                                                                                                    <  1:481160/150‑1 (MQ=255)
                            cTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCGCGCTAC‑TGTCGGGGTGACGGCGGCTGATGACAAATg                                                                                                                    >  1:1094319/1‑150 (MQ=255)
                             tCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGCCGGGGTGACGGCGGCTGATGACAAATGc                                                                                                                   <  2:1342279/150‑1 (MQ=255)
                              cccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGATGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCa                                                                                                                  >  1:941312/1‑150 (MQ=255)
                               ccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGTGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                                                                 <  2:892614/150‑1 (MQ=255)
                                       cccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                                                                         >  2:124593/1‑150 (MQ=255)
                                         cGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGa                                                                                                       >  1:1516472/1‑150 (MQ=255)
                                         cGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGACACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAAAGCAAAAAAGCCTGa                                                                                                       >  1:622055/1‑150 (MQ=255)
                                           cATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTATCTGCGTAGTAATGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATg                                                                                                     <  1:153333/150‑1 (MQ=255)
                                           cATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATg                                                                                                     >  1:1393927/1‑150 (MQ=255)
                                            aTTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgc                                                                                                    <  2:1190910/150‑1 (MQ=255)
                                            aTTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgc                                                                                                    <  1:910811/150‑1 (MQ=255)
                                                 ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTAc                                                                                               >  1:386628/1‑150 (MQ=255)
                                                   aaaTCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGc                                                                                             >  1:2834/1‑150 (MQ=255)
                                                    aaTCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCt                                                                                            <  1:422419/150‑1 (MQ=255)
                                                       cACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TTTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATTCGCTACGCTTAt                                                                                         >  2:600484/1‑150 (MQ=255)
                                                         ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTTATGACAAATGCAAAATAGCCTCATGCGCTACGCTTATCa                                                                                       >  2:1309453/1‑150 (MQ=255)
                                                           ccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAgg                                                                                     <  1:1066218/150‑1 (MQ=255)
                                                              caaCCCGTGATCGACTTTAGCGGTAATTACTAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGGGATTATCCGTCACCCCGATAC‑TGGCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCt                                                                                  <  2:1692509/149‑1 (MQ=255)
                                                                 cgcgTGATCGACTTTAGCGGTAATTACGATGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAG‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCATATCAGGCCTACg                                                                               >  1:157920/1‑150 (MQ=255)
                                                                      gATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTa                                                                          >  2:396672/1‑150 (MQ=255)
                                                                       aTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                         >  1:708271/1‑150 (MQ=255)
                                                                        tctaCTTTAGCGGTAATTACGACCATTACCTGCGTTGTAAAGGGATCGATAAAGTGATTATCCGCCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAcc                                                                        <  1:763993/147‑1 (MQ=38)
                                                                         cGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCt                                                                       >  1:159852/1‑150 (MQ=255)
                                                                          gACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTg                                                                      >  1:871341/1‑150 (MQ=255)
                                                                           aCTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGc                                                                     <  2:1124739/150‑1 (MQ=255)
                                                                            cTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                         >  2:1162923/1‑145 (MQ=255)
                                                                            cTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                         <  1:1162923/145‑1 (MQ=255)
                                                                            cTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                                                                    >  1:681481/1‑150 (MQ=255)
                                                                                aGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatat                                                                <  2:2834/150‑1 (MQ=255)
                                                                                     tAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa                                                           >  1:349339/1‑150 (MQ=255)
                                                                                      aaTTACGAAGATTACCTGCGTAGTTAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGaa                                                          >  2:1707462/1‑150 (MQ=255)
                                                                                          aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTg                                                      >  1:1146640/1‑150 (MQ=255)
                                                                                                   aCCTGCGTAGTAAAGGTATCGAGTAAGTGATTATGCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACttt                                             <  2:892382/150‑1 (MQ=255)
                                                                                                    ccTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTg                                            >  1:807299/1‑150 (MQ=255)
                                                                                                     cTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGt                                           <  2:708271/150‑1 (MQ=255)
                                                                                                       gCGTAGAAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCAGATGACAAATGCAAAATAGCCTCCTGCGCTACGCTTATCAGCCCTACGCAGTAACTGCAATAAATTGAATCTGCTGGACTTTGTAg                                         >  1:478064/1‑150 (MQ=255)
                                                                                                         gtagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGc                                       >  2:1686063/1‑150 (MQ=255)
                                                                                                           agtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCg                                     <  2:430147/150‑1 (MQ=255)
                                                                                                            gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGTTAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCgg                                    >  2:105960/1‑150 (MQ=255)
                                                                                                            gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGTTAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCgg                                    >  2:105962/1‑150 (MQ=255)
                                                                                                            gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCgg                                    <  1:908902/150‑1 (MQ=255)
                                                                                                            gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCgg                                    >  1:1450601/1‑150 (MQ=255)
                                                                                                            gtaAAGGGATCGAGTAAGTGATTATCCGTCACCACGATAC‑TGTCGGGGTGAAGGCGGCTGATGACAAAAGCAACATAGCCTGATTGGCTACGCTTATCAGGCCTACGCATTACATGCAAGATATTGTATCTGCAGATCGTTGTAGGCAgg                                    >  2:1604025/1‑150 (MQ=25)
                                                                                                            gtaAAGGGATCGAGTAAGTGATTAGCCGTCACCCCGATAC‑TGTAGGGGAGACGCAGGCTGATGACAAAAGCAAAATAGCCTGAGGCGCTACGCTAATCCGGCCTACGCCGTACCTGCACTATATATCAACTGCAGGCCTTTGGAGGCCgg                                    >  2:1597490/1‑150 (MQ=255)
                                                                                                                aGGGATCGAGTAAGTGTTTATCCGTCACCCCGATCC‑TGTCCGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATaa                                <  1:944144/150‑1 (MQ=255)
                                                                                                                 gggATCGAGTAATTTATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAg                               <  2:543902/150‑1 (MQ=255)
                                                                                                                 gggATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAg                               <  2:1001115/150‑1 (MQ=255)
                                                                                                                 gggATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGATGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAg                               <  1:808109/150‑1 (MQ=255)
                                                                                                                  ggATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGCCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAgg                              >  1:1147768/1‑150 (MQ=255)
                                                                                                                   gATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACATATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCATTATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                             >  2:1034593/1‑150 (MQ=255)
                                                                                                                   gATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                             >  2:42230/1‑150 (MQ=255)
                                                                                                                   gATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                             <  1:402419/150‑1 (MQ=255)
                                                                                                                       gAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTTAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTa                         >  2:209653/1‑150 (MQ=255)
                                                                                                                       gAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTa                         >  1:738242/1‑150 (MQ=255)
                                                                                                                          tAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcg                      <  2:1393927/150‑1 (MQ=255)
                                                                                                                            tgttaTTATCCGTCACCCC‑TTACGTGTGGGGTTTACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgcc                    <  2:807299/146‑1 (MQ=21)
                                                                                                                            attGATTATCCGTCACCCCGATAA‑TGTGGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgcc                    <  2:1467651/148‑1 (MQ=255)
                                                                                                                            aGTGATTATCCGTCACCCCGATAC‑TGTCGTGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgcc                    >  2:1242231/1‑150 (MQ=255)
                                                                                                                            aGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAGTGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgcc                    >  2:1432666/1‑150 (MQ=255)
                                                                                                                               gATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCa                 >  1:833121/1‑150 (MQ=255)
                                                                                                                                 gccgccGTCACCCCGACAG‑TATCGGGGTGCCGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCGACGCAGTACCGGCAATTTATTGAATCTGCAGGACTTTGTAGGGCGGATAAGGCGTTTACGCCGCAtg               >  1:823463/5‑149 (MQ=21)
                                                                                                                                   cgccGTCACCCCGACAG‑TATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAg             >  2:1461708/3‑150 (MQ=255)
                                                                                                                                    tCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATTTATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAgg            >  2:1714262/1‑150 (MQ=255)
                                                                                                                                    tCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATATGGCGTTAACGCCGCATCAgg            >  2:195441/1‑150 (MQ=255)
                                                                                                                                    tCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAgg            >  2:305033/1‑150 (MQ=255)
                                                                                                                                        ccACCCCTAAAC‑TGAGGGGTTCACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAtt        <  2:1273279/149‑1 (MQ=21)
                                                                                                                                         cACCGCGATAC‑TGGCGGGGTGAGGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCTGATGCGGCGTTAACGCCTTATCCGGcctac       <  2:1378051/150‑5 (MQ=21)
                                                                                                                                         cACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa       >  1:1056657/1‑150 (MQ=255)
                                                                                                                                         cACCCCGACAG‑TATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa       >  1:1631551/1‑150 (MQ=255)
                                                                                                                                           ccccGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc     >  1:680707/1‑150 (MQ=255)
                                                                                                                                           ccccGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCCCATCAGGCATTAAc     >  2:45325/1‑150 (MQ=255)
                                                                                                                                           ccccGACAG‑TATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc     >  1:1502913/1‑150 (MQ=255)
                                                                                                                                            cccGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATATGGCGTTAACGCCGCATCAGGCATTAACa    <  1:195441/150‑1 (MQ=255)
                                                                                                                                              cGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACaaa  <  1:1235677/150‑1 (MQ=255)
                                                                                                                                              cGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACaaa  >  1:940220/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAA  >  NC_000913/857431‑857722

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: ‑