Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,728,676 C→T 5.3% intergenic (+14/+16) PP_3288 → / ← PP_3289 universal stress protein family protein/acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,728,6760CT 5.3% 1025.9 / 28.8 321intergenic (+14/+16)PP_3288/PP_3289universal stress protein family protein/acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (150/153);  new base T (6/11);  total (156/165)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 323 total aligned reads displayed.

ATATCGACGTACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGACCATGCACTGAA  >  NC_002947/3728529‑3728825
                                                                                                                                                   |                                                                                                                                                     
 tatCGACGTACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAg                                                                                                                                                    <  1:844232/150‑1 (MQ=255)
  atCGACGTACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAgg                                                                                                                                                   >  2:5289292/1‑150 (MQ=255)
   tCGACGTACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCCTCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAggg                                                                                                                                                  >  1:5278992/1‑150 (MQ=255)
    cGACGTACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGc                                                                                                                                                 <  1:6672680/150‑1 (MQ=255)
     gACGTACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGcc                                                                                                                                                >  2:412516/1‑150 (MQ=255)
     gACGTACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGcc                                                                                                                                                >  1:4472478/1‑150 (MQ=255)
       cGTACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCtt                                                                                                                                              <  1:4877761/150‑1 (MQ=255)
       cGTACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCtt                                                                                                                                              <  1:1808942/150‑1 (MQ=255)
         tACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAATCACAAGGCCCTTGCGGGCCTTGt                                                                                                                                            <  2:2859562/150‑1 (MQ=255)
         tACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCTTGCGGGCCTTGt                                                                                                                                            <  2:3690004/150‑1 (MQ=255)
         tACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCTTGCGGGCCTTGt                                                                                                                                            <  1:224777/150‑1 (MQ=255)
         tACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGt                                                                                                                                            <  2:4879233/150‑1 (MQ=255)
         gaCTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCTTGCGGGCCTTGt                                                                                                                                            <  2:6320609/149‑1 (MQ=255)
          aCTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCTTGCGGGCCTTGtg                                                                                                                                           <  1:741069/150‑2 (MQ=255)
            tGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCtt                                                                                                                                         >  1:2141450/1‑150 (MQ=255)
            tGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCtt                                                                                                                                         <  2:4716773/150‑1 (MQ=255)
               tCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTc                                                                                                                                      >  2:3521577/1‑150 (MQ=255)
                 aTGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCaa                                                                                                                                    <  2:4296450/150‑1 (MQ=255)
                   gggCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAgg                                                                                                                                  >  1:2563630/1‑150 (MQ=255)
                     gCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCCTCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGcc                                                                                                                                >  2:4444209/1‑150 (MQ=255)
                     gCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGcc                                                                                                                                >  2:1011347/1‑150 (MQ=255)
                     gCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGcc                                                                                                                                >  2:1604055/1‑150 (MQ=255)
                      cAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCg                                                                                                                               <  1:6071474/150‑1 (MQ=255)
                        gTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCa                                                                                                                             <  2:1935938/150‑1 (MQ=255)
                        gTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAg                                                                                                                            <  1:5893271/150‑1 (MQ=255)
                         ttATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAg                                                                                                                            <  2:2563630/150‑1 (MQ=255)
                         ttATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAg                                                                                                                            <  1:1735410/150‑1 (MQ=255)
                         ttATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCACAAGGGCCTTGTCTTTTCAAGGCCGCAg                                                                                                                            >  1:3042480/1‑150 (MQ=255)
                            tCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTc                                                                                                                         <  2:1543464/150‑1 (MQ=255)
                              accaccGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCaa                                                                                                                       <  2:392391/150‑1 (MQ=255)
                               ccaccGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAg                                                                                                                      >  1:3946589/1‑150 (MQ=255)
                               ccaccGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAg                                                                                                                      >  1:4181585/1‑150 (MQ=255)
                                 accGTGACATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGtt                                                                                                                    >  1:6002239/1‑150 (MQ=255)
                                    gTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGt                                                                                                                 >  1:4789257/1‑150 (MQ=255)
                                    gTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGt                                                                                                                 <  1:2654630/150‑1 (MQ=255)
                                     tGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTg                                                                                                                <  1:6248224/150‑1 (MQ=255)
                                        cATCGGATTGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGATGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGAAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAg                                                                                                             <  1:242110/150‑1 (MQ=255)
                                         aTCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGa                                                                                                            >  1:4731000/1‑150 (MQ=255)
                                            ggCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAAtc                                                                                                         <  2:2694287/150‑1 (MQ=255)
                                                ggTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGc                                                                                                     <  1:3171974/150‑1 (MQ=255)
                                                  ttCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTc                                                                                                   <  2:2028200/150‑1 (MQ=255)
                                                     aTCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTcaca                                                                                                >  1:2167423/1‑150 (MQ=255)
                                                     aTCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTcaca                                                                                                >  2:5544039/1‑150 (MQ=255)
                                                      tCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTGGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACAt                                                                                               >  2:3274167/1‑150 (MQ=255)
                                                       cGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATc                                                                                              <  2:579024/150‑1 (MQ=255)
                                                            gTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATcgccgc                                                                                         <  2:1259410/150‑1 (MQ=255)
                                                              aCTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGccgccg                                                                                       <  1:6439611/150‑1 (MQ=255)
                                                               cTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTTTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCgg                                                                                      <  1:4268305/150‑1 (MQ=255)
                                                                tGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCggg                                                                                     >  2:2918515/1‑150 (MQ=255)
                                                                      gcgcGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCg                                                                               <  2:468668/150‑1 (MQ=255)
                                                                         cgTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGa                                                                            <  1:3438574/150‑1 (MQ=255)
                                                                          gTTATGCTTCGGTTGTCCAGGAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAATGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAg                                                                           <  1:3691790/150‑1 (MQ=255)
                                                                          gTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAg                                                                           >  1:4053184/1‑150 (MQ=255)
                                                                          gTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAg                                                                           <  1:3353107/150‑1 (MQ=255)
                                                                          gTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAg                                                                           <  2:661912/150‑1 (MQ=255)
                                                                          gTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAg                                                                           >  2:725781/1‑150 (MQ=255)
                                                                           ttATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAgg                                                                          <  1:5201688/150‑1 (MQ=255)
                                                                            tATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGt                                                                         <  1:1079444/150‑1 (MQ=255)
                                                                             aTGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTg                                                                        >  2:6730700/1‑150 (MQ=255)
                                                                                 aTCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGc                                                                    >  1:1253682/1‑150 (MQ=255)
                                                                                  tCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGcc                                                                   >  2:5143307/1‑150 (MQ=255)
                                                                                  tCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGcc                                                                   >  2:4229626/1‑150 (MQ=255)
                                                                                     gTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAg                                                                >  2:1228575/1‑150 (MQ=255)
                                                                                         tCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGcc                                                            <  2:5735276/150‑1 (MQ=255)
                                                                                           cagcagTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCtt                                                          >  1:6762132/1‑150 (MQ=255)
                                                                                              cagTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGAc                                                       <  1:5454411/150‑1 (MQ=255)
                                                                                                     cTGGTGATTTCTCCAGAGCGCTCCCACGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCgg                                                <  2:1379518/150‑1 (MQ=255)
                                                                                                     cTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCgg                                                >  1:2619637/1‑150 (MQ=255)
                                                                                                        gTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGc                                             >  2:3602206/1‑150 (MQ=255)
                                                                                                            tttCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATcg                                         >  2:4559855/1‑150 (MQ=255)
                                                                                                            tttCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATcg                                         >  2:568166/1‑150 (MQ=255)
                                                                                                              tctcCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATcgcg                                       <  2:686747/150‑1 (MQ=255)
                                                                                                                tcCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGcc                                                  >  2:2304733/1‑137 (MQ=255)
                                                                                                                tcCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTa                                     <  2:6305110/150‑1 (MQ=255)
                                                                                                                tcCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTa                                     >  1:3536308/1‑150 (MQ=255)
                                                                                                                   agagCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCCTCGCGTAGtt                                  >  2:6099257/1‑150 (MQ=255)
                                                                                                                   agagCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGtt                                  >  2:5997509/1‑150 (MQ=255)
                                                                                                                   agagCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGtt                                  <  1:1409388/150‑1 (MQ=255)
                                                                                                                      gcgcTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCt                               <  1:5398543/150‑1 (MQ=255)
                                                                                                                          tCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGcc                           >  1:2444700/1‑150 (MQ=255)
                                                                                                                            ccAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAg                         >  2:2652648/1‑150 (MQ=255)
                                                                                                                            ccAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAg                         >  1:3314975/1‑150 (MQ=255)
                                                                                                                             cAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAgg                        >  1:1171191/1‑150 (MQ=255)
                                                                                                                                   tGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTg                  <  1:1555713/150‑1 (MQ=255)
                                                                                                                                   tGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTg                  <  2:2167362/150‑1 (MQ=255)
                                                                                                                                    gAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTgg                 >  1:5131538/1‑150 (MQ=255)
                                                                                                                                     aaagacaAGGCCCTTGCGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGga                >  2:3030266/5‑150 (MQ=255)
                                                                                                                                     aaGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGga                <  2:1864270/150‑1 (MQ=255)
                                                                                                                                      aagacaAGGCCCTTGCGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGgag               >  2:1799525/4‑150 (MQ=255)
                                                                                                                                      aGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGgag               >  1:3765668/1‑150 (MQ=255)
                                                                                                                                      aGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGgag               <  1:1011347/150‑1 (MQ=255)
                                                                                                                                         acaAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGAcc            >  1:2803306/1‑150 (MQ=255)
                                                                                                                                         acaAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGAcc            >  2:2506612/1‑150 (MQ=255)
                                                                                                                                          caAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGACCa           <  1:193941/150‑1 (MQ=255)
                                                                                                                                          caAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGACCa           <  1:4398679/150‑1 (MQ=255)
                                                                                                                                           aaGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGACCAt          >  2:5415413/1‑150 (MQ=255)
                                                                                                                                            aGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGACCATg         <  2:4053184/150‑1 (MQ=255)
                                                                                                                                                ccGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGACCATGCACt     >  2:6574483/1‑150 (MQ=255)
                                                                                                                                                 cGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGTCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGACCATGCACTg    >  2:3176498/1‑150 (MQ=255)
                                                                                                                                                   cAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGACCATGCACTGaa  >  2:3330903/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATATCGACGTACTGGTCATGGGCAGTTATCACCACCGTGGCATCGGCTGGTTCATCGGCAGTACTGCCGAGCGCGTTATGCATCGGTTGTCCAGCAGTGTGCTGGTGATTTCTCCAGAGCGCTCCCAAGGCTGAAGCACAAGGCCCGCAAGGGCCTTGTCTTTTCAAGGCCGCAGGTCAAGTTCGTGCAGAATCTCGCTCACATCGCCGCCGGGGTAGCGGGAGGTGAAGCCCAGTGCCTTGACCAGCCGGTGCATCGCGTAGTTGCTGGCCAGGTCCCTGGAGACCATGCACTGAA  >  NC_002947/3728529‑3728825

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: