Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,140,698 C→A 5.5% M508I (ATG→ATT ptsA ← putative PTS multiphosphoryl transfer protein PtsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,140,6980CA 5.5% 300.0 / 4.8 90M508I (ATG→ATTptsAputative PTS multiphosphoryl transfer protein PtsA
Reads supporting (aligned to +/- strand):  ref base C (42/43);  new base A (2/3);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

CTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTGCCTCTGCGGGAAT  >  NC_000913/4140551‑4140847
                                                                                                                                                   |                                                                                                                                                     
ctctTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCAt                                                                                                                                                     >  1:1453924/1‑150 (MQ=255)
 tctTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCAtt                                                                                                                                                    >  1:25918/1‑150 (MQ=255)
 tctTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCAtt                                                                                                                                                    >  1:1210785/1‑150 (MQ=255)
 tctTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCCTCAAtg                                                                                                                                                    >  2:255927/1‑149 (MQ=255)
  ctTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTg                                                                                                                                                   <  2:831811/150‑1 (MQ=255)
  ctTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTg                                                                                                                                                   <  1:445195/150‑1 (MQ=255)
  ctTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTg                                                                                                                                                   <  1:477764/150‑1 (MQ=255)
     cGCAGCATTGATCGATGATGAACATCACCGACCGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTgagg                                                                                                                                                >  2:1253420/1‑150 (MQ=255)
       cagcaTTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTgaggag                                                                                                                                              >  1:1045074/1‑150 (MQ=255)
          caTTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGatc                                                                                                                                           <  2:1124014/150‑1 (MQ=255)
            ttGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCCTCAATGGAAGAGatctt                                                                                                                                         <  1:255927/150‑3 (MQ=255)
             tGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGatcatc                                                                                                                                        >  1:1330468/1‑150 (MQ=255)
              gATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCg                                                                                                                                       >  2:1372977/1‑150 (MQ=255)
              gATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCg                                                                                                                                       >  2:309169/1‑150 (MQ=255)
              gATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCg                                                                                                                                       >  2:477862/1‑150 (MQ=255)
               aTCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCCTCAATGGAAGAGatcttatg                                                                                                                                      <  2:561853/150‑6 (MQ=255)
                tCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATTTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCggg                                                                                                                                     <  2:1237534/150‑1 (MQ=255)
                 cGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCCTCAATGGAAGAGatcttatggg                                                                                                                                    <  2:638329/150‑8 (MQ=255)
                    tgatgaACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCa                                                                                                                                 >  1:769107/1‑150 (MQ=255)
                    tgatgaACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCa                                                                                                                                 >  1:472973/1‑150 (MQ=255)
                        gaACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGAt                                                                                                                             <  1:962576/150‑1 (MQ=255)
                              cACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTTCTCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAg                                                                                                                       <  2:315239/150‑1 (MQ=255)
                               aCCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTACATTGAGGAGATCATCGGGATCATGATTTTCAgg                                                                                                                      >  1:470075/1‑150 (MQ=255)
                                           tCCAGCATGATGCCGAGCTCGATTTTCTCATCAAACGCAATGTGTTCGTTACGTAGATGACGTTTGGCTTACGCCAGTTTTTCTTACACCCATAAGATCTCTTCCATTGAGGAGATCATCGGCATCATGATTTTCAGGCTGACGTGAGCg                                                                                                          >  2:472349/1‑150 (MQ=255)
                                               gCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAgg                                                                                                      <  1:1253420/150‑1 (MQ=255)
                                                   gatgCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGcgcg                                                                                                  >  2:107458/1‑150 (MQ=255)
                                                     tgCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCgga                                                                                                >  2:845975/1‑150 (MQ=255)
                                                     tgCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCgga                                                                                                <  2:1223753/150‑1 (MQ=255)
                                                           gCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCg                                                                                          >  1:532858/1‑150 (MQ=255)
                                                            cTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGa                                                                                         >  1:1511660/1‑150 (MQ=255)
                                                             tGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGAc                                                                                        >  1:162905/1‑150 (MQ=255)
                                                             tGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGATGCCGTGAGCGGAGGCGCGGAGGATCGAc                                                                                        >  1:1056978/1‑150 (MQ=255)
                                                              ggATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGAcc                                                                                       <  2:1045074/150‑1 (MQ=255)
                                                                aTTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGt                                                                                     <  2:969853/150‑1 (MQ=255)
                                                                  tttCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAg                                                                                   >  1:1414580/1‑150 (MQ=255)
                                                                  tttCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAg                                                                                   <  1:477862/150‑1 (MQ=255)
                                                                   ttCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGc                                                                                  >  2:242783/1‑150 (MQ=255)
                                                                   ttCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGc                                                                                  <  1:1185229/150‑1 (MQ=255)
                                                                    tctcATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCt                                                                                 >  2:55665/1‑150 (MQ=255)
                                                                       catcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCg                                                                              >  2:904933/1‑150 (MQ=255)
                                                                        atcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGt                                                                             <  1:1537120/150‑1 (MQ=255)
                                                                        atcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCACTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGt                                                                             >  1:912481/1‑150 (MQ=255)
                                                                         tcaAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTg                                                                            <  2:1220053/150‑1 (MQ=255)
                                                                          caAACGGAATGTGTTCGTTACGTAGGTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGGGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTgg                                                                           <  1:961897/150‑1 (MQ=255)
                                                                           aaaCGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGt                                                                          >  1:290792/1‑150 (MQ=255)
                                                                           aaaCGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGt                                                                          <  2:1407138/150‑1 (MQ=255)
                                                                            aaCGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAAGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTa                                                                         >  1:1288780/1‑150 (MQ=255)
                                                                               ggAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAaac                                                                      >  2:424075/1‑150 (MQ=255)
                                                                                  aTGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCCGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAaacaaa                                                                   <  1:460865/150‑2 (MQ=255)
                                                                                  aTGTGTTCGTTACGTAGTAGCTGTTTGGCTTCCGCCAGTATTTCTTTCACCCATAAGATCTCATCCATTCAGGAGATCATCGCGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAaacaac                                                                   <  1:460862/150‑1 (MQ=255)
                                                                                   tgtgTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAacg                                                                  >  2:622869/1‑150 (MQ=255)
                                                                                      gtTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAacgacg                                                               <  2:306300/150‑1 (MQ=255)
                                                                                      gtTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAacgacg                                                               >  1:581524/1‑150 (MQ=255)
                                                                                        tCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGAcgcg                                                             <  1:242783/150‑1 (MQ=255)
                                                                                        tCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGAcgcg                                                             >  1:1028963/1‑150 (MQ=255)
                                                                                              cGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTActct                                                       <  2:162905/150‑1 (MQ=255)
                                                                                               gTAGTTGCTGTTTGGCTTCCTCCAGTTTTTCTTTCACCCATTAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCtt                                                      <  1:1117274/150‑1 (MQ=255)
                                                                                                 aGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCa                                                    >  2:589423/1‑150 (MQ=255)
                                                                                                  gTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCAt                                                   <  1:981279/150‑1 (MQ=255)
                                                                                                     gCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATaaa                                                <  1:309169/150‑1 (MQ=255)
                                                                                                       tGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATa                                              >  1:798049/1‑150 (MQ=255)
                                                                                                       tGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATa                                              <  2:25918/150‑1 (MQ=255)
                                                                                                       tGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATa                                              <  2:172121/150‑1 (MQ=255)
                                                                                                        gTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATAc                                             >  1:338506/1‑150 (MQ=255)
                                                                                                        gTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATAc                                             <  2:318296/150‑1 (MQ=255)
                                                                                                        gTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATAc                                             <  1:450251/150‑1 (MQ=255)
                                                                                                          ttGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGc                                           <  1:1372977/150‑1 (MQ=255)
                                                                                                            ggCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCAc                                         >  1:311032/1‑150 (MQ=255)
                                                                                                            ggCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCAc                                         <  2:1509368/150‑1 (MQ=255)
                                                                                                              cTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACgg                                       <  2:1210785/150‑1 (MQ=255)
                                                                                                              cTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACgg                                       >  1:235836/1‑150 (MQ=255)
                                                                                                                tCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGgcg                                     >  2:898598/1‑150 (MQ=255)
                                                                                                                 ccGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCTTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGgcgc                                    >  2:435855/1‑150 (MQ=255)
                                                                                                                 ccGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGgcgc                                    <  2:770407/150‑1 (MQ=255)
                                                                                                                  cGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGcgcg                                   >  1:457516/1‑150 (MQ=255)
                                                                                                                  cGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGcgcg                                   >  2:1407901/1‑150 (MQ=255)
                                                                                                                        tttttCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAAcc                             <  2:769107/150‑1 (MQ=255)
                                                                                                                             ctttcACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGa                        >  1:1167050/1‑150 (MQ=255)
                                                                                                                                 cACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACgg                    >  1:728556/1‑150 (MQ=255)
                                                                                                                                  aCCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACggg                   <  1:1440939/150‑1 (MQ=255)
                                                                                                                                    ccATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGtt                 <  2:1271186/150‑1 (MQ=255)
                                                                                                                                     ccgatgATCTCCTCAATGGAAGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGttt                >  1:576608/6‑150 (MQ=255)
                                                                                                                                      aTAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCTTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTg               <  1:435855/150‑1 (MQ=255)
                                                                                                                                      aTAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTg               <  2:1349290/150‑1 (MQ=255)
                                                                                                                                      aTAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTg               <  1:1407901/150‑1 (MQ=255)
                                                                                                                                      aTAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTg               <  2:799421/150‑1 (MQ=255)
                                                                                                                                        aaGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTGcc             <  1:41314/150‑1 (MQ=255)
                                                                                                                                          gATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTGCctc           <  2:940559/150‑1 (MQ=255)
                                                                                                                                               ctTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTGCCTCTGCgg      <  1:588795/150‑1 (MQ=255)
                                                                                                                                               ctTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTGCCTCTGCgg      >  2:1512404/1‑150 (MQ=255)
                                                                                                                                                  ccATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTGCCTCTGCGGGaa   >  2:1201939/1‑150 (MQ=255)
                                                                                                                                                  ccATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTGCCTCTGCGGGaa   >  1:423081/1‑150 (MQ=255)
                                                                                                                                                   cATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTGCCTCTGCGGGAAt  >  2:873814/1‑150 (MQ=255)
                                                                                                                                                   cATTGACGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCCGAGGAAAGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTGCCTCTGCGGGAAt  >  2:594136/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCATGATTTTCAGGCTGCCGTGAGCGGAGGCGCGGAGGATCGACCGTAGCTGCGTGGTAAACAACGACGCGTACTCTTCATAAATACGCACGGCGCGATAACCGAGGAACGGGTTTGCCTCTGCGGGAAT  >  NC_000913/4140551‑4140847

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