Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,041,127 T→A 6.6% Q17H (CAA→CAT)  PP_5293 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,041,1270TA 6.6% 229.2 / 4.0 76Q17H (CAA→CAT) PP_5293hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base A (3/2);  total (40/36)
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.66e-01

TGAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTT  >  NC_002947/6040984‑6041275
                                                                                                                                               |                                                                                                                                                    
tGAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCAATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCg                                                                                                                                                 >  6:52906/1‑149 (MQ=255)
  aaaGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTg                                                                                                                                                    <  3:94712/144‑1 (MQ=255)
      ttCAGCGCCTCTTTGCCCTGGGCGTCGATAAAATGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTc                                                                                                                                           >  1:218117/1‑149 (MQ=255)
       tCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCg                                                                                                                                          <  3:37505/149‑1 (MQ=255)
                ctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGc                                                                                                                                 >  4:48548/1‑149 (MQ=255)
                 tttGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCg                                                                                                                                >  6:267177/1‑149 (MQ=255)
                 tttGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCg                                                                                                                                >  4:35786/1‑149 (MQ=255)
                    gCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCggg                                                                                                                              <  2:186228/148‑1 (MQ=255)
                     cccTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCgggg                                                                                                                             <  2:218117/148‑1 (MQ=255)
                     cccTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCgggg                                                                                                                             <  7:79933/148‑1 (MQ=255)
                     cccTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGt                                                                                                                            >  7:302606/1‑149 (MQ=255)
                          ggCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATc                                                                                                                        <  8:189053/148‑1 (MQ=255)
                             gTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGa                                                                                                                      >  8:69727/1‑147 (MQ=255)
                                 aTAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCACGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGtt                                                                                                                >  7:9551/1‑149 (MQ=255)
                                 aTAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGtt                                                                                                                >  7:72204/1‑149 (MQ=255)
                                  tAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTc                                                                                                               >  7:290882/1‑149 (MQ=255)
                                    aaaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGAAATCCGGCAATGGGc                                                                                                                                 <  1:130661/129‑1 (MQ=255)
                                    aaaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGAAATCCGGCAATGGGc                                                                                                                                 >  2:130661/1‑129 (MQ=255)
                                       gggTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATc                                                                                                          <  3:35786/149‑1 (MQ=255)
                                              gCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTc                                                                                                      <  5:52906/146‑1 (MQ=255)
                                              gCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCACTc                                                                                                      <  3:54241/146‑1 (MQ=255)
                                               cTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGc                                                                                                  <  8:302606/149‑1 (MQ=255)
                                                tCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGcc                                                                                                 >  3:162210/1‑149 (MQ=255)
                                                        gAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcc                                                                                          <  5:342151/148‑1 (MQ=255)
                                                         aGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctc                                                                                        <  7:69727/149‑1 (MQ=255)
                                                         aGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctc                                                                                        >  3:244621/1‑149 (MQ=255)
                                                          ggTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctcc                                                                                       >  5:161165/1‑149 (MQ=255)
                                                             gCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTg                                                                                    >  3:1913/1‑149 (MQ=255)
                                                              ccATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTgg                                                                                   >  2:27524/1‑149 (MQ=255)
                                                                aTAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGcc                                                                                 >  1:104546/1‑149 (MQ=255)
                                                                 tAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCg                                                                                <  4:162210/149‑1 (MQ=255)
                                                                 tAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCg                                                                                <  4:1777/149‑1 (MQ=255)
                                                                  aGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCgg                                                                               >  5:306127/1‑149 (MQ=255)
                                                                   gCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGt                                                                              <  7:47094/149‑1 (MQ=255)
                                                                    cACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGtt                                                                             >  1:327412/1‑149 (MQ=255)
                                                                     aCCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGttt                                                                            <  6:253246/149‑1 (MQ=255)
                                                                     aCCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGttt                                                                            <  2:327412/149‑1 (MQ=255)
                                                                        aGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtcag                                                                         <  6:161165/149‑1 (MQ=255)
                                                                           gggCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtcagtca                                                                      <  8:9551/149‑1 (MQ=255)
                                                                             gCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGACAGc                                                                    <  3:48548/149‑1 (MQ=255)
                                                                               ccAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCa                                                                  <  1:135400/149‑1 (MQ=255)
                                                                                  gcagGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtg                                                               <  4:244621/149‑1 (MQ=255)
                                                                                     ggTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtgctg                                                            <  5:267177/149‑1 (MQ=255)
                                                                                      gTGCCCAAGCGGCACTATCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctcct                                                                                      >  8:262417/1‑122 (MQ=255)
                                                                                         cccAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCTGTTTCAGTCAGCCAGTGctgcctg                                                        >  8:275885/1‑149 (MQ=255)
                                                                                          ccAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGctgcctgc                                                       <  1:147798/149‑1 (MQ=255)
                                                                                           cAAGCGGCACTACCAGGCCCAACACCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGAGCGCCAGGTTCGAATCGCTCAGGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTCCctgct                                                      >  5:347048/1‑149 (MQ=255)
                                                                                                   aCTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCGATTGCCGGAATTAGGAATTGAGCCTGGTGATCGACTGGTTCGAATCGCTCATGCACCGCACCTCCTTGGCCGGTTTCAGTCAGACAGTGCTGCCTGCTGCAGCTc                                               >  5:238604/1‑148 (MQ=255)
                                                                                                     tACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGa                                             >  7:185643/1‑148 (MQ=255)
                                                                                                       ccAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGCATCGCTCATGCCCCGCACCTCCTTGGCGGGTTTCTGTAAGCCAGTGCTGCCTGCTGCAGCTCGaaaa                                          >  4:213593/1‑149 (MQ=255)
                                                                                                         aGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGTTGCAGCTCGAAAAtc                                        <  8:126588/149‑1 (MQ=255)
                                                                                                           gCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtctc                                      <  6:306127/149‑1 (MQ=255)
                                                                                                            cccAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCg                                     <  2:104546/149‑1 (MQ=255)
                                                                                                             ccAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtctctt                                    >  2:66720/1‑147 (MQ=255)
                                                                                                             ccAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGt                                    >  7:100205/1‑149 (MQ=255)
                                                                                                              cAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGt                                    >  7:219136/1‑148 (MQ=255)
                                                                                                              cAACAGCGCGGAGAGGTGGCAGAACATCGCACATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGtt                                   <  8:185643/149‑1 (MQ=255)
                                                                                                                aCAGCGCGGAGAGGTGGCAGAACATCGCCAATTGCCGGATTTCGGCATTGGGCGGGGTGATCTACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTTGCCGGTTTCAGTCAGCCATTGCTGCCAGCTGCAGCTCGAAAATCTCGtac                                  >  6:4789/1‑146 (MQ=255)
                                                                                                                 cAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCAt                                <  8:100205/149‑1 (MQ=255)
                                                                                                                  aGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATg                               <  4:1913/149‑1 (MQ=255)
                                                                                                                   gcgcGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCAGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGc                              >  4:2745/1‑149 (MQ=255)
                                                                                                                             gTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtg                                                               >  7:26683/1‑106 (MQ=255)
                                                                                                                             gTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtg                                                               <  8:26683/106‑1 (MQ=255)
                                                                                                                                gCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtca                                                                          >  2:237774/1‑92 (MQ=255)
                                                                                                                                  aGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACGGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAg               >  6:147581/1‑149 (MQ=255)
                                                                                                                                   gAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGc              >  4:12377/1‑149 (MQ=255)
                                                                                                                                       atctcgcATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCtt                          <  6:93459/127‑1 (MQ=255)
                                                                                                                                       aTCGGCCATTGCCGGATTTCGTCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCGTTGGCCGGTTGCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGc          <  6:347048/149‑1 (MQ=255)
                                                                                                                                         cGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCAt        >  4:109324/1‑149 (MQ=255)
                                                                                                                                           cccATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGt                                    <  2:49695/119‑1 (MQ=255)
                                                                                                                                           cccATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGt                                    >  1:49695/1‑119 (MQ=255)
                                                                                                                                            ccAATGCCGAAATCCGGCAATGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctcct                                                                                      <  5:189878/68‑1 (MQ=25)
                                                                                                                                            ccAATGCCGAAATCCGGCAATGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctcct                                                                                      >  6:189878/1‑68 (MQ=25)
                                                                                                                                            ccAATGCCGAAATCCGGCAATGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGc     >  1:166207/1‑149 (MQ=255)
                                                                                                                                             cATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCg    >  4:206197/1‑149 (MQ=255)
                                                                                                                                              aTTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGt   >  8:22339/1‑149 (MQ=255)
                                                                                                                                               ttGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGtt  <  2:6059/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTT  >  NC_002947/6040984‑6041275

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