Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,432 T→G 13.9% L159R (CTG→CGG)  PP_3661 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,160,4320TG13.9% 248.7 / 7.1 145L159R (CTG→CGG) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base T (57/62);  new base G (17/7);  total (76/69)
Fisher's exact test for biased strand distribution p-value = 4.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

CGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCC  >  NC_002947/4160283‑4160581
                                                                                                                                                     |                                                                                                                                                     
cGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCt                                                                                                                                                       <  1:2935076/150‑1 (MQ=255)
   gtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATcggcgcggcg                                                                                                                                                    >  2:1940331/1‑150 (MQ=255)
    tCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGgcgc                                                                                                                                                   >  1:2318388/1‑150 (MQ=255)
     cGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCa                                                                                                                                                  <  2:3563713/150‑1 (MQ=255)
       gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                                <  2:3869061/150‑1 (MQ=255)
          cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTgcgccgcgcc                                                                                                                                             >  1:1072276/1‑150 (MQ=255)
             gttgTCACCAGTACCCTGCTGTTCGGCCTGTGGTGTGGCATGAAAGTTTTCAAGCTGTACCGCGATACCGCCTTGCTGACCGCCGCGGGATGCGCAATCTGCGGCGCCGCCGCCGTGGTGGCATTCGAATCGGCGCTGCGCTGCGccccc                                                                                                                                          <  2:2130642/147‑1 (MQ=255)
             gTAGTCACCAGTACCCTGTTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccgat                                                                                                                                          <  2:2644177/150‑4 (MQ=255)
             gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGccccc                                                                                                                                          >  1:2606621/1‑150 (MQ=255)
                 tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCaga                                                                                                                                      >  2:1685375/1‑148 (MQ=255)
                 tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCaaa                                                                                                                                      >  2:575598/1‑148 (MQ=255)
                   aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGcccccccaaaa                                                                                                                                    >  2:603196/1‑145 (MQ=255)
                   aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCCCGcccccccgaaa                                                                                                                                    >  1:3867375/1‑145 (MQ=255)
                   aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAaaaa                                                                                                                                    >  1:4014501/1‑150 (MQ=255)
                   aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAaaaa                                                                                                                                    >  1:360120/1‑150 (MQ=255)
                   aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGcccccccaaaa                                                                                                                                    >  1:3213771/1‑145 (MQ=255)
                   aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATCCGAATCGGCGCGGCGCCCCGcccccccccaa                                                                                                                                    >  1:2683293/1‑145 (MQ=255)
                    ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                                   <  2:2536380/150‑1 (MQ=255)
                    ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCGGCCGCCGTGCTGGCATTGGAATCGGCGCTGCGCAGGGCCCCCCATaaaa                                                                                                                                   <  2:2681179/150‑1 (MQ=255)
                     cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATCAAAg                                                                                                                                  >  2:1607799/1‑150 (MQ=255)
                             tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCCCCGCCCCCCCAAAAAGcccccagg                                                                                                                          >  1:1615413/1‑146 (MQ=255)
                                  aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTc                                                                                                                     <  2:3343436/150‑1 (MQ=255)
                                   tgggCCTGTGGTGTGGCTTGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACGGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCg                                                                                                                    <  2:398076/148‑1 (MQ=255)
                                      gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGGGCAGCGCCCCCCAAAAAAACCCCAAGGGAGTCGGaa                                                                                                                 >  1:851719/1‑150 (MQ=255)
                                      gCCTGTGGTGGGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                 <  2:2707822/150‑1 (MQ=255)
                                        cTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCTGCGCAGCGCCCCCCCGAAAAGCGCCAAGGCAGTCGGGAGc                                                                                                               >  1:3004498/1‑150 (MQ=255)
                                         tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCAGCGCCCCCCCAAAAAGCGCCAAGGCAGACGGAAGCg                                                                                                              >  2:391467/1‑150 (MQ=255)
                                             gtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCAAAAAAGCGCCCAGGCAGTCCGAAGCGTAGt                                                                                                          >  1:3696301/1‑150 (MQ=255)
                                              tgtggCATGAAAGTGTTAAAGCTCGACCGCGATACCGCCTTGGTGACCGCCGCCGGCAGCGCCATGTGGGGCGCCGCCGCCGTGGTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTg                                                                                                         <  2:4522191/150‑1 (MQ=255)
                                                tggCATTAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                                       <  2:407877/150‑1 (MQ=255)
                                                tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCCAAAACCGCCAAGGCAGTCGGAAGCGTAGTggg                                                                                                       >  1:3982341/1‑148 (MQ=255)
                                                  gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCGAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGt                                                                                                     >  2:309814/1‑150 (MQ=255)
                                                   cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGtt                                                                                                    >  2:646152/1‑150 (MQ=255)
                                                   cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGtt                                                                                                    <  2:1228576/150‑1 (MQ=255)
                                                    aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCGAAAAGCGCCCTGGCAGTCGGAAGCGGAGTGCGGTTc                                                                                                   >  1:3899502/1‑150 (MQ=255)
                                                     tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCg                                                                                                  >  2:3739100/1‑150 (MQ=255)
                                                     tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGACAAGCGCCATGGCAGTCGGAAGCGTAGGGCTGttgg                                                                                                  >  2:3739029/1‑148 (MQ=255)
                                                      gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAGCGCCAGGGCAGTCCGAACCCTAGTGCTGTTCgg                                                                                                 >  2:2446151/1‑150 (MQ=255)
                                                      gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATGGCCGTAGGAAGCGGAGTGCTGTTCgg                                                                                                 >  2:392063/1‑150 (MQ=255)
                                                       aatgtgtACACGCTCGCCCGCGAGCCCGCCTTGCTGACCGCCGCGGGCAGCGCCTGCGGCGGCGCCGCCGCCGTGCTGTCATTCGCAGCGGCGCGGCCCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                                <  1:2501280/147‑1 (MQ=255)
                                                       aaaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                                <  2:2389293/150‑1 (MQ=255)
                                                        aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTa                                                                                               <  2:3219170/150‑1 (MQ=255)
                                                          gtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAcc                                                                                             >  2:2968009/1‑150 (MQ=255)
                                                            gtcctaGCTCGACCGCGATACCGCTTTGCTGACCGCCGCCGGCAGCGCATTCTGGGGGGCCGCCGCCGTGGTGGCATTCGATCCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCtt                                                                                           <  1:1392481/145‑1 (MQ=255)
                                                            gtTCACGCTCGCCCGCGAGCCCGCCTTGCTGCCCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGGGCTGGCATTCGACTCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCtt                                                                                           <  1:3234496/150‑1 (MQ=255)
                                                              taaaGCCCGACCGCGATACGGCTTTGCTGCCCGCCGCCGGCAGCGCATTGTGGGGCGCCGCCGCCGTGCTGTCTTTCTATTGGGGGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGt                                                                                         <  1:3578227/148‑1 (MQ=255)
                                                               cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                        >  2:1585072/1‑150 (MQ=255)
                                                               cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                        >  2:483339/1‑150 (MQ=255)
                                                                 aGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATTGCAGTCGGAATCGTATTGCTGTTCGGTACCTTTTCGa                                                                                      >  2:2251503/1‑150 (MQ=255)
                                                                    tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                   >  2:3418696/1‑150 (MQ=255)
                                                                    tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCGAAAAGCGCCATGGCAGGCGGAAGCGTAGTGGTGTTCGGTACCGTGTCGATGt                                                                                   >  1:2234964/1‑150 (MQ=255)
                                                                    tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCAAAAATCGCCAATGCAGTCGGGAGCGGAGTGCTGTTCGGGACCTTGTAGAGGt                                                                                   >  1:1367575/1‑150 (MQ=255)
                                                                      gcccGCGATCCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCGGCGGCGCCGCCGCGGTGGTGGCATTCGACGCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                 <  1:4133321/148‑1 (MQ=255)
                                                                      gACCGCGATCCCGCCTGGCGGCCCGCCGCGGGCCGCGCCTTGTGGGGCGCCGCCGCGGGGCTGGCTTTCGAGGCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                 <  1:179202/150‑1 (MQ=255)
                                                                        ccGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAACTTGTCGATGTTCCt                                                                               >  2:809615/1‑150 (MQ=255)
                                                                         cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGTTCCtt                                                                              >  2:1821264/1‑150 (MQ=255)
                                                                             atcccGCGTGGCTGCCCGCCGCCGGCCGCGCCTGCTGCGGCGCCGCCGCGGTGGTGTAATGCTAATGGGCGGTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAcc                                                                          <  1:3420854/147‑1 (MQ=255)
                                                                             aTACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAcc                                                                          <  2:505991/150‑1 (MQ=255)
                                                                               aCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCa                                                                        >  2:2294759/1‑150 (MQ=255)
                                                                                  gctttGCTGACCGCCGCCGGAAGCGCATTCTGCGGCGCCGCGGCGGTGGTGTCATCCGAAGCGGCGCTCCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                     <  1:846006/147‑1 (MQ=255)
                                                                                  gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCAAAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCTTGTAGATGTTCCTTTAACCAcag                                                                     >  2:1248331/1‑148 (MQ=255)
                                                                                  gCCTGGCTGACCGCCGCCGGCAGCGCCAGCTGCGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                     <  1:4514952/150‑1 (MQ=255)
                                                                                   ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTgg                                                                    <  2:3632773/150‑1 (MQ=255)
                                                                                   ccTTGCTGACCGCCGCCGGCAGCCACATCTGCGGCGCCGCCGCCGTGCTGGAATTCTAATCGGCGCTGCGCAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTgg                                                                    <  2:390906/150‑1 (MQ=255)
                                                                                     ttGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCg                                                                  <  2:397858/150‑1 (MQ=255)
                                                                                        cTGCCCGCCGCGGGCAGCGCCCTCGGCGGCGCCGCCGCGGTGCTGGCATGCGCAGGGGCGCTGCGCACCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATc                                                               <  1:4276155/150‑1 (MQ=255)
                                                                                           aCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAAAAAAGCCCCAAGGCAGTCGGACGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCGTTACACACAGGCGATaaaa                                                            >  1:4443776/1‑149 (MQ=255)
                                                                                              gccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAACGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGcc                                                         <  2:3357694/150‑1 (MQ=255)
                                                                                                cgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAACCACTGGCGATCAAttaccg                                                       >  2:2265766/1‑145 (MQ=255)
                                                                                                  ccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                                     <  2:1444858/150‑1 (MQ=255)
                                                                                                         cgcCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTTCCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCGGCa                                              >  2:4019864/1‑150 (MQ=255)
                                                                                                            cATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGATGGTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                           >  2:292606/1‑150 (MQ=255)
                                                                                                               cggcggcgCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                        <  2:2821528/148‑1 (MQ=255)
                                                                                                               cTGCGGCGCCGCCGCCGTGCTGTCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                        <  2:851719/150‑1 (MQ=255)
                                                                                                                     cgccgccgcCGTGCGGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAt                                  <  2:3185669/150‑1 (MQ=255)
                                                                                                                      gccgccgccGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                                 <  2:392871/150‑1 (MQ=255)
                                                                                                                        cgccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATAAATTCCGGCTGGCTGCATCTGGATACCATggg                               >  2:732056/1‑150 (MQ=255)
                                                                                                                         gccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAATCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTAGATTTTCCATTACCCACAGGCGATCAAATCCGGCTGGCTGCATCTGGATAACAAGGgc                              >  2:1399273/1‑150 (MQ=255)
                                                                                                                         gccgccgGGGGGGAATTCTAAGCGGCGCTGCGCCGCCCCCCCCATAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGgc                              <  1:603196/150‑1 (MQ=255)
                                                                                                                              ggTGCTGGCATTCGAGCCGGCGGGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                         <  1:4177330/149‑1 (MQ=255)
                                                                                                                               gTGCTGGCTTTCTATTCGGCGCTGCGCCGCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                        <  1:3739029/150‑1 (MQ=255)
                                                                                                                               gTGCTGGCATTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                        <  2:4187125/150‑1 (MQ=255)
                                                                                                                               gTGCTGGATTTCGAAGCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                        <  1:3739100/150‑1 (MQ=255)
                                                                                                                                tGCTGGCATTAGATGCGGCGCTGGGATGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggtc                       <  2:3057633/150‑3 (MQ=255)
                                                                                                                                 gCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgct                      <  2:4510366/150‑1 (MQ=255)
                                                                                                                                   tGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCAGGATACCATGGGCGCCGGgctgc                    >  2:1080440/1‑150 (MQ=255)
                                                                                                                                   tGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCGTTACCCACAGGCGATCAATGCCGGCTGGCAGCATCTGGATACCAAGGGCGCCGGgctgc                    >  2:2492826/1‑150 (MQ=255)
                                                                                                                                   tGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGGAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATAACATGGGCGCCGGgctgc                    >  2:2903680/1‑150 (MQ=255)
                                                                                                                                   tGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAATGCAGTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGAAGTTCCTTTACCCACAGGCGATAAAATCCCGCTGGCTGCATATGGATAACATGGGCGCCGggcggc                    >  2:2903684/1‑150 (MQ=255)
                                                                                                                                      ctttCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgc                 <  2:3410411/148‑1 (MQ=255)
                                                                                                                                      cATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgc                 <  2:2606621/150‑1 (MQ=255)
                                                                                                                                       aTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCATTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgct                >  2:2134634/1‑150 (MQ=255)
                                                                                                                                        gccgccgcggcgcggcgCCGCGCCCCCCCTAACAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCGGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg               <  1:1415477/143‑1 (MQ=255)
                                                                                                                                           gattCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc            <  1:1106471/147‑1 (MQ=255)
                                                                                                                                            aaTCGGCGCTGCGCAGCGCCCCCCATAAAACCGCCATGGCAGTCGGAAGCCGAGTGCTGTTCGGTACCATGTCGATGTTCCTTTAAACAAAGGCGATCAATGCCGGCTGGCTGCATCTCGATACAAACGGCGCCGGGCGGCTGCTGggcg           >  1:2778254/1‑150 (MQ=255)
                                                                                                                                              tCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcggc         >  2:22186/1‑150 (MQ=255)
                                                                                                                                                 gcgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt      >  1:1974655/1‑150 (MQ=255)
                                                                                                                                                   gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt    >  2:3747752/1‑150 (MQ=255)
                                                                                                                                                    cTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTc   <  2:4101069/150‑1 (MQ=255)
                                                                                                                                                     tGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc  >  2:640558/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCC  >  NC_002947/4160283‑4160581

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