Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,038,537 A→T 6.5% D29V (GAC→GTC)  PP_2646 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,038,5370AT 6.5% 1150.8 / 45.5 359D29V (GAC→GTC) PP_2646hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (170/162);  new base T (10/13);  total (184/175)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 381 total aligned reads displayed.

CTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGCT  >  NC_002947/3038388‑3038686
                                                                                                                                                     |                                                                                                                                                      
      ggTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTcca                                                                                                                                                  >  1:8077236/1‑150 (MQ=255)
          gCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgc                                                                                                                                              <  1:4893206/150‑1 (MQ=255)
          gCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgc                                                                                                                                              >  1:2579788/1‑150 (MQ=255)
            ccAGGCCATGAAGACCCTCTATCTCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgcca                                                                                                                                           <  2:5764672/150‑1 (MQ=255)
              aGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCACGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgccagc                                                                                                                                         >  2:4328390/1‑150 (MQ=255)
              aGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCTGGCGCTGGACGTCGTccagggcctggc                                                                                                                                         <  2:505293/150‑5 (MQ=39)
                 ccATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCCGCGGc                                                                                                                                      >  1:6342151/1‑150 (MQ=255)
                           ccTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCCGGCCCTGGACGACGTCCAGCGCCAGCGGCCCATCCGGgt                                                                                                                            >  1:3029624/1‑150 (MQ=255)
                           ccTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGgt                                                                                                                            <  2:3465896/150‑1 (MQ=255)
                           ccTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGgt                                                                                                                            >  2:1873809/1‑150 (MQ=255)
                           ccTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGgt                                                                                                                            >  2:1456351/1‑150 (MQ=255)
                           ccTCTAGATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGgt                                                                                                                            <  1:4750829/150‑1 (MQ=255)
                             tctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGgtgt                                                                                                                          <  2:4141535/150‑1 (MQ=255)
                             tctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGgtgt                                                                                                                          <  1:2339325/150‑1 (MQ=255)
                               tataTCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTg                                                                                                                        >  2:7485798/1‑150 (MQ=255)
                                   tCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATcc                                                                                                                    >  2:4529644/1‑150 (MQ=255)
                                   tCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATcc                                                                                                                    >  2:5361042/1‑150 (MQ=255)
                                      gccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACg                                                                                                                 >  2:4564388/1‑150 (MQ=255)
                                      gccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACg                                                                                                                 >  1:4939316/1‑150 (MQ=255)
                                      gccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACg                                                                                                                 >  2:8080970/1‑150 (MQ=255)
                                       ccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACgg                                                                                                                >  2:363277/1‑150 (MQ=255)
                                       ccgTTTCGGCCAGGGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGGTGCGGGTGTTGATCCACgg                                                                                                                >  2:2976204/1‑150 (MQ=255)
                                          tttCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTgg                                                                                                             >  2:7336868/1‑150 (MQ=255)
                                                  aGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGt                                                                                                     >  1:5411193/1‑150 (MQ=255)
                                                    tGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGttt                                                                                                   >  2:6555994/1‑150 (MQ=255)
                                                       ttCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCTGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCtt                                                                                                <  2:6486643/150‑1 (MQ=255)
                                                       ttCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCtt                                                                                                <  2:7843536/150‑1 (MQ=255)
                                                       ttCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCtt                                                                                                <  1:1367749/150‑1 (MQ=255)
                                                       ttCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCtt                                                                                                <  1:5787279/150‑1 (MQ=255)
                                                       gtCTGAGCGCGGCAAGGAGAGCAGACCATGCGAGTACTGAGATGTGCAGGGCGCAACGACGCCGACCCCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCtt                                                                                                <  2:5669702/149‑1 (MQ=255)
                                                        tCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGTGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                               <  2:6364875/150‑1 (MQ=255)
                                                          tGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTTATCCACCGTGGCAGCCAGTTTCTtggg                                                                                             >  1:1071196/1‑149 (MQ=255)
                                                          tGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtggt                                                                                             <  1:5372869/150‑1 (MQ=255)
                                                            aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGGTTCTtggtgg                                                                                           >  1:517832/1‑150 (MQ=255)
                                                              cGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTgggg                                                                                         >  2:1504028/1‑150 (MQ=255)
                                                               gATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGa                                                                                        <  1:8499371/150‑1 (MQ=255)
                                                                aTGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGaa                                                                                       <  1:3277110/150‑1 (MQ=255)
                                                                  gCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAAt                                                                                     >  2:8220022/1‑150 (MQ=255)
                                                                   cAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                                    <  1:3068286/150‑1 (MQ=255)
                                                                   cAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                                    <  1:2740662/150‑1 (MQ=255)
                                                                   cAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                                    <  1:179022/150‑1 (MQ=255)
                                                                   cAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                                    <  2:8633887/150‑1 (MQ=255)
                                                                   cAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                                    <  1:8812529/150‑1 (MQ=255)
                                                                      ggAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGaa                                                                                 <  1:13399/150‑1 (MQ=255)
                                                                            cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTgg                                                                           >  2:5879801/1‑150 (MQ=255)
                                                                            cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTgg                                                                           >  2:7136598/1‑150 (MQ=255)
                                                                            cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTgg                                                                           <  2:3208853/150‑1 (MQ=255)
                                                                            cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTgg                                                                           <  1:7225528/150‑1 (MQ=255)
                                                                             aGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTggg                                                                          >  2:1765823/1‑150 (MQ=255)
                                                                                ccATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCAGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACCGTGGCAGCCAGTTTCTTGGTGGGGAAAACGAAAGCTGGgcgc                                                                       >  1:5546414/1‑150 (MQ=255)
                                                                                  aTGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGgcgcgc                                                                     >  1:5041816/1‑150 (MQ=255)
                                                                                    gCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGa                                                                   <  1:6691991/150‑1 (MQ=255)
                                                                                       aGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCa                                                                <  1:4248702/150‑1 (MQ=255)
                                                                                       aGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCa                                                                >  2:2414338/1‑150 (MQ=255)
                                                                                            tGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCAGGACGACGTCCAGCGCCACCGGCCGATCCGGGTGTTGATACACGGTGGCAGCCAGTTTATTGGTGGGGAAATCGAAAGCTTGGCGCGCGAGCATGgcg                                                           >  1:8339951/1‑150 (MQ=255)
                                                                                            tGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCACGCCCAGGACGACCTCCAGCGCCAGCGGCCCCACCGGGTGTTGATCCACCGTGGCAGCCAGTTTCGTGGTGGGGAAATCGACAGCTGGGCGCGCGCGCAGGgcg                                                           >  1:6618031/1‑150 (MQ=255)
                                                                                             gATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGgcgc                                                          <  2:5781915/150‑1 (MQ=255)
                                                                                             gATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCCCCACGCCCGGGACGACGTCCAGCGCCCGCCGCCGATCCGGGGGTTGATCCACGGGGGCAGCCCCTTTGTTGGTGGGGAACTCGAAAGCTGGGCGCGCGAGCAGGgcgc                                                          >  1:4080852/1‑150 (MQ=255)
                                                                                              atatGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGcc                                                         >  2:3410847/1‑150 (MQ=255)
                                                                                                 tgtgCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGAt                                                      <  2:7974213/150‑1 (MQ=255)
                                                                                                     cAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGGCCACCGCCCGCGGCCCATCCCGGTGTTGATCCACGGTGGCCGCCACTTTCTTTGTGGGGAAAACGAAAGCTGGGCGCGCGAGCCTGGCGCCGATCAgg                                                  >  1:4463624/1‑150 (MQ=255)
                                                                                                       gggCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGCCGACGTCCAGCGCCCGCGGCCGATCCGGGTGTTGTTACACGGTGGCAGCCAGTTTATTGGTGGGGGAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTg                                                >  1:6275498/1‑150 (MQ=255)
                                                                                                        ggCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGc                                               <  1:206794/150‑1 (MQ=255)
                                                                                                        ggCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGc                                               <  2:6182443/150‑1 (MQ=255)
                                                                                                          cgcAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGt                                             <  2:6816654/150‑1 (MQ=255)
                                                                                                             aaCTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTAc                                          >  2:3371528/1‑150 (MQ=255)
                                                                                                             aaCTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATACACGGTGGCAGCCACTTTCTTGGTGGGGAAATAGAAATCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTAc                                          >  1:982554/1‑150 (MQ=255)
                                                                                                               cTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTAccc                                        <  1:7883362/150‑1 (MQ=255)
                                                                                                                tACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCg                                       >  1:8608572/1‑150 (MQ=255)
                                                                                                                   gCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCa                                    <  1:7551807/150‑1 (MQ=255)
                                                                                                                       acacCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACt                                <  2:4494072/150‑1 (MQ=255)
                                                                                                                        cacCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTg                               <  2:5559020/150‑1 (MQ=255)
                                                                                                                        cacCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTg                               <  2:4851005/150‑1 (MQ=255)
                                                                                                                         acCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTgg                              >  1:8046940/1‑150 (MQ=255)
                                                                                                                         acCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTgg                              >  2:114496/1‑150 (MQ=255)
                                                                                                                          cccGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGc                             <  1:2227560/150‑1 (MQ=255)
                                                                                                                           ccGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCa                            >  2:2131224/1‑150 (MQ=255)
                                                                                                                           ccGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCa                            <  1:8273264/150‑1 (MQ=255)
                                                                                                                           ccGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCa                            <  2:8509679/150‑1 (MQ=255)
                                                                                                                           ccGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCa                            <  2:3061821/150‑1 (MQ=255)
                                                                                                                           ccGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCa                            <  1:6774042/150‑1 (MQ=255)
                                                                                                                           ccGGCGCTGCCGCCAGGCCCTGGAAGACGTCCAGCGCCAGCGGCCGATACGGGTGTTGATACACGGTGGCAGCCAGTTTCTTGGTGGGGAAATAGAAAGCTGGGCGCGCGAGCATTGCGCCGATCTGGTGCGTTAAACGCCAGACTGGCa                            >  1:4401345/1‑150 (MQ=255)
                                                                                                                              gcgcTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAACTCGAAAGCTGGGCGCGCGGGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCt                         <  2:2768012/150‑1 (MQ=255)
                                                                                                                                    gccgctggCGCTGGACGTCGTCCAGGGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACgg                   >  1:3591103/7‑150 (MQ=255)
                                                                                                                                    ccgccAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGACATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACgg                   >  2:628041/1‑150 (MQ=255)
                                                                                                                                     cgccAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGc                  >  1:3754883/1‑150 (MQ=255)
                                                                                                                                       gctggCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCa                 >  1:6419770/4‑150 (MQ=37)
                                                                                                                                       ccAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCaa                <  1:5996402/150‑1 (MQ=255)
                                                                                                                                        cAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAg               <  1:5594776/150‑1 (MQ=255)
                                                                                                                                           gCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAg            <  2:5411929/150‑1 (MQ=255)
                                                                                                                                            cccTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAgg           <  2:2576775/150‑1 (MQ=255)
                                                                                                                                             gcTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAgg           <  2:8564981/149‑1 (MQ=38)
                                                                                                                                             ccTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGc          >  2:1297095/1‑150 (MQ=255)
                                                                                                                                                 gacgacGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGaa      >  1:1669450/1‑150 (MQ=255)
                                                                                                                                                  acgacgTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAAc     <  1:6330578/150‑1 (MQ=255)
                                                                                                                                                  acgacgTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAAc     <  2:3924934/150‑1 (MQ=255)
                                                                                                                                                  acgacgTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAAc     <  1:7043510/150‑1 (MQ=255)
                                                                                                                                                    gacgTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGc   >  2:4735200/1‑150 (MQ=255)
                                                                                                                                                    gacgTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGc   >  1:2923404/1‑150 (MQ=255)
                                                                                                                                                    gacgTCCAGCGCCAGCGGCCGATCCGGGTGTTGATACACGGTGGCAGCCATTTTCTTGGTGGGGAAATCGAAAGCTGGGAGCGAGAGCATGGCGCCGATCTGGTTAGTTACCCGCCAAAATGGCAGCTGCACGGCAAGCAGGCCGAACg  >  2:4805279/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGCT  >  NC_002947/3038388‑3038686

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