Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,770,225 T→G 5.6% L92R (CTG→CGG)  PP_5600 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,770,2250TG 5.6% 575.6 / 16.9 178L92R (CTG→CGG) PP_5600hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (100/68);  new base G (4/6);  total (104/74)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 187 total aligned reads displayed.

TGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGGCAGC  >  NC_002947/3770076‑3770374
                                                                                                                                                     |                                                                                                                                                         
tgatgCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCt                                                                                                                                                           >  1:1527357/1‑150 (MQ=255)
   tgCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCt                                                                                                                                                        >  2:1747327/1‑150 (MQ=255)
    gccgccTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTc                                                                                                                                                       >  1:4879346/1‑150 (MQ=255)
       gccTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGcc                                                                                                                                                    >  1:2498994/1‑150 (MQ=255)
        ccTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCg                                                                                                                                                   <  1:1057149/150‑1 (MQ=255)
        ccTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCGGCTCGCCg                                                                                                                                                   <  1:3596146/150‑1 (MQ=255)
            gTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGa                                                                                                                                               >  1:1446637/1‑150 (MQ=255)
              gCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGc                                                                                                                                             >  2:5052127/1‑150 (MQ=255)
               cACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGcc                                                                                                                                            >  2:54223/1‑150 (MQ=255)
               cACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTGTCGGCCTGCTGACCCTGCTCGCCGGCGAGcc                                                                                                                                            >  2:54416/1‑150 (MQ=255)
                 cGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCgg                                                                                                                                          <  1:2828516/150‑1 (MQ=255)
                 cGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCAgg                                                                                                                                          <  2:3351433/150‑1 (MQ=255)
                 cGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCGGCTCGCCGGCGAGCAgg                                                                                                                                          <  2:4364706/150‑1 (MQ=255)
                  ggCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCggg                                                                                                                                         <  2:2649005/150‑1 (MQ=255)
                  ggCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCGGCTCGCCGGCGAGCAggg                                                                                                                                         <  2:5072204/150‑1 (MQ=255)
                      gACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGc                                                                                                                                    <  2:4779348/150‑1 (MQ=255)
                      gACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGc                                                                                                                                    <  1:2182617/150‑1 (MQ=255)
                            agcGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTa                                                                                                                           >  2:3037075/1‑150 (MQ=255)
                             gcGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTAc                                                                                                                          <  2:3475738/150‑1 (MQ=255)
                                      tCCATGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCcaaccaac                                                                                                                 <  1:1616940/150‑1 (MQ=255)
                                      tCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCcaaccaac                                                                                                                 <  1:1004256/150‑1 (MQ=255)
                                            gTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGc                                                                                                           <  1:4917988/150‑1 (MQ=255)
                                                gCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACg                                                                                                       >  2:3749583/1‑150 (MQ=255)
                                                 cAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTTCCCAACCAACCTCCGCCACGc                                                                                                      >  1:3941619/1‑150 (MQ=255)
                                                  aaGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGcc                                                                                                     >  2:2315236/1‑150 (MQ=255)
                                                   aGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTATTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCt                                                                                                    <  1:3281898/150‑1 (MQ=255)
                                                    gACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTg                                                                                                   >  2:2184036/1‑150 (MQ=255)
                                                    gACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCGGCTCGCCGGCGAGCAgggtc                                                                                                                                       <  1:3032395/118‑3 (MQ=255)
                                                      cTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCCGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGcc                                                                                                 >  1:4282260/1‑150 (MQ=255)
                                                       ttCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGccc                                                                                                >  2:1641078/1‑150 (MQ=255)
                                                       ttCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGccc                                                                                                >  2:3437385/1‑150 (MQ=255)
                                                       ttCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCGCCGCCACCCCGGccc                                                                                                >  1:621027/1‑150 (MQ=255)
                                                        tCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGc                                                                                                                                                     <  2:106017/100‑1 (MQ=255)
                                                         cAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTATTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCaa                                                                                              <  1:1850817/150‑1 (MQ=255)
                                                          cgCGAGCCGCCCCTGTCTGCGGGGGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGGTGGCCGTCTTTGGCGGCCTGCTGACCCTGCTCGCGGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAag                                                                                             <  1:1066198/149‑1 (MQ=255)
                                                          aGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAag                                                                                             <  1:3885549/150‑1 (MQ=255)
                                                          aGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAag                                                                                             >  2:5142197/1‑150 (MQ=255)
                                                           gCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAagc                                                                                            >  2:451509/1‑150 (MQ=255)
                                                            cGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAagca                                                                                           <  1:5393255/150‑1 (MQ=255)
                                                              aTCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAagcagc                                                                                         <  1:3539205/150‑1 (MQ=255)
                                                               tCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCCGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAgcg                                                                                        >  1:2773576/1‑150 (MQ=255)
                                                                 cGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCt                                                                                      >  2:461128/1‑150 (MQ=255)
                                                                 cGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCGGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCt                                                                                      <  1:1342985/150‑1 (MQ=255)
                                                                  gACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCCGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTg                                                                                     >  1:3710818/1‑150 (MQ=255)
                                                                    cccTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGc                                                                                   <  2:4949053/150‑1 (MQ=255)
                                                                     ccTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCGTACccaa                                                                                                                      >  1:178871/1‑114 (MQ=255)
                                                                      cTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCa                                                                                 >  2:1502796/1‑150 (MQ=255)
                                                                       tGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCAc                                                                                <  2:3721475/150‑1 (MQ=255)
                                                                         tCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGcc                                                                                                 <  1:4225381/131‑1 (MQ=255)
                                                                         tCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGcc                                                                                                 >  2:4225381/1‑131 (MQ=255)
                                                                             ccGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCa                                                                          >  2:3822077/1‑150 (MQ=255)
                                                                             ccGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGACCCGGGGAGC‑‑‑CCCGTACCCAACAAACCTCCGCCACGCGTGCCCACGCAGCGCTGGCCACGGGTCa                                                                          >  1:2065750/1‑150 (MQ=255)
                                                                                     gACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCcgcg                                                                 >  1:1565863/1‑150 (MQ=255)
                                                                                          cGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTTCCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCctc                                                             <  2:3941619/150‑1 (MQ=255)
                                                                                          cGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCctc                                                             >  2:2045574/1‑150 (MQ=255)
                                                                                           gTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCctct                                                            >  1:3546332/1‑150 (MQ=255)
                                                                                             gATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTcc                                                          >  2:1581089/1‑150 (MQ=255)
                                                                                             gATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTcc                                                          >  2:4122461/1‑150 (MQ=255)
                                                                                              aTCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTccc                                                         >  2:752616/1‑150 (MQ=255)
                                                                                               tCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACAGGTCAATCCCCGCGCCTCTCCcg                                                        >  1:5178057/1‑150 (MQ=255)
                                                                                               tCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCACCGCCACGCCAGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcg                                                        >  1:807440/1‑150 (MQ=255)
                                                                                                 gCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgct                                                      >  2:2016184/1‑150 (MQ=255)
                                                                                                  ccAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgctt                                                     >  1:4287618/1‑150 (MQ=255)
                                                                                                  ccAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgctt                                                     >  2:1599284/1‑150 (MQ=255)
                                                                                                   cAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttc                                                    <  2:4586874/150‑1 (MQ=255)
                                                                                                    aGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcg                                                   >  1:2672583/1‑150 (MQ=255)
                                                                                                     gCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgc                                                  <  1:2264972/150‑1 (MQ=255)
                                                                                                      ccTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgcg                                                 >  1:3562791/1‑149 (MQ=255)
                                                                                                       cTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgctt                                                >  1:1630782/1‑150 (MQ=255)
                                                                                                          ttCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTtctc                                             <  2:5106246/150‑1 (MQ=255)
                                                                                                          ttCAGCGCCACACTGCTGGCCGTATTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTtctc                                             <  1:3401716/150‑1 (MQ=255)
                                                                                                           tCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCACCCAACCTCCGCCACGCCTGCCCCAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTcg                                            >  2:1707761/1‑150 (MQ=255)
                                                                                                             aGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTCCCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTcgcg                                          >  2:653510/1‑150 (MQ=255)
                                                                                                             aGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCtctc                                                           >  1:3944381/1‑133 (MQ=255)
                                                                                                              gcgcCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCCATCCCCGCGCCTCTCCCGCTTCGCGTCTCGCGt                                         >  2:3628350/1‑150 (MQ=255)
                                                                                                                gcCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCcgc                                                                  <  2:3863569/123‑1 (MQ=255)
                                                                                                                 ccACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCACGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCGTCTCGCGTTCt                                      >  2:2472920/1‑150 (MQ=255)
                                                                                                                 ccACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCt                                      >  2:2472894/1‑150 (MQ=255)
                                                                                                                  cacaCTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTg                                     <  2:1446637/150‑1 (MQ=255)
                                                                                                                  cacaCTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTg                                     >  1:3846917/1‑150 (MQ=255)
                                                                                                                          tgGCCGTCTTTTTCGGCCTGCTTACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGcc                             >  1:4077741/1‑150 (MQ=255)
                                                                                                                          tgGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGcc                             >  1:4272784/1‑150 (MQ=255)
                                                                                                                            gCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCAc                           >  2:3247560/1‑150 (MQ=255)
                                                                                                                            gCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCAc                           >  2:1721413/1‑150 (MQ=255)
                                                                                                                            gCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCAc                           <  2:4879346/150‑1 (MQ=255)
                                                                                                                                 cTTTATCGGCCTGCTGACACTGCACGCCGGCGAGCCGGTCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTCGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTACTCGCGTTCTGACATCGCCACGAccc                      >  2:5332740/1‑150 (MQ=255)
                                                                                                                                     ttCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTg                  >  2:2964863/1‑150 (MQ=255)
                                                                                                                                      tCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTgc                 >  2:1847874/1‑150 (MQ=255)
                                                                                                                                          ccTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCGTCTCGCGTTCTGACATCGCCACGTCCCCGTgccgcc             >  2:2922133/1‑150 (MQ=255)
                                                                                                                                              gctgcccTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGag         >  1:4582315/5‑149 (MQ=255)
                                                                                                                                              gctgcccTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGAt         >  1:4582223/5‑150 (MQ=255)
                                                                                                                                              ctgACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCACCCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGAt         <  2:2498994/150‑1 (MQ=255)
                                                                                                                                               tgACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATc        <  2:1183899/150‑1 (MQ=255)
                                                                                                                                               ngACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGGTCTGACATCGCCACGTCCCCGTGCCGCCTGATc        >  1:2554881/2‑150 (MQ=255)
                                                                                                                                                 gcccGGCTCGCCGGCGAGCAGGGTCAGCAGGCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGa          >  2:329818/2‑150 (MQ=255)
                                                                                                                                                 cccgtGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCGTACCCAACCAACCTCCGCCACGCCTGCCCACGCAGCGCTGGCCACTGGTCACTCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCgg      <  1:4555514/146‑1 (MQ=255)
                                                                                                                                                    cTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGGCa    >  1:172409/1‑149 (MQ=255)
                                                                                                                                                     tGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCcgcg                                                                 <  2:4284559/87‑1 (MQ=255)
                                                                                                                                                     tGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCcgcg                                                                 <  2:498308/87‑1 (MQ=255)
                                                                                                                                                     tGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGGCAGc  >  2:5202164/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                         
TGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGC‑‑‑CCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGGCAGC  >  NC_002947/3770076‑3770374

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