Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,006,511 T→G 9.1% F1017V (TTC→GTC)  PP_2627 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,006,5110TG 9.1% 235.2 / 14.2 88F1017V (TTC→GTC) PP_2627membrane protein
Reads supporting (aligned to +/- strand):  ref base T (46/34);  new base G (5/3);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

GCTGACGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATC  >  NC_002947/3006364‑3006660
                                                                                                                                                   |                                                                                                                                                     
gCTGACGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTGCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTc                                                                                                                                                     >  1:962256/1‑150 (MQ=255)
         tACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCGGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGAc                                                                                                                                            >  1:1755435/1‑150 (MQ=255)
         tACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGAc                                                                                                                                            >  2:1876363/1‑150 (MQ=255)
         tACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGAc                                                                                                                                            >  2:201204/1‑150 (MQ=255)
          aCAGTTTTTTCCAGCAACGCCTGGCAGACCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACa                                                                                                                                           >  2:1166457/1‑150 (MQ=255)
            aGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAAc                                                                                                                                         <  2:456051/150‑1 (MQ=255)
                 tttCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGa                                                                                                                                    <  1:1686094/150‑1 (MQ=255)
                 tttCCAGCAACGCCTGGCAGACCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGa                                                                                                                                    <  1:1166457/150‑1 (MQ=255)
                  ttCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAg                                                                                                                                   >  1:983819/1‑150 (MQ=255)
                    ccAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGc                                                                                                                                 <  2:1689315/150‑1 (MQ=255)
                    ccAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGc                                                                                                                                 >  1:167904/1‑150 (MQ=255)
                     cAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGcc                                                                                                                                >  1:644899/1‑150 (MQ=255)
                       gCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCtt                                                                                                                              >  2:228372/1‑150 (MQ=255)
                       gCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCtt                                                                                                                              >  1:779903/1‑150 (MQ=255)
                        cAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTg                                                                                                                             <  1:632345/150‑1 (MQ=255)
                         aaCGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTGTCCTTCAGGAAAACCTTGAgataccga                                                                                                                            <  2:711981/150‑8 (MQ=39)
                          aCGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCa                                                                                                                           >  1:1652848/1‑150 (MQ=255)
                            gCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGt                                                                                                                         <  1:1441926/150‑1 (MQ=255)
                            gCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGt                                                                                                                         >  1:724905/1‑150 (MQ=255)
                               tGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGCACAGGCGCCTGCAGCAGTTTAATGCTCGGTATACGAATGTATGCCGGAAGGACAACCTTGAGGCCTTGCAGTCCg                                                                                                                      >  1:1867494/1‑150 (MQ=255)
                                   aGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCt                                                                                                                                                   >  2:372356/1‑117 (MQ=255)
                                   aGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCt                                                                                                                                                   <  1:372356/117‑1 (MQ=255)
                                         cTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACgg                                                                                                            <  2:1391538/150‑1 (MQ=255)
                                          tACCCCTTTGTGGTGCAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTATTCAGGATTATTTGAGTAGTTGATGAGTTTAATGATCTGTATATCAAGGTTTTGCCGAAGGCCAACCTGGAGTCCTTGCAGTCCGCTCAGCACggg                                                                                                           <  1:1876363/150‑1 (MQ=255)
                                          tACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGTGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACggg                                                                                                           <  1:376257/150‑1 (MQ=255)
                                             cccTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAACGGGCGCCTGCAGCAGTTTCATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAACACGGGCaa                                                                                                        >  1:1342134/1‑150 (MQ=255)
                                             cccTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCaa                                                                                                        >  1:499500/1‑150 (MQ=255)
                                                tttGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATcc                                                                                                     <  2:779903/150‑1 (MQ=255)
                                                    tggtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGa                                                                                                 <  1:123478/150‑1 (MQ=255)
                                                    tggtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGa                                                                                                 <  1:201204/150‑1 (MQ=255)
                                                      gtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGAtt                                                                                               >  2:1922762/1‑150 (MQ=255)
                                                           ggCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGc                                                                                          <  1:350936/150‑1 (MQ=255)
                                                              aCCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTACCGa                                                                                       >  2:279275/1‑150 (MQ=255)
                                                              aCCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGATATCTCAAGGATTTCCTGAAGGAAAACCATGAGGCCTGGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGa                                                                                       >  1:890261/1‑150 (MQ=255)
                                                                ccGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCACGCAATTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCATTCCCGGATTCGTGCCGACg                                                                                     >  1:1861127/1‑150 (MQ=255)
                                                                  gATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGt                                                                                                                         >  1:886161/1‑112 (MQ=255)
                                                                  gATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGt                                                                                                                         <  2:886161/112‑1 (MQ=255)
                                                                  gATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTTATTCGTGCCGACGTg                                                                                   >  1:391675/1‑150 (MQ=255)
                                                                     gCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATc                                                                                <  2:499500/150‑1 (MQ=255)
                                                                      cttcttTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCg                                                                               >  1:1604610/1‑150 (MQ=255)
                                                                       ttcttTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCtt                                                                                                                                      <  1:600013/94‑1 (MQ=255)
                                                                       ttcttTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCtt                                                                                                                                      >  2:600013/1‑94 (MQ=255)
                                                                               aTGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCg                                                                      >  1:1604480/1‑150 (MQ=255)
                                                                                tGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGa                                                                     <  2:1652848/150‑1 (MQ=255)
                                                                                    ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCa                                                                                                               <  1:1729202/104‑1 (MQ=255)
                                                                                    ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCa                                                                                                               >  2:1729202/1‑104 (MQ=255)
                                                                                    ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTATCTCGTATGCCGTCTTAAGGTTTGCCTGAATGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTg                                                                 >  1:976558/1‑150 (MQ=255)
                                                                                    ttCGAGGCGTTCTTCGGGCAGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTg                                                                 >  1:1193126/1‑150 (MQ=255)
                                                                                          gCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAg                                                           >  2:1070324/1‑150 (MQ=255)
                                                                                           cggtcgtcGGGCCGAAAGGGCGCCTGCAGGCGTTTGGTGGTCGGTATGCCGCGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgc                                                          <  2:1755435/144‑1 (MQ=255)
                                                                                            gTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTGTCCTTCAGGAAAACCTTGa                                                                                                                                    >  2:413218/1‑75 (MQ=25)
                                                                                            gTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTGTCCTTCAGGAAAACCTTGa                                                                                                                                    <  1:413218/75‑1 (MQ=37)
                                                                                             ttcttcGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgcgc                                                        >  2:1529777/1‑150 (MQ=255)
                                                                                                    ggCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCcgcg                                                 >  2:1456794/1‑150 (MQ=255)
                                                                                                     gCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAg                                                                                                                 >  1:130290/1‑85 (MQ=255)
                                                                                                     gCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAg                                                                                                                 <  2:130290/85‑1 (MQ=255)
                                                                                                      ccGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGaa                                               >  2:1233277/1‑150 (MQ=255)
                                                                                                       cGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGaaa                                              <  1:1386388/150‑1 (MQ=255)
                                                                                                                gcCTGCAGCAGTTTAATGATCGGCATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCg                                     >  2:951352/1‑150 (MQ=255)
                                                                                                                 ccTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGa                                    <  1:1233277/150‑1 (MQ=255)
                                                                                                                  cTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGAt                                   >  1:1244450/1‑150 (MQ=255)
                                                                                                                         agTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTg                            <  2:1867494/150‑1 (MQ=255)
                                                                                                                         agTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTg                            <  2:1668781/150‑1 (MQ=255)
                                                                                                                         agTTTAATGATCGGTATCTCAAGGCTTTCCTGAAGGACACCCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGGATGCGCGAAACCGTCTTCGATCAGCGTg                            <  2:962256/150‑1 (MQ=255)
                                                                                                                           tttAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGc                                                                   >  2:522300/1‑109 (MQ=255)
                                                                                                                           tttAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGc                                                                   <  1:522300/109‑1 (MQ=255)
                                                                                                                           tttAATGATCGGTATATCAAGGTTTTCCTGAAGGACAACCTTGAGGCCtt                                                                                                                              <  2:1923764/50‑1 (MQ=255)
                                                                                                                           tttAATGATCGGTATATCAAGGTTTTCCTGAAGGACAACCTTGAGGCCtt                                                                                                                              >  1:1923764/1‑50 (MQ=255)
                                                                                                                             tAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCaa                        >  2:25517/1‑150 (MQ=255)
                                                                                                                                  aTCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGNTaa                   <  2:1193126/150‑1 (MQ=255)
                                                                                                                                   tCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTg                            >  2:1920858/1‑140 (MQ=255)
                                                                                                                                   tCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTg                            <  1:1920858/140‑1 (MQ=255)
                                                                                                                                   tCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAg                  >  2:1616509/1‑150 (MQ=255)
                                                                                                                                    cGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAgt                 <  1:626433/150‑1 (MQ=255)
                                                                                                                                     ggTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAgtg                >  1:1545574/1‑150 (MQ=255)
                                                                                                                                     aaggcctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAgtg                >  1:285545/6‑150 (MQ=255)
                                                                                                                                      gTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAgtgt               >  1:1804070/1‑150 (MQ=255)
                                                                                                                                        aTCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtgt             <  2:1604480/150‑2 (MQ=255)
                                                                                                                                          ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTgg                                                                >  2:392605/1‑97 (MQ=38)
                                                                                                                                          ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTgg                                                                <  1:392605/97‑1 (MQ=38)
                                                                                                                                          ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtggag           >  1:1258575/1‑147 (MQ=255)
                                                                                                                                          ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg           >  1:1328121/1‑150 (MQ=255)
                                                                                                                                           tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAgt                 >  2:1134762/1‑143 (MQ=255)
                                                                                                                                           tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAgt                 <  1:1134762/143‑1 (MQ=255)
                                                                                                                                           tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt          >  1:1007450/1‑150 (MQ=255)
                                                                                                                                             aaGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTc        >  1:298254/1‑150 (MQ=255)
                                                                                                                                             aaGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTc        >  1:10153/1‑150 (MQ=255)
                                                                                                                                              aGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCa       <  2:1007450/150‑1 (MQ=255)
                                                                                                                                                  tttCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCAt   <  1:1695236/150‑1 (MQ=255)
                                                                                                                                                   ttCCTGAAGGACAACCTTGAGGCTTTGGATTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATc  >  1:991019/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTGACGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATC  >  NC_002947/3006364‑3006660

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