Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,624,250 T→A 6.4% E413V (GAA→GTA)  PP_3195 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,624,2500TA 6.4% 189.3 / 2.7 62E413V (GAA→GTA) PP_3195hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base A (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAGGTGCGGTGTGCCCATTTGA  >  NC_002947/3624104‑3624397
                                                                                                                                                  |                                                                                                                                                   
cccGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCg                                                                                                                                                   <  3:295348/149‑1 (MQ=255)
 ccGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCgg                                                                                                                                                  >  2:251007/1‑149 (MQ=255)
      aCCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCg                                                                                                                                                   <  4:89955/143‑1 (MQ=255)
      aCCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCg                                                                                                                                                   >  3:89955/1‑143 (MQ=255)
               tGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGc                                                                                                                                     >  4:340677/1‑148 (MQ=255)
                     ttaTGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATGGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTa                                                                                                                               <  6:27285/146‑1 (MQ=255)
                           ccAGGTATCGATGGCATCGATACGTGGCGCTCTGTGGCGCAGAAACGATAACGGCGACAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCg                                                                                                                        <  2:81112/149‑1 (MQ=255)
                           ccAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTACGGCGATCCAGTCGGCGATGTAGTGCCCg                                                                                                                        >  7:262369/1‑149 (MQ=255)
                           caaGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCTGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTACGGCGATCCAGTCGGCGATGTAGTGCCCg                                                                                                                        >  7:262366/3‑149 (MQ=255)
                            cAGGTATCGCTGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCgg                                                                                                                       >  7:189961/1‑149 (MQ=255)
                             aGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGCGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCgg                                                                                                                       >  6:384353/1‑148 (MQ=255)
                              ggTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCa                                                                                                                     <  7:65297/149‑1 (MQ=255)
                              ggTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCa                                                                                                                     >  4:17362/1‑149 (MQ=255)
                                    gATGGCGTCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCAGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATg                                                                                                               >  1:308194/1‑149 (MQ=255)
                                    gATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATg                                                                                                               >  3:255048/1‑149 (MQ=255)
                                      tGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCt                                                                                                             >  1:403895/1‑149 (MQ=255)
                                         cATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGc                                                                                                              >  1:407762/1‑145 (MQ=255)
                                         cATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGc                                                                                                              <  2:407762/145‑1 (MQ=255)
                                          aTCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGc                                                                                                              >  2:81717/1‑144 (MQ=255)
                                          aTCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGc                                                                                                              <  1:81717/144‑1 (MQ=255)
                                          aTCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTa                                                                                                         <  3:17362/149‑1 (MQ=255)
                                          aTCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTa                                                                                                         >  5:374624/1‑149 (MQ=255)
                                          aTCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTa                                                                                                         <  1:251007/149‑1 (MQ=255)
                                           tCGATACCTGGCGCTCTCTGGCGCAGAAACGAGACCGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCGGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTAt                                                                                                        <  5:285477/149‑1 (MQ=255)
                                              aTACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGa                                                                                                     <  4:255048/149‑1 (MQ=255)
                                                 ccTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAg                                                                                                   <  6:125015/148‑1 (MQ=255)
                                                  cTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCt                                                                                                 >  1:50915/1‑149 (MQ=255)
                                                   tGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTc                                                                                                >  2:257319/1‑149 (MQ=255)
                                                     gcgcTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTc                                                                                                <  6:8821/147‑1 (MQ=255)
                                                     gcgcTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTc                                                                                                <  3:340677/147‑1 (MQ=255)
                                                           ggtGGGGCATAAACGTCAACGGCGCCATTCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCTAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAAtgctg                                                                                        <  2:308194/147‑1 (MQ=255)
                                                               gcgcAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCAc                                                                                     <  2:174619/148‑1 (MQ=255)
                                                                 gcAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAAt                                                                                            <  4:373906/139‑1 (MQ=255)
                                                                 gcAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAAt                                                                                            >  3:373906/1‑139 (MQ=255)
                                                                 gcAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGcc                                                                                  >  3:229088/1‑149 (MQ=255)
                                                                  cAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCa                                                                                 >  3:121698/1‑149 (MQ=255)
                                                                     aaaCGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTACCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCt                                                                              >  5:330233/1‑149 (MQ=255)
                                                                         ttgggcGGCGCCGATTGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCAc                                                                          >  1:172091/6‑149 (MQ=255)
                                                                            aaCGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGt                                                                       >  6:85735/1‑149 (MQ=255)
                                                                                    cAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAAAGCTGCACGCCACCTTCACGGTCAGTGaa                                                                >  5:190536/1‑148 (MQ=255)
                                                                                                      ccAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTc                                             <  5:384353/149‑1 (MQ=255)
                                                                                                       cAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTACCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCg                                            <  6:330233/149‑1 (MQ=255)
                                                                                                           gCCGGCAGCGCCTATCAGGCCGTGTGGCCGCCAAGGGTTACGGCGATCCAGTCGGCGATGTAGTGCCCGGCATCGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCt                                        <  8:262369/149‑1 (MQ=255)
                                                                                                           gCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCt                                        <  8:189961/149‑1 (MQ=255)
                                                                                                           gCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTACGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCt                                        <  8:262366/149‑1 (MQ=255)
                                                                                                             cGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGt                                      <  3:353542/149‑1 (MQ=255)
                                                                                                                    gcCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGtt                                     >  8:222604/1‑143 (MQ=255)
                                                                                                                    gcCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGtt                                     <  7:222604/143‑1 (MQ=255)
                                                                                                                       gATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCa                                                                     <  6:334017/108‑1 (MQ=255)
                                                                                                                         tCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCt                          >  6:204430/1‑149 (MQ=255)
                                                                                                                             gCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCa                                                        >  8:36602/1‑115 (MQ=255)
                                                                                                                             gCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCa                                                        <  7:36602/115‑1 (MQ=255)
                                                                                                                             gCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCg                        >  4:158257/1‑147 (MQ=255)
                                                                                                                             gCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTa                      <  2:172091/149‑1 (MQ=255)
                                                                                                                              ccGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGTAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAg                     <  6:374624/149‑1 (MQ=255)
                                                                                                                                gtgtGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAGGt                   >  5:294230/1‑149 (MQ=255)
                                                                                                                                   tgGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAGGTGCg                >  6:95177/1‑149 (MQ=255)
                                                                                                                                       cgccAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAGGTGCGgtgt            <  2:403895/149‑1 (MQ=255)
                                                                                                                                        gccAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAGGTGCGGtgtg           <  4:121698/149‑1 (MQ=255)
                                                                                                                                          cAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAGGTGCGGTGTGcc         >  5:413527/1‑149 (MQ=255)
                                                                                                                                            aGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAGGTGCGGTGTGCCCa       >  1:45592/1‑149 (MQ=255)
                                                                                                                                             gggTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAGGTGCGGTGTGCCCAt      <  4:229088/149‑1 (MQ=255)
                                                                                                                                                tttCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCggg                                          <  8:88582/110‑1 (MQ=255)
                                                                                                                                                tttCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCggg                                          >  7:88582/1‑110 (MQ=255)
                                                                                                                                                  tCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAGGTGCGGTGTGCCCATTTGa  <  5:85735/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTTATCGAAGCTCGAATGCTGCACGCCACCTTCACGGTCAGTGAAAATCTTCAGTTCGCGTTTCGGGCTGTTGACCAGTTGCTCGTAGGTGCGGTGTGCCCATTTGA  >  NC_002947/3624104‑3624397

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