Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,729,090 T→A 11.8% L113Q (CTG→CAG)  PP_3292 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,729,0900TA11.8% 431.2 / 37.9 155L113Q (CTG→CAG) PP_3292membrane protein
Reads supporting (aligned to +/- strand):  ref base T (59/76);  new base A (7/11);  total (68/87)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 159 total aligned reads displayed.

GGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGC  >  NC_002947/3728942‑3729239
                                                                                                                                                     |                                                                                                                                                     
    gggCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCTCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTgg                                                                                                                                                  <  1:3820187/150‑1 (MQ=255)
    gggCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTgg                                                                                                                                                  <  2:4858579/150‑1 (MQ=255)
    gggCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTgg                                                                                                                                                  <  2:4858505/150‑1 (MQ=255)
     ggCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTggg                                                                                                                                                 <  2:414650/150‑1 (MQ=255)
       cGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGca                                                                                                                                               <  2:107740/150‑1 (MQ=255)
         gTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGTCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGcaca                                                                                                                                             <  2:1009558/150‑1 (MQ=255)
         gTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGcaca                                                                                                                                             >  2:1676624/1‑150 (MQ=255)
         gTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGcaca                                                                                                                                             <  1:2131948/150‑1 (MQ=255)
          tGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGccag                                                                                                                                             <  1:5318035/150‑4 (MQ=255)
          tGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCAGCAGGGccagg                                                                                                                                            <  2:246254/150‑5 (MQ=255)
           gccgCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAgg                                                                                                                                           >  2:342808/4‑150 (MQ=255)
           gCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGccagg                                                                                                                                            <  1:2377092/150‑5 (MQ=255)
           gCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAgg                                                                                                                                           <  1:5376329/150‑1 (MQ=255)
             aGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCa                                                                                                                                         >  1:2428453/1‑150 (MQ=255)
             aGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCa                                                                                                                                         >  2:3960725/1‑150 (MQ=255)
              gCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGCCAGGCAc                                                                                                                                        <  2:69487/150‑1 (MQ=255)
              gCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGCCAGGCAc                                                                                                                                        <  2:69500/150‑1 (MQ=255)
                  ggCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGCCAGGCAcaaca                                                                                                                                    <  2:4869124/150‑5 (MQ=255)
                  ggCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGCCAGGCAcaaca                                                                                                                                    <  2:4164075/150‑5 (MQ=255)
                    cGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTtg                                                                                                                                  <  1:3598609/150‑1 (MQ=255)
                     gTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTtgg                                                                                                                                 <  1:324179/150‑1 (MQ=255)
                     gTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTtgg                                                                                                                                 <  2:1137040/150‑1 (MQ=255)
                       gAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTtggct                                                                                                                               <  1:5021909/150‑1 (MQ=255)
                        aaGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTtggctg                                                                                                                              <  1:680980/150‑1 (MQ=255)
                         aGCGCACCGTTGGCAAAGCCATTTTCAACCGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTtggctgg                                                                                                                             <  2:3462829/150‑1 (MQ=255)
                          gcgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAgg                                                                                                                                           <  2:1261860/135‑1 (MQ=255)
                          gcgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAgg                                                                                                                                           >  1:1261860/1‑135 (MQ=255)
                           cgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCa                                                                                                                           <  1:2230867/150‑1 (MQ=255)
                           cgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCa                                                                                                                           <  2:93025/150‑1 (MQ=255)
                              aCCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATTGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAggg                                                                                                                        <  2:1455506/150‑1 (MQ=255)
                                cGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGgcg                                                                                                                      >  1:2881103/1‑150 (MQ=255)
                                       aaaGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTg                                                                                                               >  1:2634973/1‑150 (MQ=255)
                                         aGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGAc                                                                                                             <  2:4572362/150‑1 (MQ=255)
                                              ttttCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTtg                                                                                                                                  >  2:3818170/1‑124 (MQ=255)
                                               tttCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGAt                                                                                                       <  2:3910855/150‑1 (MQ=255)
                                                 tCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCg                                                                                                     >  2:631367/1‑150 (MQ=255)
                                                   aaCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGTGAATCGAACATCCTGCTGGTGGTTCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTg                                                                                                   <  1:1696340/150‑1 (MQ=255)
                                                   aaCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTg                                                                                                   >  2:2143407/1‑150 (MQ=255)
                                                    aCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTgg                                                                                                  <  1:956997/150‑1 (MQ=255)
                                                    aCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTgg                                                                                                  <  1:5109632/150‑1 (MQ=255)
                                                     cGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTATTGGCTGGCAGGGCGCTTAGTGACGCAGATCGTGGc                                                                                                 <  2:3716722/150‑1 (MQ=255)
                                                         ggaacgccAGTCTGGCCTGCCAGCCTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAAc                                                                                             >  1:4097973/8‑150 (MQ=255)
                                                         cAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAAc                                                                                             >  1:2218357/1‑150 (MQ=255)
                                                          aGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACt                                                                                            >  2:1450879/1‑150 (MQ=255)
                                                             agaCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCg                                                                                         <  1:3691667/150‑1 (MQ=255)
                                                              gaCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCATGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGa                                                                                        <  1:37970/150‑1 (MQ=255)
                                                              gaCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGa                                                                                        <  1:3857730/150‑1 (MQ=255)
                                                                 aGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCCGCTGGGCCCAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGc                                                                                     >  1:5201078/1‑150 (MQ=255)
                                                                     tGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAAc                                                                                 >  2:5233685/1‑150 (MQ=255)
                                                                                 cTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCaa                                                                     <  1:2101747/150‑1 (MQ=255)
                                                                                 cTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCaa                                                                     >  2:3105241/1‑150 (MQ=255)
                                                                                  tGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGTGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAg                                                                    <  1:330467/150‑1 (MQ=255)
                                                                                  tGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAg                                                                    <  2:2130249/150‑1 (MQ=255)
                                                                                   ggCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGt                                                                   <  1:2547633/150‑1 (MQ=255)
                                                                                    gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTg                                                                  >  1:5198611/1‑150 (MQ=255)
                                                                                    gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTg                                                                  <  1:1333330/150‑1 (MQ=255)
                                                                                     cGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTgg                                                                 <  1:1119440/150‑1 (MQ=255)
                                                                                      gTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGa                                                                <  2:3154105/150‑1 (MQ=255)
                                                                                      gTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGa                                                                <  1:2143407/150‑1 (MQ=255)
                                                                                       ttCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACAGCGAGGCCAACCTGCAGGCTCAAGTGGat                                                               >  1:1410356/1‑150 (MQ=255)
                                                                                          cTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGAtata                                                            >  1:4379166/1‑150 (MQ=255)
                                                                                          cTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCGTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACAGCGAGGCCAACCTGCAGGCTCAAGTGGAtata                                                            >  1:1923585/1‑150 (MQ=255)
                                                                                             cGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATAccc                                                         <  1:794382/150‑1 (MQ=255)
                                                                                                tCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGaa                                                      <  1:342808/150‑1 (MQ=255)
                                                                                                  aGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGa                                                    >  2:4040624/1‑150 (MQ=255)
                                                                                                       cTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaa                                               <  2:2103019/150‑1 (MQ=255)
                                                                                                       cTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaa                                               <  1:1807775/150‑1 (MQ=255)
                                                                                                       cTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCATTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaa                                               <  1:1046462/150‑1 (MQ=255)
                                                                                                        tCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGa                                                                                        <  2:1937808/108‑1 (MQ=255)
                                                                                                        tCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCCACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaaa                                              >  1:2382540/1‑150 (MQ=255)
                                                                                                        tCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCGTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGa                                                                                        >  1:1937808/1‑108 (MQ=255)
                                                                                                         cGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaaaa                                             <  1:5187591/150‑1 (MQ=255)
                                                                                                               tCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAACAACTCAcc                                       >  1:4811235/1‑150 (MQ=255)
                                                                                                                 gAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACcgc                                     >  2:3355628/1‑150 (MQ=255)
                                                                                                                  aaCATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCCACCTGCCGGCTCAAGTGGATATACCCGAAGATGAAAAAATCACcgcg                                    >  1:4250771/1‑150 (MQ=255)
                                                                                                                  aaCATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATAGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATAACCGAAGATGAAAAACTCACcgcg                                    >  1:4250565/1‑150 (MQ=255)
                                                                                                                   aCATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCgcgc                                   <  1:86639/150‑1 (MQ=255)
                                                                                                                          gctgGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCCACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCgg                            >  1:4159457/1‑150 (MQ=255)
                                                                                                                          gctgGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCgg                            <  2:130702/150‑1 (MQ=255)
                                                                                                                           ctgGTGGTGCTGTTGTGCCTGGCCCAGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGata                                                              >  1:591797/1‑115 (MQ=255)
                                                                                                                               tggtgCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAAc                       <  1:464903/150‑1 (MQ=255)
                                                                                                                                  tgCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCTGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCTAGTGGATATACCCGAAGTTGAAAAACTCACCGCGCTTCTCGGCAAACCgg                    >  2:3627078/1‑150 (MQ=255)
                                                                                                                                  tgCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCgg                    <  2:3751543/150‑1 (MQ=255)
                                                                                                                                   gCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCTAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgt                   >  2:4087750/1‑150 (MQ=255)
                                                                                                                                   gCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgt                   >  2:571726/1‑150 (MQ=255)
                                                                                                                                     tgttgtGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgtgt                 >  1:1534787/1‑150 (MQ=255)
                                                                                                                                       ttgtGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGc               >  2:2815843/1‑150 (MQ=255)
                                                                                                                                         gtgCCTGTGCCCAGCAGGGCCAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCAc             <  2:501630/150‑1 (MQ=255)
                                                                                                                                          tgCCTGTGCCCAGCAGGGCCAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACg            >  1:2410738/1‑150 (MQ=255)
                                                                                                                                          tgCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACg            >  1:1621604/1‑150 (MQ=255)
                                                                                                                                          tgCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACg            >  1:3218583/1‑150 (MQ=255)
                                                                                                                                            cctgtgCCCAGCAGGGCCAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCAc             >  2:684936/6‑147 (MQ=255)
                                                                                                                                            cctgtgCCCAGCAGGGCCAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTg          >  2:580974/6‑150 (MQ=255)
                                                                                                                                            cctgtgCCCAGCAGGGCCAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTg          >  2:2827190/6‑150 (MQ=255)
                                                                                                                                            ccTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTg          >  1:4217684/1‑150 (MQ=255)
                                                                                                                                             cTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGc         >  2:2435085/1‑150 (MQ=255)
                                                                                                                                                 gCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTg      >  2:651929/1‑150 (MQ=255)
                                                                                                                                                   ccTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGgc    >  2:2762046/1‑150 (MQ=255)
                                                                                                                                                   ccTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGgc    >  2:5333247/1‑150 (MQ=255)
                                                                                                                                                     tgctgGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGgcgc  >  1:2354254/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGC  >  NC_002947/3728942‑3729239

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: