Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,668,348 T→G 6.0% K113N (AAA→AAC)  PP_3231 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,668,3480TG 6.0% 206.1 / 3.9 67K113N (AAA→AAC) PP_3231hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base G (2/2);  total (33/34)
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.78e-01

ACATCCAGGCCGCGGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGC  >  NC_002947/3668217‑3668493
                                                                                                                                    |                                                                                                                                                 
acaTCCAGGCCGCGGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATA‑GGGCGAAGGCTTTGACTTCTCTATTGCGCCACATCCAGGTAGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTgcgc                                                                                                                                  >  8:182788/1‑149 (MQ=255)
acaTCCAGGCCGCGGAATTCCAAGCGGTGGTCGGTGGGCAGCACGCGGTTGAGACCGAGCAGCCATT‑GGGCGAAGGATTTGCCTTATCTATTGCGCCACCTCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCAAGACGAGTgccc                                                                                                                                  >  1:72984/1‑147 (MQ=255)
   tCCAGGCCGCGGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATGCAGGTTGGATAACGGCTGATAATATGGCGTTTCCACGCGGAGTgcgcg                                                                                                                                 >  4:55782/1‑146 (MQ=255)
    ccAGGCCGCGGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAACCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAg                                                                                                                              >  3:94945/1‑149 (MQ=255)
      aGGCCGCGGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGcc                                                                                                                            >  1:182404/1‑149 (MQ=255)
        gCCGCGGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTc                                                                                                                          <  5:106206/149‑1 (MQ=255)
         ccGCGGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCg                                                                                                                         <  6:356309/149‑1 (MQ=255)
           gcgGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAATATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCgca                                                                                                                       >  2:363441/1‑147 (MQ=255)
            cgGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGTTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAATCCTCGGAc                                                                                                                      <  6:208065/149‑1 (MQ=255)
            cgGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGAc                                                                                                                      >  1:195418/1‑149 (MQ=255)
            cgGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGAc                                                                                                                      <  6:167191/149‑1 (MQ=255)
            cgGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGAc                                                                                                                      <  2:146733/149‑1 (MQ=255)
                aTTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTAc                                                                                                                  <  1:176755/149‑1 (MQ=255)
                aTTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTAGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTa                                                                                                                   <  4:152859/148‑1 (MQ=255)
                 ttCCACGGGGTGGTCGGTGGGCAGCTGGCGGTTGTGCCCGTGCAGCCATT‑GGTCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTAc                                                                                                                  <  7:241213/148‑1 (MQ=255)
                     aCGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATc                                                                                                             >  6:246860/1‑149 (MQ=255)
                       gcgGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGa                                                                                                           <  6:291859/149‑1 (MQ=255)
                       gcgGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGa                                                                                                           <  3:9252/149‑1 (MQ=255)
                          gtggtCGGTGGGCTGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACAt                                                                                                        <  1:134589/149‑1 (MQ=255)
                          gtggtCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACAt                                                                                                        >  5:198781/1‑149 (MQ=255)
                                ggTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCgg                                                                                                  >  6:77178/1‑149 (MQ=255)
                                       agcagGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGGGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCgg                                                                                           <  1:363441/149‑1 (MQ=255)
                                       agcagGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCgg                                                                                           <  5:246860/149‑1 (MQ=255)
                                          agGCGGTTGTGCCCGTGCAGCCATTGGGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGc                                                                                          >  1:4285/1‑148 (MQ=255)
                                            gCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGt                                                                                      >  4:2623/1‑149 (MQ=255)
                                                  gtgCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGcc                                                                                 <  6:154457/148‑1 (MQ=255)
                                                   tgCCCTTGAAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTc                                                                               <  4:200211/149‑1 (MQ=255)
                                                      ccGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGAc                                                                            >  3:203115/1‑149 (MQ=255)
                                                              ccATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTAc                                                                                                                  >  6:308668/1‑103 (MQ=255)
                                                              ccATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTAc                                                                                                                  <  5:308668/103‑1 (MQ=255)
                                                              ccATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGcgacga                                                                    >  7:109029/1‑149 (MQ=255)
                                                                 tt‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGt                                                                 <  2:342225/149‑1 (MQ=255)
                                                                      gCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTg                                                               >  8:113383/1‑147 (MQ=255)
                                                                      gCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGaa                                                             <  1:124431/149‑1 (MQ=255)
                                                                      gCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGaa                                                             <  3:340789/149‑1 (MQ=255)
                                                                         aaGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGcc                                                          <  8:109029/149‑1 (MQ=255)
                                                                            gCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATg                                                       <  4:94945/149‑1 (MQ=255)
                                                                              tttGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTa                                                                                                                   <  8:128061/87‑1 (MQ=255)
                                                                              tttGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTa                                                                                                                   >  7:128061/1‑87 (MQ=255)
                                                                              tttGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCt                                                     >  6:198550/1‑149 (MQ=255)
                                                                               ttGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTg                                                    >  2:363629/1‑149 (MQ=255)
                                                                                  aCTTCGCTATTGCGCCACATCCATGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGGCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTc                                                 <  7:182788/149‑1 (MQ=255)
                                                                                  aCTTCGATATTTCGCCACATCCAGGTTGGAAAACGGCTTATAATATGGCGTTTCCAGGGCGAGTGCGCCAAGCCTCGGACGTACTGATCGAGATGCCCGGCATCCGGCAAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGtt                                                  <  6:229422/148‑1 (MQ=255)
                                                                                    ttCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGa                                               >  6:333813/1‑149 (MQ=255)
                                                                                    ttCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGa                                               >  5:129056/1‑149 (MQ=255)
                                                                                    ttCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGGTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATg                                                                  <  8:360086/130‑1 (MQ=255)
                                                                                    ttCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGGTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATg                                                                  <  8:355286/130‑1 (MQ=255)
                                                                                    ttCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGGTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATg                                                                  >  7:355286/1‑130 (MQ=255)
                                                                                    ttCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGGTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATg                                                                  >  7:360086/1‑130 (MQ=255)
                                                                                     tCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGAt                                              >  5:116653/1‑149 (MQ=255)
                                                                                       gCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATc                                             <  4:203115/148‑1 (MQ=255)
                                                                                         tATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAgg                                          >  6:333819/1‑149 (MQ=255)
                                                                                                ccACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCggg                                   >  8:267750/1‑149 (MQ=255)
                                                                                                       cAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCAGGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGc                            <  2:72984/149‑1 (MQ=255)
                                                                                                       cAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCg                                                                                                   <  1:59408/78‑1 (MQ=255)
                                                                                                       cAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCg                                                                                                   >  2:59408/1‑78 (MQ=255)
                                                                                                          gTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCAc                         >  3:65577/1‑149 (MQ=255)
                                                                                                             ggAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCg                      >  6:58864/1‑149 (MQ=255)
                                                                                                                aaaCGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGgc                   <  5:77178/149‑1 (MQ=255)
                                                                                                                aaaCGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGgc                   <  2:182404/149‑1 (MQ=255)
                                                                                                                aaaCGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATCCCCGGCATCCCGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGgc                   <  3:55782/149‑1 (MQ=255)
                                                                                                                     gCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGAt              <  7:267750/149‑1 (MQ=255)
                                                                                                                       tGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGa               >  4:10786/1‑146 (MQ=255)
                                                                                                                         aaaaTATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGa           <  7:113383/148‑1 (MQ=255)
                                                                                                                           aaTATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACt        >  4:286993/1‑149 (MQ=255)
                                                                                                                              atGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCgc     >  2:94631/1‑149 (MQ=255)
                                                                                                                                 gCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCgctgc  >  6:172824/1‑149 (MQ=255)
                                                                                                                                    |                                                                                                                                                 
ACATCCAGGCCGCGGAATTCCACGCGGTGGTCGGTGGGCAGCAGGCGGTTGTGCCCGTGCAGCCATT‑GGGCGAAGGCTTTGACTTCGCTATTGCGCCACATCCAGGTTGGAAAACGGCTGAAAATATGGCGTTTCCACGCCGAGTGCGCCAAGCCTCGGACGTACTGATCGACATGCCCGGCATCCGGCCAGTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGC  >  NC_002947/3668217‑3668493

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