Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,727,796 T→A 6.2% Y21N (TAT→AAT) ‡ PP_3288 → universal stress protein family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,727,7960TA 6.2% 1107.6 / 39.7 331Y21N (TAT→AAT) ‡PP_3288universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base T (173/136);  new base A (5/16);  total (179/152)
Fisher's exact test for biased strand distribution p-value = 5.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

GTAGGTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAA  >  NC_002947/3727649‑3727945
                                                                                                                                                   |                                                                                                                                                     
tagtatGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTAt                                                                                                                                                     <  2:6796095/145‑1 (MQ=255)
gTAGGTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTAt                                                                                                                                                     <  1:2957264/150‑1 (MQ=255)
   ggTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGAc                                                                                                                                                  <  2:2865232/150‑1 (MQ=255)
    gTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGAcc                                                                                                                                                 >  1:1949888/1‑150 (MQ=255)
       ccTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCggg                                                                                                                                              <  1:6868580/150‑1 (MQ=255)
         ttGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCTGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGgct                                                                                                                                            >  2:1289013/1‑150 (MQ=255)
          tGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGgctg                                                                                                                                           >  1:1471644/1‑150 (MQ=255)
             gTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCggggttcgg                                                                                                                                        <  2:5741404/150‑7 (MQ=255)
                gATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTgcgc                                                                                                                                     >  1:5165777/1‑150 (MQ=255)
                gATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTgcgc                                                                                                                                     >  1:211468/1‑150 (MQ=255)
                 aTAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCt                                                                                                                                    <  1:5046068/150‑1 (MQ=255)
                 aTAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCt                                                                                                                                    <  1:4125561/150‑1 (MQ=255)
                 aTAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCt                                                                                                                                    <  1:3623734/150‑1 (MQ=255)
                  tAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTg                                                                                                                                   >  1:1201810/1‑150 (MQ=255)
                   aaaCATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTgg                                                                                                                                  <  1:5413836/150‑1 (MQ=255)
                   aaaCATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCAGCAGCCCGGTCATAGACCGGGGTTCGgcgcatc                                                                                                                                  <  2:6433738/150‑4 (MQ=38)
                    aaCATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggc                                                                                                                                 <  1:6778280/150‑1 (MQ=255)
                       aTCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaa                                                                                                                              >  2:4666990/1‑150 (MQ=255)
                            gTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaaggcca                                                                                                                         >  2:6258534/1‑150 (MQ=255)
                               cGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGccaaggccaagg                                                                                                                      <  1:2162734/150‑1 (MQ=255)
                                gCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAggg                                                                                                                     <  2:5015419/150‑1 (MQ=255)
                                  aaCGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGaa                                                                                                                   <  1:3975672/150‑1 (MQ=255)
                                      tgtAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGgc                                                                                                               <  1:5859101/150‑1 (MQ=255)
                                           cTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGc                                                                                                          <  2:1785237/150‑1 (MQ=255)
                                            tGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGca                                                                                                         >  1:5978623/1‑150 (MQ=255)
                                            tGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGca                                                                                                         >  1:2776447/1‑150 (MQ=255)
                                                aGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACAtc                                                                                                     >  2:2904332/1‑150 (MQ=255)
                                                 gATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACAtcg                                                                                                    <  1:462957/150‑1 (MQ=255)
                                                      tctACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCt                                                                                               >  1:1381096/1‑150 (MQ=255)
                                                       ctacctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCtt                                                                                              >  2:4758934/1‑150 (MQ=255)
                                                        tacctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCttt                                                                                             <  1:1178147/150‑1 (MQ=255)
                                                        tacctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCttt                                                                                             <  2:809674/150‑1 (MQ=255)
                                                         acctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTc                                                                                            <  1:5301987/150‑1 (MQ=255)
                                                         acctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTc                                                                                            <  1:1397065/150‑1 (MQ=255)
                                                            tacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGAc                                                                                         >  2:757391/1‑150 (MQ=255)
                                                            tacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGAc                                                                                         >  2:3405080/1‑150 (MQ=255)
                                                                tATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACg                                                                                     >  1:6207031/1‑150 (MQ=255)
                                                                     cAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGaa                                                                                <  2:3912064/150‑1 (MQ=255)
                                                                         cTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGc                                                                            <  2:6866702/150‑1 (MQ=255)
                                                                         cTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGc                                                                            >  2:6545031/1‑150 (MQ=255)
                                                                          tCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCt                                                                           >  2:2662596/1‑150 (MQ=255)
                                                                            agagGTGTGCCATGGAAACTCCGCAAAGTCTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTGTCGACTCCGTGGAAGGGCTTg                                                                         >  1:2687691/1‑150 (MQ=255)
                                                                                gtgtGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAc                                                                     >  2:494456/1‑150 (MQ=255)
                                                                                   tgCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggc                                                                  <  1:2392662/150‑1 (MQ=255)
                                                                                     ccATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggcgg                                                                >  1:4910674/1‑150 (MQ=255)
                                                                                     ccATGGAAACGCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGAGGCGCCGAACCCCGGTCTAGGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggcgg                                                                <  2:5367447/150‑1 (MQ=255)
                                                                                       aTGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGc                                                              >  1:6060784/1‑150 (MQ=255)
                                                                                       aTGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGc                                                              >  1:4308990/1‑150 (MQ=255)
                                                                                       aTGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGc                                                              >  1:3094360/1‑150 (MQ=255)
                                                                                        tGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGTCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGcc                                                             >  2:1172794/1‑150 (MQ=255)
                                                                                        tGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGcc                                                             >  2:3253979/1‑150 (MQ=255)
                                                                                           aaaCTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTgg                                                          >  2:5479662/1‑150 (MQ=255)
                                                                                            aaCTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGt                                                         <  2:1466868/150‑1 (MQ=255)
                                                                                              cTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGa                                                       >  1:1006572/1‑150 (MQ=255)
                                                                                               tCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGaa                                                      <  2:5578266/150‑1 (MQ=255)
                                                                                               tCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCATCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGaa                                                      <  2:3222593/150‑1 (MQ=255)
                                                                                                ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                     <  1:485451/150‑1 (MQ=255)
                                                                                                ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                     <  2:6060784/150‑1 (MQ=255)
                                                                                                ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                     <  1:3878270/150‑1 (MQ=255)
                                                                                                ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                     <  1:3751694/150‑1 (MQ=255)
                                                                                                ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                     <  1:4755278/150‑1 (MQ=255)
                                                                                                 cGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATg                                                    >  2:6771661/1‑150 (MQ=255)
                                                                                                  gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGa                                                   <  1:6650578/150‑1 (MQ=255)
                                                                                                  gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGa                                                   >  1:2255564/1‑150 (MQ=255)
                                                                                                  gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGa                                                   >  1:784918/1‑150 (MQ=255)
                                                                                                  gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGa                                                   <  1:443899/150‑1 (MQ=255)
                                                                                                   cAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGAc                                                  <  1:3409308/150‑1 (MQ=255)
                                                                                                    aaaGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTCTGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGAcc                                                 >  1:1846966/1‑150 (MQ=255)
                                                                                                      agtctgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACCTCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAg                                               <  2:2539986/147‑1 (MQ=255)
                                                                                                      aGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAg                                               >  2:3713748/1‑150 (MQ=255)
                                                                                                          tgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGACTGGCCAGgcgc                                           >  1:6607790/1‑150 (MQ=255)
                                                                                                          tgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGCCCAGgcgc                                           >  2:6611017/1‑150 (MQ=255)
                                                                                                            ctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                         >  1:5478794/1‑150 (MQ=255)
                                                                                                            ctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                         >  1:6383771/1‑150 (MQ=255)
                                                                                                            ctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                         >  2:460314/1‑150 (MQ=255)
                                                                                                              gctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGc                                       >  2:2183141/1‑150 (MQ=255)
                                                                                                               ctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCa                                      >  1:2415995/1‑150 (MQ=255)
                                                                                                                tgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGa                                                   <  2:4153870/136‑1 (MQ=255)
                                                                                                                tgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAg                                     >  1:3282066/1‑150 (MQ=255)
                                                                                                                tgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCGGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAg                                     >  1:1717066/1‑150 (MQ=255)
                                                                                                                 gATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGt                                    >  2:3991569/1‑150 (MQ=255)
                                                                                                                      ccTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATgc                               >  2:1620770/1‑150 (MQ=255)
                                                                                                                      ccTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATgc                               <  2:1095222/150‑1 (MQ=255)
                                                                                                                       cTCCCCCCTGATGCGCCCAACCACGGTCAATGACCGGGCTGCAGACCTGGCAAAGGCCAAGGGAAAGGAGCTGCTCAACGTCGCTATCCACTACGAGGAAGGTCTAGCCACGGCGGGCCTGGTGAAAGACCAGGCGCTGGCAGTCATgcg                              >  2:5690448/1‑150 (MQ=255)
                                                                                                                           ccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAg                          <  2:6376866/150‑1 (MQ=255)
                                                                                                                           ccccTGATGAGCCGAACCCCGGTCTATGACCGTGATGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAg                          >  2:5176622/1‑150 (MQ=255)
                                                                                                                            cccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                         >  2:3602562/1‑150 (MQ=255)
                                                                                                                            cccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                         >  2:1307275/1‑150 (MQ=255)
                                                                                                                            cccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                         <  2:375803/150‑1 (MQ=255)
                                                                                                                            cccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                         <  2:4415342/150‑1 (MQ=255)
                                                                                                                                  tGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACg                   <  1:4591059/150‑1 (MQ=255)
                                                                                                                                    cgcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTg                 <  2:1908628/150‑1 (MQ=255)
                                                                                                                                     gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATg                                                                                                                 >  1:6714463/1‑53 (MQ=255)
                                                                                                                                        gAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagc             <  1:2074769/150‑1 (MQ=255)
                                                                                                                                            cccGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagcagca         >  2:4311329/1‑150 (MQ=255)
                                                                                                                                            cccGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagcagca         >  2:4895381/1‑150 (MQ=255)
                                                                                                                                            cccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagcagca         >  2:3025124/1‑150 (MQ=255)
                                                                                                                                               ggTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCg      <  1:5792511/150‑1 (MQ=255)
                                                                                                                                                 tCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTg    <  2:6637150/150‑1 (MQ=255)
                                                                                                                                                 tCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTg    <  2:2311673/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTAGGTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAA  >  NC_002947/3727649‑3727945

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