Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,041,135 T→A 7.5% I15F (ATC→TTC) ‡ PP_5293 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,041,1350TA 7.5% 219.1 / 5.3 77I15F (ATC→TTC) ‡PP_5293hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base A (4/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGAAAT  >  NC_002947/6040990‑6041280
                                                                                                                                                 |                                                                                                                                                 
ttCAGCGCCTCTTTGCCCTGGGCGTCGATAAAATGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTc                                                                                                                                                >  1:218117/1‑149 (MQ=255)
 tCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCg                                                                                                                                               <  3:37505/149‑1 (MQ=255)
          ctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGc                                                                                                                                      >  4:48548/1‑149 (MQ=255)
           tttGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCg                                                                                                                                     >  4:35786/1‑149 (MQ=255)
           tttGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCg                                                                                                                                     >  6:267177/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                                                                                                                                      >  2:130661/1‑129 (MQ=255)
                              aaaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGAAATCCGGCAATGGGc                                                                                                                                      <  1:130661/129‑1 (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)
                                                                                                                                             cGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGaaa   >  5:243366/1‑149 (MQ=255)
                                                                                                                                              ggATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGAAAt  >  1:73446/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGAAAT  >  NC_002947/6040990‑6041280

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: ‑