Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,350,501 T→A 9.3% K11N (AAA→AAT) ‡ PP_1176 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,350,5010TA 9.3% 208.6 / 10.8 75K11N (AAA→AAT) ‡PP_1176hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/34);  new base A (4/3);  total (38/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 = 7.60e-01

AGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGT  >  NC_002947/1350357‑1350647
                                                                                                                                                |                                                                                                                                                   
agCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTgg                                                                                                                                                 <  1:178511/149‑1 (MQ=255)
        cTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCAGGAAACCGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCtt                                                                                                                                          <  2:269175/148‑1 (MQ=255)
          gACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgtt                                                                                                                                       >  2:381628/1‑149 (MQ=255)
           aCGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgtt                                                                                                                                       <  4:151729/148‑1 (MQ=255)
           aCGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgt                                                                                                                                        <  2:382721/147‑1 (MQ=255)
             gagaACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGAt                                                                                                                                    >  3:43268/1‑149 (MQ=255)
             gagaACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCtttttatt                                                                                                                                    <  5:388590/149‑4 (MQ=255)
                     gATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGcc                                                                                                                            >  4:31113/1‑149 (MQ=255)
                       ttCGGCTATGTGAGCGTTCTTAAGCGCCGTAAAACGTTATCGCGATGTCCAAAAACACGCAAAGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCg                                                                                                                          <  1:45058/149‑1 (MQ=255)
                       ttCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCg                                                                                                                          >  5:138014/1‑149 (MQ=255)
                          cGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGt                                                                                                                        >  8:55078/1‑148 (MQ=255)
                               agtgAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaat                                                                                                                    <  8:156672/146‑1 (MQ=255)
                                gtgAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaatcaa                                                                                                                >  5:193725/1‑149 (MQ=255)
                                 tgAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTcaatcaat                                                                                                               <  7:152577/149‑1 (MQ=255)
                                      gTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAg                                                                                                          >  3:133453/1‑149 (MQ=255)
                                        tCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCt                                                                                                        >  6:324380/1‑149 (MQ=255)
                                        tCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCt                                                                                                        >  5:340114/1‑149 (MQ=255)
                                           tAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTACCAGGAAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCACAGCTCaa                                                                                                     >  7:395743/1‑149 (MQ=255)
                                            aaGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCaa                                                                                                     <  2:242417/148‑1 (MQ=255)
                                               cgcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGcc                                                                                                 <  8:209092/149‑1 (MQ=255)
                                               cgcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGcc                                                                                                 <  6:340114/149‑1 (MQ=255)
                                               cgcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGAAATAAAAATGCTGACATTTTGCAGCCTTGTGGATGCTCAGACCGGTCAATCAATCCCAGATCAACGcc                                                                                                 >  4:85772/1‑149 (MQ=255)
                                                gacGGAAAACGTTATCGCGATGTCCAAAATCAGGCAATGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGTTCAATCAATCCCAGCTCAACGccc                                                                                                <  4:242345/147‑1 (MQ=255)
                                                  cGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGATTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCa                                                                                              >  5:61663/1‑149 (MQ=255)
                                                   ggAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCAc                                                                                             >  1:76519/1‑149 (MQ=255)
                                                       aaCGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGt                                                                                         <  4:43268/149‑1 (MQ=255)
                                                        aCGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTc                                                                                        <  1:402759/149‑1 (MQ=255)
                                                          gTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTg                                                                                      <  5:324380/149‑1 (MQ=255)
                                                                  gATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGa                                                                              <  6:193725/149‑1 (MQ=255)
                                                                  gATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGa                                                                              <  6:138014/149‑1 (MQ=255)
                                                                       ccAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATcaca                                                                         >  4:102479/1‑149 (MQ=255)
                                                                        cAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCacac                                                                        >  3:396211/1‑149 (MQ=255)
                                                                        cAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCacac                                                                        >  6:137646/1‑149 (MQ=255)
                                                                         aaaaaCAGGCAAGTTCATCAGCCGTAGGAAACGTCCTGTTCTGCTTCTCCAAGGCAAAAAAAAGGCTGCCATTAGGCACCATTGTTGATGCTCAGACCGGTCAATCAATCACAGCTCAACGAAACACCACTCTGATAAGCGATCACAcg                                                                       >  2:185846/1‑149 (MQ=255)
                                                                         aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCCCAcg                                                                       >  3:266108/1‑149 (MQ=255)
                                                                          aaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcgc                                                                      <  1:34728/149‑1 (MQ=255)
                                                                            aaCAGGCTAGGTCGTCTGCCGTAGGAAACGTCCTGTTCTGCTTCTTCCATGCAATAAAAAGGCTGCCATTTTGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCAACCCGTCTGATACTCGATCAGACGCGt                                                                    <  4:133453/149‑1 (MQ=255)
                                                                               atgCAAGGGCGGCAGCCTTATGAAAGGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCACGGTCAATCAATCCCAGCTCATCGCCCCACCCGTCTGATACTCGATCACACGCGtctc                                                                 <  3:85772/147‑1 (MQ=255)
                                                                                ggCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCa                                                                <  3:31113/149‑1 (MQ=255)
                                                                                 gCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCa                                                                >  7:372393/1‑148 (MQ=255)
                                                                                  cAAGGTCGGCAGCCGTACGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATCATCGTTCACGCGCGTCTCaac                                                              >  8:388198/1‑148 (MQ=255)
                                                                                   aaGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCa                                                                >  5:247167/1‑146 (MQ=255)
                                                                                    aGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaaga                                                            <  1:104245/149‑1 (MQ=255)
                                                                                      gTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaagaag                                                          <  3:102479/149‑1 (MQ=255)
                                                                                         ggCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGttc                                                       <  4:396211/149‑1 (MQ=255)
                                                                                             gCCGTAGGAAACGTCCTGTTCTGCGTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTtcttct                                                   <  8:297009/149‑1 (MQ=255)
                                                                                             gCCGTAGGAAACGTCCTGTGCTGCTTCTTCCAGTCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCAGGTCAATCAATCCCATCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTtcttct                                                   <  7:289537/149‑1 (MQ=255)
                                                                                              ccGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTtctc                                                  >  6:347495/1‑149 (MQ=255)
                                                                                              ccGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTtctc                                                  <  7:55078/149‑1 (MQ=255)
                                                                                              ccGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTtctc                                                  >  6:47554/1‑149 (MQ=255)
                                                                                                     aaaCGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCttcttc                                           <  5:137646/149‑1 (MQ=255)
                                                                                                         gTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGt                                       <  2:76519/149‑1 (MQ=255)
                                                                                                            cTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTcc                                     <  5:47554/148‑1 (MQ=255)
                                                                                                              gTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCAt                                   >  2:113669/1‑148 (MQ=255)
                                                                                                               ttctgcttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGa                                 >  2:34442/1‑149 (MQ=255)
                                                                                                                    cttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGAtct                              <  6:247167/147‑1 (MQ=255)
                                                                                                                           cAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATc                      >  8:143045/1‑148 (MQ=255)
                                                                                                                              gCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATg                  >  2:422140/1‑149 (MQ=255)
                                                                                                                              gCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATg                  >  3:73128/1‑149 (MQ=255)
                                                                                                                              gCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATg                  >  1:400931/1‑149 (MQ=255)
                                                                                                                                   aacaaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATg                                  >  3:138138/4‑128 (MQ=255)
                                                                                                                                   aacaaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATg                                  <  4:138138/125‑1 (MQ=255)
                                                                                                                                   aaaaaGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCCAATCCCAACTCAACGCCCCACCCGTc                                                                                        <  1:226648/75‑1 (MQ=255)
                                                                                                                                   aaaaaGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAAc             >  3:183374/1‑149 (MQ=255)
                                                                                                                                   aaaaaGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAACTCAACGCCCCACCCGTc                                                                                        >  2:226648/1‑74 (MQ=255)
                                                                                                                                    acaaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCAATGATCTCGCTCATCCATGGGAACg            >  6:262740/3‑149 (MQ=255)
                                                                                                                                     caaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACgg           >  6:423862/2‑149 (MQ=255)
                                                                                                                                     caaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACgg           >  5:131765/2‑149 (MQ=255)
                                                                                                                                     aaaGGCTGCCATTTGACAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACgg           >  4:356208/1‑149 (MQ=255)
                                                                                                                                             ccATTTTGCAGCCTTGTTGATGCTCAGCCCGGTCAATCACTCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTa    <  8:395743/148‑1 (MQ=255)
                                                                                                                                             ccATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTa    <  6:61663/148‑1 (MQ=255)
                                                                                                                                             ccATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTCCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTa    <  1:113669/148‑1 (MQ=255)
                                                                                                                                              cATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCGACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATgg                 <  2:179923/134‑1 (MQ=255)
                                                                                                                                              cATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATgg                 >  1:179923/1‑134 (MQ=255)
                                                                                                                                              cATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGAGCACACGCGTCTCAAACAAGTTGTTCTCGTGCTTCAAGTCCATGATCTGGCTCATCCATGCGAACAGATTGTTAgg  >  8:361466/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGT  >  NC_002947/1350357‑1350647

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