Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,432 T→G 22.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,4320TG22.9% 311.2 / 21.1 186L159R (CTG→CGG) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base T (57/81);  new base G (31/16);  total (89/97)
Fisher's exact test for biased strand distribution p-value = 4.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

GTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCC  >  NC_002947/4160286‑4160581
                                                                                                                                                  |                                                                                                                                                      
gtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGTCCGCCGCCGGCATCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTGGTATCGGCTCTgcg                                                                                                                                                     <  2:4959025/150‑1 (MQ=255)
  cGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCa                                                                                                                                                   <  2:209441/150‑1 (MQ=255)
    gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                                 <  2:6727744/150‑1 (MQ=255)
    gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                                 <  2:8951268/150‑1 (MQ=255)
           tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGCATCGGCGCGGCGCCGCG‑cccccc                                                                                                                                         >  2:9042385/1‑150 (MQ=255)
            aGTCCCCAGTACCCTGCTGTGCGGCCTGTGGTGTGGCATGAAAGTTTTCATGCTCGCCCGCTATACCGCGTGGCTGACGGCCGCGGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCTATTCGGCGCTgcgccgcg‑ccgattc                                                                                                                                        <  2:1270538/150‑6 (MQ=255)
             gTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCAt                                                                                                                                       >  2:4796574/1‑150 (MQ=255)
              tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑cccccccac                                                                                                                                      >  1:7787564/1‑147 (MQ=255)
                aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATaaa                                                                                                                                    <  2:2775216/150‑1 (MQ=255)
                 ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATaaaa                                                                                                                                   >  2:25836/1‑150 (MQ=255)
                 ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATaaaa                                                                                                                                   >  2:26091/1‑150 (MQ=255)
                  cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGTGGGCAGCG‑CGCCCCATAAAAg                                                                                                                                  >  1:3714161/1‑150 (MQ=255)
                  cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCCAAAAg                                                                                                                                  >  1:4458515/1‑150 (MQ=255)
                   aGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAgc                                                                                                                                 <  2:5836456/150‑1 (MQ=255)
                       cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCa                                                                                                                             >  2:3934612/1‑150 (MQ=255)
                        ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAAGCGCCAt                                                                                                                            >  2:2484626/1‑150 (MQ=255)
                        ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCaa                                                                                                                            >  2:1788605/1‑149 (MQ=255)
                        ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCCAAAAACCGccca                                                                                                                            >  1:4650955/1‑148 (MQ=255)
                        ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAAAAAAGCGCCaa                                                                                                                            >  2:7948709/1‑149 (MQ=255)
                          tgctgTTGGGCTTGTGGTGTGGCATGAAAGTGTTCATGCTCGACGGCGATACCGCGTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTGGAATCGGCGCTGCGCTGCG‑CCCCCCATAAAAGCGCCATgg                                                                                                                          <  2:8104995/150‑1 (MQ=255)
                          tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCAGAAAAGCGCCAAgg                                                                                                                          >  1:624982/1‑150 (MQ=255)
                           gctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCAAAAAGCGcccaggc                                                                                                                         >  1:2298556/1‑145 (MQ=255)
                            gtgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCGGTGGTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCa                                                                                                                        <  2:6899163/149‑1 (MQ=255)
                                   gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGAAGCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                 <  2:7082745/150‑1 (MQ=255)
                                   gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                 <  2:7082765/150‑1 (MQ=255)
                                   gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                 >  2:3688902/1‑150 (MQ=255)
                                   gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                 >  2:8987654/1‑150 (MQ=255)
                                    ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCAAAAAAAACCCAAGGCAGTCGGAAg                                                                                                                >  1:4060853/1‑150 (MQ=255)
                                    ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCACTCGGAAg                                                                                                                >  2:67571/1‑150 (MQ=255)
                                    ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGAGCCAAGGCAGTCGGAAg                                                                                                                >  1:7617022/1‑150 (MQ=255)
                                      tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCAAAAAGCGCCAAGGCAGGCGGAAGCg                                                                                                              >  1:7303541/1‑150 (MQ=255)
                                      tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCC‑CCCCCCCAAAAACCCCCCAGGCACGCGGaaccc                                                                                                              >  1:6239333/1‑147 (MQ=255)
                                      tgtggtgtggATTGAATGTTTTCAAGCCCGACGGCGATACGGCGTTGCTGACCCCGGCGGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCTTTCTATGCGGGGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCg                                                                                                              <  1:2796172/150‑1 (MQ=255)
                                      tgtggtggggCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATGTGCGGCGCCGCCGCCGTGCTGGCTTTCGACTGGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCg                                                                                                              <  1:2795606/150‑1 (MQ=255)
                                        tggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCACCG‑CCCCCCCGAAAAGAGCCCAGGCACGCGGAAGCCTa                                                                                                            >  1:7372600/1‑150 (MQ=255)
                                          gtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGt                                                                                                          <  2:6316682/150‑1 (MQ=255)
                                           tgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CTCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTg                                                                                                         <  1:4085946/150‑1 (MQ=255)
                                           tgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CTCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTg                                                                                                         <  1:4086083/150‑1 (MQ=255)
                                             tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAGAAAACCGCCAAGGCCGTCGGAAGCGTAGTGcg                                                                                                       >  2:1881515/1‑149 (MQ=255)
                                             gggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATCAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                                       <  2:3171868/149‑1 (MQ=255)
                                                cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCGCCCCCCCAAAAAACCCCCAGGGC‑GTCGGAAGCGGAGTGCTGtt                                                                                                    >  1:3121188/1‑150 (MQ=255)
                                                 aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTTGTGCTGTTc                                                                                                   >  2:2836051/1‑150 (MQ=255)
                                                 aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                                   >  2:7846012/1‑150 (MQ=255)
                                                  ttattgtgtTCAAGCTCGACCGCGATACCGCCGTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAACCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCg                                                                                                  <  1:1495443/145‑1 (MQ=255)
                                                  tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCCAAAAAGCGCCCTGGCCGTCGGAAGCGTAGTGGTGTTCg                                                                                                  >  2:3343990/1‑150 (MQ=255)
                                                   gataGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCgg                                                                                                 <  1:6864294/147‑1 (MQ=255)
                                                    taaGTTTTCTTGCTGGACCGCGATACCGCGTTGGTGACCGCCGCCGGCAGCGCCATGTGCGGCGCCGCCGCCGTTGTGGAATTCGAATGGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                                <  1:8311536/149‑1 (MQ=255)
                                                    aaaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGGCGGAAGCGTAGTGCTGTTCGGt                                                                                                >  2:1631200/1‑150 (MQ=255)
                                                      aGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAc                                                                                              <  2:5640638/150‑1 (MQ=255)
                                                       gtgtTCACGCTCGCCCGCGCGCCCGCCTTGCTGCCCGCCGCCGGCCGCGCCATGGGCGGCGCCGCCGCGGTGCGGGCTTTCGACGCGGGGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAcc                                                                                             <  1:7969214/150‑1 (MQ=255)
                                                          ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTg                                                                                          >  2:9056703/1‑150 (MQ=255)
                                                          ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAGAAAAGCGCCAGGGCAGTCGGGAGCGTAGTGCTGTTCGGTAccaag                                                                                          >  1:8201450/1‑147 (MQ=255)
                                                            cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                        >  2:3811238/1‑150 (MQ=255)
                                                             aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCTTGTCg                                                                                       >  2:1398508/1‑150 (MQ=255)
                                                               gCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                     >  2:7923211/1‑150 (MQ=255)
                                                                 tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGt                                                                                   >  2:2509218/1‑150 (MQ=255)
                                                                 ccgcccGCGATCCCGCGTTGCGGCCCGCGGCGGGCTGCGCCTTGGGGGGCGCCGCCGCGGGGGGGGATTTCGACGCGGGGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                   <  1:4620497/146‑1 (MQ=255)
                                                                    cccGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTcc                                                                                <  1:5964590/149‑1 (MQ=255)
                                                                        cgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTa                                                                            >  2:5469830/1‑150 (MQ=255)
                                                                         gATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCG‑CCCCCCCAAAAATCGCCAAGGCCGGCGGGAGAGTAGTTGTGTTGGGGACAATGTCGATGTTGCGTTaa                                                                           >  1:4846818/1‑149 (MQ=255)
                                                                             ccGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGGAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCAc                                                                       <  2:5762557/150‑1 (MQ=255)
                                                                             ccGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCAc                                                                       <  2:2818361/150‑1 (MQ=255)
                                                                              cGCCTTGGTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACt                                                                      <  2:8529861/150‑1 (MQ=255)
                                                                              cGCCTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCCGCCGCCGTGCGGGCATTCGAAGCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACt                                                                      <  2:3271148/150‑1 (MQ=255)
                                                                              cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACt                                                                      >  2:9167805/1‑150 (MQ=255)
                                                                                   ggCTGACCGCCGCCGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGATTTCTAAGCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGa                                                                 <  1:8961372/149‑1 (MQ=255)
                                                                                      tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTGGGGACCTTGTCGATTTTGCGTTTCCCACAGGCGCTCa                                                              >  1:7581979/1‑150 (MQ=255)
                                                                                          cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGc                                                          <  2:609993/150‑1 (MQ=255)
                                                                                           gccgcggGCAGCGCCATCTGCGGCGCCGCCGCCGTGGTGGCATTCTAATCGGCGGTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGcc                                                         <  2:4211767/150‑1 (MQ=255)
                                                                                              gccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggc                                                      <  2:117003/150‑1 (MQ=255)
                                                                                                cgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCC‑CCCCCCAGAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAACCCCAGGCGATCAATGCCggctg                                                    >  2:1309121/1‑150 (MQ=255)
                                                                                                 ggCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctgg                                                   <  2:7738364/150‑1 (MQ=255)
                                                                                                   cAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggct                                                 <  2:1835220/150‑1 (MQ=255)
                                                                                                         cATCGGGGGCGCCGCCGCCGTGCTGGCATTCGTCGCGGGGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                           <  1:1309121/150‑1 (MQ=255)
                                                                                                          agggggggcgCCGCCGGGGGGCTGGCTTGCGACGGGGCGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg                                          <  1:3345078/144‑1 (MQ=255)
                                                                                                           tCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCGGGTCGGGACCTTGTCACTGTTCCGTTACCCACGGGCGCTCAATGCCGGCGGGCGGCATCTgg                                         >  2:2257622/1‑150 (MQ=255)
                                                                                                            cggcggcgCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                        <  2:7193036/148‑1 (MQ=255)
                                                                                                             tGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGAAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACACAATGGCGATCAATGCCGGCTGGCTGCATCTGGAt                                       >  2:7032580/1‑150 (MQ=255)
                                                                                                                 gcgcCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCa                                   <  2:8625317/150‑1 (MQ=255)
                                                                                                                   gccgccgccGTTCTGGAATTCGAATGGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                                 <  1:1896685/150‑1 (MQ=255)
                                                                                                                   gccgccgccGTGCTGGCATTCGATTCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                                 <  1:8600246/150‑1 (MQ=255)
                                                                                                                            gTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                        <  2:591252/150‑1 (MQ=255)
                                                                                                                             ggCTGGCATCCGACTCGGCGCTGCGCAGCG‑CCCCCCATACCAGCGCCATGGCCGTCGGAAGCGTCGTGCGGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgc                       <  2:3668191/149‑1 (MQ=255)
                                                                                                                                    aTTCGAATCGGCGCTGCGCCGCG‑CCCCCCATAAAATCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgct                >  1:8391613/1‑150 (MQ=255)
                                                                                                                                     ttCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAATCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAACCACAGGCGATAAAATCCGGCTGGCTGCATCTGGATAACAAGGGCGCCGGGctgctgctg               >  2:4227513/1‑150 (MQ=255)
                                                                                                                                      tatggggggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTACTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg              >  1:6121921/8‑150 (MQ=255)
                                                                                                                                      tCGATTGGGGGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg              <  1:1337946/150‑1 (MQ=255)
                                                                                                                                         aaTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg           <  2:8321236/150‑1 (MQ=255)
                                                                                                                                          gggggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAAGGGCGCCGGGCTGCTGCTGggcgg          >  1:1001851/4‑150 (MQ=255)
                                                                                                                                          cgcggcgcggcgCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg          <  1:4978001/148‑1 (MQ=255)
                                                                                                                                          aTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg          >  2:3118490/1‑150 (MQ=255)
                                                                                                                                          aTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg          <  2:624952/150‑1 (MQ=255)
                                                                                                                                            gggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa        >  1:435756/2‑150 (MQ=255)
                                                                                                                                            gggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa        >  1:435511/2‑150 (MQ=255)
                                                                                                                                             ggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAACCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCAc       >  1:3138459/1‑150 (MQ=255)
                                                                                                                                                gcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt    >  2:7273526/1‑150 (MQ=255)
                                                                                                                                                gcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt    >  2:4437646/1‑150 (MQ=255)
                                                                                                                                                gcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt    >  1:6656659/1‑150 (MQ=255)
                                                                                                                                                 cTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTc   <  1:4922752/150‑1 (MQ=255)
                                                                                                                                                  tGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc  >  1:3752437/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCC  >  NC_002947/4160286‑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: ‑