Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,576 T→C 10.7% intergenic (+21/+220) 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,5760TC10.7% 287.8 / 23.1 114intergenic (+21/+220)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base T (62/39);  new base C (4/8);  total (66/48)
Fisher's exact test for biased strand distribution p-value = 7.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-02

Only 100 of 116 total aligned reads displayed.

ATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAA  >  NC_000913/857427‑857722
                                                                                                                                                     |                                                                                                                                                   
aTCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGAt                                                                                                                                                     <  1:587544/150‑1 (MQ=255)
aTCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGAt                                                                                                                                                     <  2:612100/150‑1 (MQ=255)
    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)
                                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)
                                              gCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGAt                                                                                                      >  2:402419/1‑150 (MQ=255)
                                               cATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTATCTGCGTAGTAATGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATg                                                                                                     <  1:153333/150‑1 (MQ=255)
                                               cATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATg                                                                                                     >  1:1393927/1‑150 (MQ=255)
                                                aTTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgc                                                                                                    <  1:910811/150‑1 (MQ=255)
                                                aTTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgc                                                                                                    <  2:1190910/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‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCa                                                                                       >  2:1662173/1‑150 (MQ=255)
                                                               ccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAgg                                                                                     >  2:1235677/1‑150 (MQ=255)
                                                               ccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAgg                                                                                     <  1:1066218/150‑1 (MQ=255)
                                                                   aaCGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTa                                                                                 <  2:64242/150‑1 (MQ=255)
                                                                     cgcgTGATCGACTTTAGCGGTAATTACGATGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAG‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCATATCAGGCCTACg                                                                               >  1:157920/1‑150 (MQ=255)
                                                                     cgcgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGCTAC‑TGGAGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCAACGCTTATCAGGCCTACg                                                                               >  1:152509/1‑150 (MQ=255)
                                                                          gATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTa                                                                          >  2:396672/1‑150 (MQ=255)
                                                                           aTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                         >  1:708271/1‑150 (MQ=255)
                                                                             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                                                                         <  1:1162923/145‑1 (MQ=255)
                                                                                cTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                         >  2:1162923/1‑145 (MQ=255)
                                                                                cTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                                                                    >  1:681481/1‑150 (MQ=255)
                                                                                    aGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatat                                                                <  2:2834/150‑1 (MQ=255)
                                                                                     gCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatata                                                               <  1:124593/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)
                                                                                                           gCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGATGATGACAAATTCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACGTTGGAg                                         >  1:471537/1‑150 (MQ=255)
                                                                                                           gCGTAGAAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCAGATGACAAATGCAAAATAGCCTCCTGCGCTACGCTTATCAGCCCTACGCAGTAACTGCAATAAATTGAATCTGCTGGACTTTGTAg                                         >  1:478064/1‑150 (MQ=255)
                                                                                                             gtagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGc                                       >  2:1686063/1‑150 (MQ=255)
                                                                                                                gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGTTAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCgg                                    >  2:105962/1‑150 (MQ=255)
                                                                                                                gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGTTAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCgg                                    >  2:105960/1‑150 (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)
                                                                                                                 taAAGGGATCTATTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGTGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGa                                   <  1:901645/150‑1 (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‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAgg                              >  2:18369/1‑150 (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)
                                                                                                                         tCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGt                           <  2:1377591/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)
                                                                                                                                aGTGATTATCCGTCACCCCGATAC‑TGTCGTGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgcc                    >  2:1242231/1‑150 (MQ=255)
                                                                                                                                aGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAGTGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgcc                    >  2:1432666/1‑150 (MQ=255)
                                                                                                                                   gATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCa                 >  1:833121/1‑150 (MQ=255)
                                                                                                                                    aTTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCAGGATAAGGCGTTAACGCCGAAt                >  1:224826/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)
                                                                                                                                            ccACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAtt        <  1:716473/149‑1 (MQ=255)
                                                                                                                                             cACCGCGATAC‑TGGCGGGGTGAGGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCTGATGCGGCGTTAACGCCTTATCCGGcctac       <  2:1378051/150‑5 (MQ=21)
                                                                                                                                             cACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa       >  1:1056657/1‑150 (MQ=255)
                                                                                                                                               ccccGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc     >  1:680707/1‑150 (MQ=255)
                                                                                                                                               ccccGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCCCATCAGGCATTAAc     >  2:45325/1‑150 (MQ=255)
                                                                                                                                                cccGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATATGGCGTTAACGCCGCATCAGGCATTAACa    <  1:195441/150‑1 (MQ=255)
                                                                                                                                                  cGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACaaa  >  1:940220/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAC‑TGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAA  >  NC_000913/857427‑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: ‑