Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 563,005 C→A 9.5% P30T (CCT→ACT) ‡ rpoA → DNA‑directed RNA polymerase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947563,0050CA 9.5% 216.7 / 11.9 74P30T (CCT→ACT) ‡rpoADNA‑directed RNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (37/30);  new base A (4/3);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

ATCAACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTAA  >  NC_002947/562857‑563152
                                                                                                                                                     |                                                                                                                                                   
aTCAACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGc                                                                                                                                                     <  6:203078/149‑1 (MQ=255)
    aCGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCctct                                                                                                                                                 >  3:224944/1‑149 (MQ=255)
       aaaaCCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCtctc                                                                                                                                                >  8:346319/1‑147 (MQ=255)
        aaaCCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGa                                                                                                                                              <  8:165816/148‑1 (MQ=255)
              gATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGt                                                                                                                                       >  5:4445/1‑149 (MQ=255)
                   tCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCTGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATCACGCTCGAGAGGCTCGAGCGtaattttg                                                                                                                                   <  3:288392/148‑8 (MQ=255)
                    cGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGc                                                                                                                                 >  1:267200/1‑149 (MQ=255)
                     gAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGc                                                                                                                                 >  1:31648/1‑148 (MQ=255)
                       gCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGc                                                                                                                                 >  5:34702/1‑146 (MQ=255)
                       gCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGc                                                                                                                                 <  6:34702/146‑1 (MQ=255)
                        ctctACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATa                                                                                                                             >  6:187317/1‑149 (MQ=255)
                         tctACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATAc                                                                                                                            >  1:364316/1‑149 (MQ=255)
                           tACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCGAAACCACGCTCGAGAGGCTCGAGCGTAATTTTGGcgc                                                                                                                               >  3:123919/1‑142 (MQ=37)
                           tACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCGAAACCACGCTCGAGAGGCTCGAGCGTAATTTTGGcgc                                                                                                                               <  4:123919/144‑3 (MQ=37)
                           tACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATAccc                                                                                                                          >  7:264056/1‑149 (MQ=255)
                                cAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTggg                                                                                                                      >  8:44172/1‑148 (MQ=255)
                                cAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTgg                                                                                                                       <  7:346319/147‑1 (MQ=255)
                                   gTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGCCCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAAc                                                                                                                  >  6:201351/1‑149 (MQ=255)
                                       gggCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACgcgc                                                                                                              >  1:118018/1‑149 (MQ=255)
                                         gCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTg                                                                                                            >  5:303297/1‑149 (MQ=255)
                                          cTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGAGGCTCGAGCg                                                                                                                                           <  7:245157/117‑1 (MQ=255)
                                          cTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGAGGCTCGAGCg                                                                                                                                           >  8:245157/1‑117 (MQ=255)
                                           tAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTgcg                                                                                                          >  5:117651/1‑149 (MQ=255)
                                             gAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGcgcc                                                                                                        <  4:45958/149‑1 (MQ=255)
                                                aaTAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCa                                                                                                     >  5:219720/1‑149 (MQ=255)
                                                     gggCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtg                                                                                                <  4:224944/147‑1 (MQ=255)
                                                         aTCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttgt                                                                                            >  4:88118/1‑149 (MQ=255)
                                                          tCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtc                                                                                           >  3:56387/1‑149 (MQ=255)
                                                          tCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtc                                                                                           <  6:303297/149‑1 (MQ=255)
                                                            cATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtcct                                                                                         <  2:267200/149‑1 (MQ=255)
                                                             aTGCAGATTTCGGTAAATGAGTTCCTGGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTATCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtcct                                                                                         >  1:204281/1‑149 (MQ=255)
                                                             aTGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtcctc                                                                                        <  2:364316/149‑1 (MQ=255)
                                                              tGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtcctc                                                                                        >  6:189321/1‑148 (MQ=255)
                                                                  gATTTCGGTAAATGAGTACCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGc                                                                                   >  4:336347/1‑149 (MQ=255)
                                                                         gTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACCCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCtg                                                                            >  3:198271/1‑149 (MQ=255)
                                                                          tAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCtgt                                                                           <  4:56387/149‑1 (MQ=255)
                                                                            aaTGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCtgtg                                                                          <  5:187317/148‑1 (MQ=255)
                                                                             aTGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCtgtg                                                                          >  1:104332/1‑147 (MQ=255)
                                                                                 gTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCagta                                                                     <  7:25920/148‑1 (MQ=255)
                                                                                  ttCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCagtagt                                                                   >  8:8043/1‑149 (MQ=255)
                                                                                         aCTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGc                                                             <  6:117651/149‑1 (MQ=255)
                                                                                          cTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGcc                                                            <  5:189321/149‑1 (MQ=255)
                                                                                          cTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGcc                                                            <  8:264056/149‑1 (MQ=255)
                                                                                           tCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGcc                                                            >  4:199293/1‑148 (MQ=255)
                                                                                             cccGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGc                                                                                   <  7:333734/123‑1 (MQ=255)
                                                                                             cccGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGc                                                                                   >  8:333734/1‑123 (MQ=255)
                                                                                             cccGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCg                                                           <  6:243641/147‑1 (MQ=255)
                                                                                              ccGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGcc                                                            >  6:263135/1‑145 (MQ=255)
                                                                                              ccGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCgaga                                                        <  2:104332/149‑1 (MQ=255)
                                                                                                   acaTTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGAt                                                   <  4:198271/149‑1 (MQ=255)
                                                                                                    caTTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATg                                                  <  6:219720/149‑1 (MQ=255)
                                                                                                     aTTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGa                                                                <  2:28054/134‑1 (MQ=255)
                                                                                                        gATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGAGCAGTAGTCGAGGCCGAGATCGATGGCg                                               >  3:16489/1‑148 (MQ=255)
                                                                                                         aTGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTa                                             >  1:67261/1‑149 (MQ=255)
                                                                                                          tgtgCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTAc                                            >  5:30025/1‑149 (MQ=255)
                                                                                                          tgtgCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTAc                                            >  6:210040/1‑149 (MQ=255)
                                                                                                              cAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTcc                                         <  3:88118/148‑1 (MQ=255)
                                                                                                                    gtcagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCg                                                                 >  4:140111/1‑118 (MQ=255)
                                                                                                                    gtcagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCg                                                                 <  3:140111/118‑1 (MQ=255)
                                                                                                                     tcagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTa                                  <  4:16489/148‑1 (MQ=255)
                                                                                                                       agtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTc                               >  7:132837/1‑149 (MQ=255)
                                                                                                                           cGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGTAGTAGTTGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCcgcg                           >  4:8275/1‑149 (MQ=255)
                                                                                                                           cGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCcgc                            <  6:4445/148‑1 (MQ=255)
                                                                                                                               ccGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCgaga                                                        <  5:342969/116‑1 (MQ=255)
                                                                                                                               ccGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCgaga                                                        >  6:342969/1‑116 (MQ=255)
                                                                                                                                 gcgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGaa                     <  7:44172/149‑1 (MQ=255)
                                                                                                                                 gcgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGaa                     <  5:201351/149‑1 (MQ=255)
                                                                                                                                  cgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                    <  5:263135/149‑1 (MQ=255)
                                                                                                                                      aaaaTTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgtgt                <  2:67261/149‑1 (MQ=255)
                                                                                                                                        aaTTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGtt               >  1:258003/1‑148 (MQ=255)
                                                                                                                                           caCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAg            >  3:198391/2‑148 (MQ=255)
                                                                                                                                           caCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAg            >  3:203065/2‑148 (MQ=255)
                                                                                                                                                tCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTc                               <  2:212704/124‑1 (MQ=255)
                                                                                                                                                tCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGAc      >  6:266865/1‑149 (MQ=255)
                                                                                                                                                    gCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTaa  >  4:197/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATCAACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCT‑GACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTAA  >  NC_002947/562857‑563152

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