Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,262,313 T→C 9.8% T158A (ACC→GCC) ‡ PP_2857 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,262,3130TC 9.8% 379.1 / 23.1 132T158A (ACC→GCC) ‡PP_2857hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (66/53);  new base C (6/7);  total (72/60)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 132 total aligned reads displayed.

TCCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTGGCA  >  NC_002947/3262165‑3262461
                                                                                                                                                    |                                                                                                                                                    
tcCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTa                                                                                                                                                     <  1:534129/150‑1 (MQ=255)
tcCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTa                                                                                                                                                     >  2:4421674/1‑150 (MQ=255)
  cGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGt                                                                                                                                                   >  2:27598/1‑150 (MQ=255)
    tACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCAACCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGa                                                                                                                                                 >  2:604822/1‑150 (MQ=255)
       gCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCt                                                                                                                                              <  2:29713/150‑1 (MQ=255)
       gCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGAcca                                                                                                                                              >  1:4786049/1‑149 (MQ=255)
         tATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCCGCTACCGCTGAAGGCCTGGACACTGGTAGTGACCTgg                                                                                                                                            >  2:3034387/1‑150 (MQ=255)
          aTTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGACCAGGc                                                                                                                                           <  1:4151946/150‑1 (MQ=255)
          aTTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCATCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGACCAGGc                                                                                                                                           <  2:1309616/150‑1 (MQ=255)
           ttAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGACCAGGcc                                                                                                                                          <  1:54830/150‑2 (MQ=255)
                ttAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGc                                                                                                                                     >  2:222196/1‑150 (MQ=255)
                ttAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTAACGGTGAAGGCCAGGTCACTGGTAGTGACCTTGCGATAGc                                                                                                                                     >  2:222116/1‑150 (MQ=255)
                ttAAGTGTATTTCTTAACTGCAGGGTCAATCCGCGAGACGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCACTCGGCTACCGGTGAAGGCCTGGTCACTTGTAGTTACCTGGCGATAGc                                                                                                                                     <  2:5388053/150‑1 (MQ=255)
                   aGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCctg                                                                                                                                  <  2:2368438/150‑1 (MQ=255)
                    gtgtCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCACGGCTACCGGTGAAGGCCAGGTCACTGGTAGTGACCTGGCGATAGCctgc                                                                                                                                 >  2:3537806/1‑150 (MQ=255)
                         gTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGaa                                                                                                                            >  1:5098419/1‑150 (MQ=255)
                          ttCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGTCTACCGGTGAAGGCCTGGTCACTACCAGTGACCAGGcc                                                                                                                                          <  2:2817821/135‑2 (MQ=255)
                                aCTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGt                                                                                                                     <  2:2593616/150‑1 (MQ=255)
                                  tGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTg                                                                                                                   >  1:4939309/1‑150 (MQ=255)
                                           aaTCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGa                                                                                                            <  1:556717/148‑1 (MQ=255)
                                           aaTCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgt                                                                                                          >  2:3503235/1‑150 (MQ=255)
                                           aaTCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgt                                                                                                          >  1:1909227/1‑150 (MQ=255)
                                             tacgcgACCATGATGATCAAGCTGATGTCCATTCGGTGTCCAACCACGGCACTATGGTTCAACTGGATCACCACTCGGCTACAGGAGTAGGCAAGGTCACTGGTAGTGACAAGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgatg                                                                                                        >  2:3611224/3‑150 (MQ=255)
                                             tCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGACTGTCAGTTGCTCCTGAgttg                                                                                                        >  2:3611225/1‑150 (MQ=255)
                                             tCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttg                                                                                                        >  2:3611222/1‑150 (MQ=255)
                                             tCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTAGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttg                                                                                                        >  2:4725031/1‑150 (MQ=255)
                                               cgcgACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttgtt                                                                                                      <  2:2511395/150‑1 (MQ=255)
                                               cgcgACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttgtt                                                                                                      <  2:722628/150‑1 (MQ=255)
                                                gcgACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGttgttg                                                                                                     <  2:28784/150‑1 (MQ=255)
                                                    cccGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCa                                                                                                 <  1:5039761/150‑1 (MQ=255)
                                                      cGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAgc                                                                                               <  1:2883219/150‑1 (MQ=255)
                                                           gCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCAc                                                                                          <  1:4707019/150‑1 (MQ=255)
                                                            cTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACg                                                                                         >  1:2076915/1‑150 (MQ=255)
                                                             tCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGt                                                                                        >  1:1253123/1‑150 (MQ=255)
                                                             tCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGt                                                                                        >  1:1253496/1‑150 (MQ=255)
                                                               aaGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTc                                                                                      >  2:230947/1‑150 (MQ=255)
                                                               aaGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTc                                                                                      >  1:3713758/1‑150 (MQ=255)
                                                                 ggTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTg                                                                                    <  2:4553066/150‑1 (MQ=255)
                                                                 ggTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTg                                                                                    >  1:706378/1‑150 (MQ=255)
                                                                  gtgtTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGt                                                                                   >  1:4223107/1‑150 (MQ=255)
                                                                  gtgtTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGt                                                                                   >  2:4934761/1‑150 (MQ=255)
                                                                    gtTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTg                                                                                 <  2:2910621/150‑1 (MQ=255)
                                                                    gtTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTg                                                                                 <  2:181268/150‑1 (MQ=255)
                                                                     ttGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGc                                                                                >  2:1518317/1‑150 (MQ=255)
                                                                       ggCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGt                                                                                                                                                   <  2:2817867/81‑1 (MQ=255)
                                                                       ggCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACAGGTGAAGGCCTGGTCACTGGTAGt                                                                                                                                                   >  1:2817867/1‑81 (MQ=255)
                                                                          cATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAg                                                                           <  2:1056342/150‑1 (MQ=255)
                                                                           aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttgt                                                                                                       <  1:633069/121‑1 (MQ=255)
                                                                           aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttgt                                                                                                       >  2:633069/1‑121 (MQ=255)
                                                                           aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAgg                                                                          <  2:1493838/150‑1 (MQ=255)
                                                                           aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAgg                                                                          >  2:4802980/1‑150 (MQ=255)
                                                                           aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAgg                                                                          >  2:1608386/1‑150 (MQ=255)
                                                                            ttCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAAt                                                                                                                           <  2:1902833/100‑1 (MQ=255)
                                                                            ttCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGa                                                                         >  1:3074784/1‑150 (MQ=255)
                                                                              cGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACg                                                                       <  2:4950050/150‑1 (MQ=255)
                                                                                     ccACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTgtcgt                                                                >  2:350400/1‑150 (MQ=255)
                                                                                       aCCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGtcgtcg                                                              >  2:144453/1‑150 (MQ=255)
                                                                                       aCCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGtcgtcg                                                              <  1:1599702/150‑1 (MQ=255)
                                                                                       aCCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTGCTGTTGCATGAGGACGCTGtcgtcg                                                              >  2:1894030/1‑150 (MQ=255)
                                                                                        ccccGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGa                                                             <  1:2711329/150‑1 (MQ=255)
                                                                                        ccccGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGa                                                             <  1:4718848/150‑1 (MQ=255)
                                                                                                ctctGGTTCAACTGGATCACCCCTCGGCTCCCGGTGAAGGCCTGGTCACTGGTCGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTGCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGc                                                     >  2:3751404/1‑150 (MQ=255)
                                                                                                 tctGGTTCAACTGTATCACCCCTCGGCTACCGGTGAAGGCCTGGTCTCTGTTAGTGACCTGGCGATATCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCttttct                                                    <  1:3194164/150‑4 (MQ=255)
                                                                                                 tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCa                                                                                                                       >  2:4000077/1‑83 (MQ=255)
                                                                                                 tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCa                                                                                                                       <  1:4000077/83‑1 (MQ=255)
                                                                                                 tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCg                                                    >  1:4732256/1‑150 (MQ=255)
                                                                                                 tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCg                                                    <  2:4654288/150‑1 (MQ=255)
                                                                                                 tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCg                                                    >  1:3363343/1‑150 (MQ=255)
                                                                                                 tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCg                                                    <  2:1936425/150‑1 (MQ=255)
                                                                                                   tGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGt                                                  >  1:4020932/1‑150 (MQ=255)
                                                                                                     gTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTtg                                                <  2:2126308/150‑1 (MQ=255)
                                                                                                         aaCTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGa                                                                            >  2:4502673/1‑118 (MQ=255)
                                                                                                         aaCTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGa                                                                            <  1:4502673/118‑1 (MQ=255)
                                                                                                         aaCTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCa                                            >  1:1055583/1‑150 (MQ=255)
                                                                                                                tCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATc                                                                                                  >  2:514199/1‑89 (MQ=255)
                                                                                                                tCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATc                                                                                                  <  1:514199/89‑1 (MQ=255)
                                                                                                                tCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGc                                     <  1:2121371/150‑1 (MQ=255)
                                                                                                                    cccTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGAt                                 <  2:5137101/150‑1 (MQ=255)
                                                                                                                         ggCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCa                            <  1:1456846/150‑1 (MQ=255)
                                                                                                                         ggCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCa                            <  2:4737841/150‑1 (MQ=255)
                                                                                                                             aCCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGc                        <  2:705823/150‑1 (MQ=255)
                                                                                                                               cGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGt                      >  1:1765456/1‑150 (MQ=255)
                                                                                                                               cGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGt                      >  1:2917450/1‑150 (MQ=255)
                                                                                                                               cGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGt                      >  2:600748/1‑150 (MQ=255)
                                                                                                                                 cataaGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTAGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTg                    >  2:3663189/4‑150 (MQ=255)
                                                                                                                                  tGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTgg                   <  2:2917450/150‑1 (MQ=255)
                                                                                                                                   gAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGc                  <  1:4421674/150‑1 (MQ=255)
                                                                                                                                       gCCAGGTCACTACCAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGc              >  2:3258483/1‑150 (MQ=255)
                                                                                                                                        ccaggTCACTACCAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCt             >  1:3645282/4‑150 (MQ=255)
                                                                                                                                        ccaggTCACTACCAGTGACCAGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGa                                                                            <  1:3265555/84‑1 (MQ=38)
                                                                                                                                        ccaggTCACTACCAGTGACCAGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGa                                                                            >  2:3265555/4‑87 (MQ=38)
                                                                                                                                         cTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCCGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTg            >  2:1700209/1‑150 (MQ=255)
                                                                                                                                         cTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTg            >  1:4936565/1‑150 (MQ=255)
                                                                                                                                         cTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTg            >  1:4025432/1‑150 (MQ=255)
                                                                                                                                             tCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGgtg        <  2:2718085/150‑1 (MQ=255)
                                                                                                                                               aCTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGtt      >  1:1890991/1‑150 (MQ=255)
                                                                                                                                                cTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTg     <  1:2245631/150‑1 (MQ=255)
                                                                                                                                                 tGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTgg    >  2:4040883/1‑150 (MQ=255)
                                                                                                                                                 tGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTgg    >  2:722427/1‑150 (MQ=255)
                                                                                                                                                  ggTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTGGc   <  1:222196/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTGGCA  >  NC_002947/3262165‑3262461

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