Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,571,650 C→T 5.9% D490D (GAC→GAT) ‡ treY → maltooligosyl trehalose synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,571,6500CT 5.9% 250.8 / 4.5 85D490D (GAC→GAT) ‡treYmaltooligosyl trehalose synthase
Reads supporting (aligned to +/- strand):  ref base C (42/38);  new base T (2/3);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

CCCGTGACCTGACCCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCAACTGACTGCGCCCAG  >  NC_002947/4571502‑4571790
                                                                                                                                                    |                                                                                                                                            
cccGTGACCTGACCCTGGGCGCGATTCGCCGGGCCTTGCAGGCGCTGGTGGCGGGTGGCCCGGTTTACCGAACCTACTTCAAGGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGAcc                                                                                                                                             >  1:290095/1‑150 (MQ=255)
 ccGTGACCTGACCCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGAc                                                                                                                                              >  1:439609/1‑148 (MQ=255)
 ccGTGACCTGACCCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGAc                                                                                                                                              <  2:439609/148‑1 (MQ=255)
  cGTGACCTGACCCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCtt                                                                                                                                           <  1:1320921/150‑1 (MQ=255)
    tgatctgaccGTGGGGGCGATGCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAACAGGCCCTGACCAATGCCCGGCAAGACCTTGc                                                                                                                                         <  2:107211/150‑1 (MQ=255)
    tgacctgaccCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGc                                                                                                                                         <  2:705816/150‑1 (MQ=255)
      acctgaccCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCg                                                                                                                                       >  2:860810/1‑150 (MQ=255)
       cctgaccCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGGTCTTGCCgg                                                                                                                                      >  1:732102/1‑149 (MQ=255)
       cctgaccCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGGTCTTGCCgg                                                                                                                                      <  2:532341/150‑2 (MQ=255)
        ctgaccCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGGTCTTGCCggg                                                                                                                                     <  2:732102/150‑3 (MQ=255)
        ctgaccCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGaa                                                                                                                                     >  1:805933/1‑150 (MQ=255)
        ctgaccCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGAAGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGaa                                                                                                                                     >  1:1858840/1‑150 (MQ=255)
         tgaccCTGGGCGCGATTCGCTGGGCCTTGCAGGCGTTGGTGGCTCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCGCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAg                                                                                                                                    <  2:201688/150‑1 (MQ=255)
          gaccCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGGTCTTGCCgggca                                                                                                                                   <  1:777352/150‑5 (MQ=255)
           accCTGGGCGCGAGTCGCCGGGCCGTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCt                                                                                                                                  <  2:290095/150‑1 (MQ=255)
            cccTGGGCGCGATTCGCAGGGCCTTGCAGGCTTTGGGTGCGCAGTACCCGGTGTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCATTATGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTg                                                                                                                                 >  2:1858785/1‑150 (MQ=255)
               tGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACt                                                                                                                              <  2:805933/150‑1 (MQ=255)
               ggggCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACt                                                                                                                              <  2:336117/149‑1 (MQ=255)
                gggCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTg                                                                                                                             <  2:195203/150‑1 (MQ=255)
                      gATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCt                                                                                                                       >  2:1271463/1‑150 (MQ=255)
                      gATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCt                                                                                                                       >  2:1090165/1‑150 (MQ=255)
                      gATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCt                                                                                                                       >  2:753746/1‑150 (MQ=255)
                          cGCCGGGCCTTGCAGGCGTTGGTGGCGGAGTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTg                                                                                                                   <  2:142876/150‑1 (MQ=255)
                            ccGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGa                                                                                                                 >  1:1480845/1‑150 (MQ=255)
                             cGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGAc                                                                                                                <  1:748839/150‑1 (MQ=255)
                                 ccTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAAc                                                                                                            <  1:1271463/150‑1 (MQ=255)
                                 ccTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAAc                                                                                                            <  1:1649836/150‑1 (MQ=255)
                                      cAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCACAGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGaa                                                                                                       <  2:268875/150‑1 (MQ=255)
                                       aGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCAGCTACTGGACCAACTTGacc                                                                                                      >  2:211598/1‑148 (MQ=255)
                                           gTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTg                                                                                                  >  1:1344102/1‑150 (MQ=255)
                                                tggCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTgg                                                                                             <  1:4719/150‑1 (MQ=255)
                                                  gcgcATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTgggc                                                                                           >  1:1235391/1‑150 (MQ=255)
                                                    gcATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCCGGGAGAAAAGTGCTTCCAGCAGGCCCTGACCAATGCCCGCCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTgggcgg                                                                                         >  2:1606667/1‑150 (MQ=255)
                                                    gcATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTgggcgg                                                                                         >  2:1573110/1‑150 (MQ=255)
                                                          cccGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCtt                                                                                                                                           >  1:1091404/1‑94 (MQ=255)
                                                          cccGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCtt                                                                                                                                           <  2:1091404/94‑1 (MQ=255)
                                                           ccGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGcc                                                                                  <  2:1235391/150‑1 (MQ=255)
                                                               tttACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTggc                                                                              <  1:578230/150‑1 (MQ=255)
                                                                    cGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCAt                                                                         <  2:616978/150‑1 (MQ=255)
                                                                     gAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATc                                                                        >  2:1722138/1‑150 (MQ=255)
                                                                     gAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATc                                                                        >  2:1722260/1‑150 (MQ=255)
                                                                      aaCCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCt                                                                       >  2:1398724/1‑150 (MQ=255)
                                                                      aaCCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCt                                                                       >  2:1398576/1‑150 (MQ=255)
                                                                         cTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGcc                                                                    >  2:1158980/1‑150 (MQ=255)
                                                                           aCTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGc                                                                                   >  2:596750/1‑133 (MQ=255)
                                                                           aCTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGc                                                                                   <  1:596750/133‑1 (MQ=255)
                                                                            cTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCAcc                                                                 >  2:1712744/1‑150 (MQ=255)
                                                                             ttCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCg                                                                >  1:1087073/1‑150 (MQ=255)
                                                                             ttCAATACCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCATGACCAATACCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGAACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCg                                                                >  1:673861/1‑150 (MQ=255)
                                                                                aaTGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGACAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCgggg                                                             <  1:583619/150‑1 (MQ=255)
                                                                                  tGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGGCGTCCCGGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGgcg                                                           >  2:66336/1‑150 (MQ=255)
                                                                                     cTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGCCCACCTTGAACCGTGGCTGGCCGGGTCAGACTGGCGGCATCCGCCACTGGGgtgcgc                                                        <  2:1858840/150‑1 (MQ=255)
                                                                                      tGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGcc                                                       >  2:1574354/1‑150 (MQ=255)
                                                                                           gcgcCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCgcaa                                                  <  2:1227499/150‑1 (MQ=255)
                                                                                           gcgcCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCgcaa                                                  <  1:1659054/150‑1 (MQ=255)
                                                                                           gcgcCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCgcaa                                                  <  1:830192/150‑1 (MQ=255)
                                                                                            cgcCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCgcaag                                                 <  1:1573110/150‑1 (MQ=255)
                                                                                               ccGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCgcaagcaa                                              >  1:1546924/1‑150 (MQ=255)
                                                                                                cGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAAc                                             <  2:203727/150‑1 (MQ=255)
                                                                                                 ggCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCCGGGCGCGCCCGCAAGCAACt                                            >  2:161159/1‑150 (MQ=255)
                                                                                                     cAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTgcgc                                        >  1:1736484/1‑150 (MQ=255)
                                                                                                      aGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGCCCCAATTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGcc                                       <  1:610244/150‑1 (MQ=255)
                                                                                                          cGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGc                                   >  2:1350385/1‑150 (MQ=255)
                                                                                                           gAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCgcaa                                                  >  2:1121899/1‑134 (MQ=255)
                                                                                                           gAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCgcaa                                                  <  1:1121899/134‑1 (MQ=255)
                                                                                                             aaCGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCtg                                >  1:407162/1‑150 (MQ=255)
                                                                                                                gTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCtgcgt                             >  2:210923/1‑150 (MQ=255)
                                                                                                                 ttcttcCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCtgcgtg                            >  1:214665/1‑150 (MQ=255)
                                                                                                                  tcttcCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCtgcgtgc                           <  1:1574354/150‑1 (MQ=255)
                                                                                                                   cttcCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGTGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCtgcgtgcg                          <  2:214665/150‑1 (MQ=255)
                                                                                                                    ttcCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTgcgc                         >  2:1448489/1‑150 (MQ=255)
                                                                                                                        agcagGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTcc                     >  2:661172/1‑150 (MQ=255)
                                                                                                                        agcagGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTcc                     <  1:210923/150‑1 (MQ=255)
                                                                                                                             gCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCt                                                                                 >  2:884982/1‑85 (MQ=255)
                                                                                                                             gCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCt                                                                                 <  1:884982/85‑1 (MQ=255)
                                                                                                                             gCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCaa                <  2:477074/150‑1 (MQ=255)
                                                                                                                              cccTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCCGCAac               >  2:210566/1‑150 (MQ=255)
                                                                                                                               ccTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCgcaa                                                  <  1:189637/114‑1 (MQ=255)
                                                                                                                               ccTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCgcaa                                                  >  2:189637/1‑114 (MQ=255)
                                                                                                                               ccTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCCGGGCGACCCCGCAAGCAACGACGCCACGCCGGCGTGCGCTTCCAGCAact              >  2:1153628/1‑150 (MQ=255)
                                                                                                                                 tGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAATCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCAactga            >  1:1178415/1‑150 (MQ=255)
                                                                                                                                   aCCAATGCCCGCCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCAactgact          <  1:1606667/150‑1 (MQ=255)
                                                                                                                                      aaTGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCAACTGACTgcg       >  1:689063/1‑150 (MQ=255)
                                                                                                                                       agcttcGGCAAGGTCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCAACTGACTgcgc      >  1:1512448/6‑150 (MQ=255)
                                                                                                                                        tGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCAACTGACTGCGcc     <  1:419191/150‑1 (MQ=255)
                                                                                                                                         gCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCAACTGACTGCGccc    <  2:1736484/150‑1 (MQ=255)
                                                                                                                                           ccGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCAACTGACTGCGCCCAg  >  2:1142607/1‑150 (MQ=255)
                                                                                                                                           ccGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCAACTGACTGCGCCCAg  <  2:330805/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CCCGTGACCTGACCCTGGGCGCGATTCGCCGGGCCTTGCAGGCGTTGGTGGCGCATTACCCGGTTTACCGAACCTACTTCAATGCCTGCGGGCGCCCGGCGCAGGACGAAACGTTCTTCCAGCAGGCCCTGACCAATGCCCGGCAAGACCTTGCCGAAGCTGACTGGCCGCTACTGGACCAACTTGAACGGTGGCTGGGCGGGCAAGCCTGGCGGCATCTGCCACCGGGGCGCGCCCGCAAGCAACTGCGCCACGCCTGCGTGCGCTTCCAGCAACTGACTGCGCCCAG  >  NC_002947/4571502‑4571790

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