Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,727,798 T→A 6.1% Y21* (TAT→TAA) ‡ 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,7980TA 6.1% 1074.6 / 39.6 328Y21* (TAT→TAA) ‡PP_3288universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base T (176/130);  new base A (4/16);  total (180/148)
Fisher's exact test for biased strand distribution p-value = 1.78e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 335 total aligned reads displayed.

GTAGGTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAAT  >  NC_002947/3727649‑3727946
                                                                                                                                                     |                                                                                                                                                    
tagtatGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTAt                                                                                                                                                      <  2:6796095/145‑1 (MQ=255)
gTAGGTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTAt                                                                                                                                                      <  1:2957264/150‑1 (MQ=255)
   ggTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGAc                                                                                                                                                   >  1:6729890/1‑150 (MQ=255)
    gTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGAcc                                                                                                                                                  >  1:1949888/1‑150 (MQ=255)
       ccTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCggg                                                                                                                                               <  1:6868580/150‑1 (MQ=255)
       ccTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCggg                                                                                                                                               >  2:1900615/1‑150 (MQ=255)
        cTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCgggg                                                                                                                                              <  2:3214586/150‑2 (MQ=255)
        cTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCgggg                                                                                                                                              <  1:6031956/150‑2 (MQ=255)
         ttGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCTGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGgct                                                                                                                                             >  2:1289013/1‑150 (MQ=255)
          tGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGgctg                                                                                                                                            >  1:1471644/1‑150 (MQ=255)
           gAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGgctgc                                                                                                                                           <  2:1147586/150‑1 (MQ=255)
             gTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCggggttcgg                                                                                                                                         <  2:5741404/150‑7 (MQ=255)
              ttGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCAGCAGCCCGGTCATAGACCggggttcggc                                                                                                                                        <  1:4940277/150‑8 (MQ=255)
                gATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTgcgc                                                                                                                                      >  1:5165777/1‑150 (MQ=255)
                gATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTgcgc                                                                                                                                      >  1:211468/1‑150 (MQ=255)
                 aTAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCt                                                                                                                                     <  1:3623734/150‑1 (MQ=255)
                 aTAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCt                                                                                                                                     <  1:4125561/150‑1 (MQ=255)
                 aTAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCt                                                                                                                                     <  1:5046068/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)
                     aCATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggcc                                                                                                                                 >  2:1692703/1‑150 (MQ=255)
                         cAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaagg                                                                                                                             <  2:4259553/150‑1 (MQ=255)
                            gTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaaggcca                                                                                                                          >  2:6258534/1‑150 (MQ=255)
                             ttCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaaggccaa                                                                                                                         <  1:4357587/150‑1 (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)
                                         aGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCt                                                                                                             <  1:5814197/150‑1 (MQ=255)
                                            tGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGca                                                                                                          >  1:2776447/1‑150 (MQ=255)
                                                aGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACAtc                                                                                                      >  2:2904332/1‑150 (MQ=255)
                                                    gCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACAtcgtcg                                                                                                  >  1:6071492/1‑150 (MQ=255)
                                                      tctACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCt                                                                                                >  1:1381096/1‑150 (MQ=255)
                                                       ctacctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCtt                                                                                               >  2:4758934/1‑150 (MQ=255)
                                                        tacctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCttt                                                                                              <  1:6459491/150‑1 (MQ=255)
                                                        tacctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCttt                                                                                              <  2:809674/150‑1 (MQ=255)
                                                        tacctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCttt                                                                                              <  1:1178147/150‑1 (MQ=255)
                                                         acctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTc                                                                                             <  1:1397065/150‑1 (MQ=255)
                                                         acctacGTATCGCAACCTCAGAGGTGTACCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTc                                                                                             <  1:857434/150‑1 (MQ=255)
                                                            tacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGAc                                                                                          >  2:757391/1‑150 (MQ=255)
                                                            tacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGAc                                                                                          >  2:3405080/1‑150 (MQ=255)
                                                             acGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCTCTTTCGACt                                                                                         >  1:6805794/1‑150 (MQ=255)
                                                                tATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACg                                                                                      >  1:6207031/1‑150 (MQ=255)
                                                                  tCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTg                                                                                    >  2:284637/1‑150 (MQ=255)
                                                                   cGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTgg                                                                                   <  2:3813983/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)
                                                                          tCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCt                                                                            >  2:3595398/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)
                                                                                     ccATGGAAACGCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGAGGCGCCGAACCCCGGTCTAGGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggcgg                                                                 <  2:5367447/150‑1 (MQ=255)
                                                                                       aTGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGc                                                               >  1:3094360/1‑150 (MQ=255)
                                                                                        tGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGcc                                                              >  2:3253979/1‑150 (MQ=255)
                                                                                        tGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGcc                                                              >  2:2285631/1‑150 (MQ=255)
                                                                                           aaaCTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTgg                                                           >  2:5479662/1‑150 (MQ=255)
                                                                                              cTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGa                                                        >  1:1006572/1‑150 (MQ=255)
                                                                                               tCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGaa                                                       <  2:5578266/150‑1 (MQ=255)
                                                                                                ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                      >  1:5400958/1‑150 (MQ=255)
                                                                                                ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                      <  2:6060784/150‑1 (MQ=255)
                                                                                                ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                      <  1:485451/150‑1 (MQ=255)
                                                                                                ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                      <  1:4755278/150‑1 (MQ=255)
                                                                                                ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                      <  1:3751694/150‑1 (MQ=255)
                                                                                                 cGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATg                                                     >  2:6771661/1‑150 (MQ=255)
                                                                                                  gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGa                                                    >  2:1379692/1‑150 (MQ=255)
                                                                                                  gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGa                                                    <  1:6650578/150‑1 (MQ=255)
                                                                                                  gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGa                                                    <  1:443899/150‑1 (MQ=255)
                                                                                                  gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGa                                                    >  1:2255564/1‑150 (MQ=255)
                                                                                                    aaaGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTCTGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGAcc                                                  >  1:1846966/1‑150 (MQ=255)
                                                                                                     aaGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCa                                                 <  1:3498796/150‑1 (MQ=255)
                                                                                                      agtctgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACCTCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAg                                                <  2:2539986/147‑1 (MQ=255)
                                                                                                      aGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAg                                                >  2:3713748/1‑150 (MQ=255)
                                                                                                          tgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGCCCAGgcgc                                            >  2:6611017/1‑150 (MQ=255)
                                                                                                          tgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGACCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgccc                                            <  1:3895677/150‑3 (MQ=255)
                                                                                                            ctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                          >  1:5478794/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                                      >  2:4642141/1‑150 (MQ=255)
                                                                                                                tgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAg                                      >  1:3282066/1‑150 (MQ=255)
                                                                                                                tgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCGGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAg                                      >  1:1717066/1‑150 (MQ=255)
                                                                                                                 gATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGt                                     >  2:3991569/1‑150 (MQ=255)
                                                                                                                 gATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGt                                     >  1:6616641/1‑150 (MQ=255)
                                                                                                                      ccTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATgc                                >  2:1620770/1‑150 (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:1307185/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)
                                                                                                                                  tGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACg                    <  1:4666990/150‑1 (MQ=255)
                                                                                                                                        gAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagc              <  1:2074769/150‑1 (MQ=255)
                                                                                                                                            cccGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagcagca          >  2:4895381/1‑150 (MQ=255)
                                                                                                                                            cccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagcagca          >  2:3025124/1‑150 (MQ=255)
                                                                                                                                                 tCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTg     <  2:6637150/150‑1 (MQ=255)
                                                                                                                                                    aTGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAAt  <  2:1074378/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTAGGTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAAT  >  NC_002947/3727649‑3727946

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: