Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 77,503 A→T 6.2% V14E (GTG→GAG)  rsmB ← ribosomal RNA small subunit methyltransferase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294777,5030AT 6.2% 215.7 / 3.0 74V14E (GTG→GAG) rsmBribosomal RNA small subunit methyltransferase B
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base T (2/3);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGGGC  >  NC_002947/77358‑77650
                                                                                                                                                  |                                                                                                                                                   
cagcagCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAg                                                                                                                                                  >  3:237387/1‑149 (MQ=255)
 agcagcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGc                                                                                                                                                 >  2:90555/1‑149 (MQ=255)
  gcagcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                                >  8:17189/1‑149 (MQ=255)
   cagcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                                <  1:322267/148‑1 (MQ=255)
   cagcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                                >  2:322267/1‑148 (MQ=255)
       tgtgCGGCCAGCAGGTCGAGTCGTGGCTGCGATCGTTCGGTACCGTACGCGAGGTCCTGGTTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCATCGAAGCCTTGGCGCTGAGCACAGCGGCCAg                                                                                                                                           <  8:35230/149‑1 (MQ=255)
              ccAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCCCgg                                                                                                                                    >  1:171713/1‑149 (MQ=255)
                agcagGTAGAGGCGTGGCTGCCAGCGTGCCGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTGATCTACATTTTACAGTTGCGCCCGCAGCGAGCTGTTCAGCGAAGCCATGCCCCAGAACCCAAAGGCCAGCGCACCggc                                                                                                                                  >  2:302169/1‑149 (MQ=255)
                         aGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAg                                                                                                                         <  5:37370/149‑1 (MQ=255)
                           gCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGc                                                                                                                       >  1:265103/1‑149 (MQ=255)
                                 cTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGTGCTCAGCGGCAAg                                                                                                                                           <  8:288363/123‑1 (MQ=255)
                                 cTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGTGCTCAGCGGCAAg                                                                                                                                           >  7:288363/1‑123 (MQ=255)
                                   gCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGtt                                                                                                                >  2:132165/1‑148 (MQ=255)
                                       gCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTg                                                                                                           <  3:296001/149‑1 (MQ=255)
                                        cGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTgg                                                                                                          <  1:132165/149‑1 (MQ=255)
                                           gCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGATCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGcc                                                                                                        >  1:63087/1‑149 (MQ=255)
                                            cGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGccc                                                                                                       <  1:90555/148‑1 (MQ=255)
                                              gTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCTCGTTCATCTACCTTGTCCAGTTGCTCCGGCAGCGAGCTGTTCAGCGAATCCTTGCCGCTGAGCACAGCGGCCAGCGCACGTTCGGCGGCGAGGCGTGGTTTCATTGGCCCAGc                                                                                                    <  1:302169/149‑1 (MQ=255)
                                               tACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCa                                                                                                   >  5:311872/1‑149 (MQ=255)
                                                 ccGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCAcc                                                                                                 <  6:220088/149‑1 (MQ=255)
                                                 ccGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCAc                                                                                                  <  7:358904/148‑1 (MQ=255)
                                                  cGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCt                                                                                                >  5:285081/1‑149 (MQ=255)
                                                   gAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCtt                                                                                               >  6:185696/1‑149 (MQ=255)
                                                         cAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCcggc                                                                                         >  3:216023/1‑149 (MQ=255)
                                                             tCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCggcggc                                                                                      >  5:325582/1‑148 (MQ=255)
                                                                tGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                                  >  2:307267/1‑149 (MQ=255)
                                                                        gCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGTTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTcgcg                                                                          >  1:76935/1‑149 (MQ=255)
                                                                        gCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGggcggcggc                                                                                                                            <  6:284225/99‑1 (MQ=255)
                                                                        gCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGggcggcggc                                                                                                                            >  5:284225/1‑99 (MQ=255)
                                                                        gCCGCGGTCGCGTTAATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGCCCAGCGAACGGGCGGCGGCCAGGCGTGGGTTCAATGGCCCAGCACCTTGCCGGCGGCCAACTTCTcgcg                                                                          >  8:335115/1‑149 (MQ=255)
                                                                         ccGCGGTCGCTTTCATGTAGCTTTTACAGTTGCGCCGGCAGCGAGCTGTTGAGCGAAGGGTTTGCGTTGATCAGAGGGGCCAGCTCACGTGCGGCGGCGAGGCGTGGGTTCATTGGACGTGGACTTTGCGGGGGGAGAACTtctggggg                                                                         <  6:43951/149‑6 (MQ=255)
                                                                         ccGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTgg                                                                                                                   <  1:105165/107‑1 (MQ=255)
                                                                         ccGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTgg                                                                                                                   >  2:105165/1‑107 (MQ=255)
                                                                         ccGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCgcgg                                                                         <  3:366689/149‑1 (MQ=255)
                                                                         ccGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGACCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCgcgg                                                                         >  3:118543/1‑149 (MQ=255)
                                                                             ggTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGggcggcgg                                                                                                                             <  3:315764/93‑1 (MQ=255)
                                                                             ggTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGggcggcgg                                                                                                                             >  4:315764/1‑93 (MQ=255)
                                                                                  gcgTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCTCAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGtt                                                                 <  4:237387/149‑1 (MQ=255)
                                                                                   cgTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTg                                                                <  5:380124/149‑1 (MQ=255)
                                                                                    gTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTg                                                                >  6:360290/1‑148 (MQ=255)
                                                                                        aTCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTg                                                                   >  3:368203/1‑141 (MQ=255)
                                                                                        aTCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTg                                                                   <  4:368203/141‑1 (MQ=255)
                                                                                        aTCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAg                                                           >  3:231052/1‑149 (MQ=255)
                                                                                            aCCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTc                                                        >  4:202197/1‑148 (MQ=255)
                                                                                            aCCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCa                                                       >  3:287948/1‑149 (MQ=255)
                                                                                                 gTCCAGTTGCGCCGGCAGCGAGCTGTTCATCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTg                                                    <  7:17189/147‑1 (MQ=255)
                                                                                                  tCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGaa                                                  <  2:171713/148‑1 (MQ=255)
                                                                                                    cAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGc                                                <  4:216023/148‑1 (MQ=255)
                                                                                                         gcgcCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAAt                                           <  5:185696/148‑1 (MQ=255)
                                                                                                             cGGCAGCGAGTTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCtt                                      <  2:76935/149‑1 (MQ=255)
                                                                                                             cGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCtt                                      >  8:17245/1‑149 (MQ=255)
                                                                                                               gCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTgc                                    >  5:149549/1‑149 (MQ=255)
                                                                                                                    gAGCTGTTCAGCGTAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTTGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCa                               <  7:335115/149‑1 (MQ=255)
                                                                                                                     aGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCcg                                                                                           >  2:72443/1‑88 (MQ=255)
                                                                                                                     aGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCcg                                                                                           <  1:72443/88‑1 (MQ=255)
                                                                                                                      gCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAgg                             <  4:231052/149‑1 (MQ=255)
                                                                                                                         gTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGcggcgg                                                                      <  3:277783/105‑1 (MQ=255)
                                                                                                                         gTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGcggcgg                                                                      >  4:277783/1‑105 (MQ=255)
                                                                                                                         gTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAg                          >  2:7177/1‑149 (MQ=255)
                                                                                                                            cAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGtt                        <  1:7177/148‑1 (MQ=255)
                                                                                                                               cGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAg                    >  5:4621/1‑149 (MQ=255)
                                                                                                                                   gCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCgg                <  2:63087/149‑1 (MQ=255)
                                                                                                                                   gCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTACAGGCgg                <  4:118543/149‑1 (MQ=255)
                                                                                                                                    ccTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCgg                >  3:191903/1‑148 (MQ=255)
                                                                                                                                    ccTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCgg                <  4:191903/148‑1 (MQ=255)
                                                                                                                                     ctggCCGCTGTGCTCAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAAt                                           <  5:224971/117‑1 (MQ=255)
                                                                                                                                     ctggCCGCTGTGCTCAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAAt                                           >  6:224971/4‑120 (MQ=255)
                                                                                                                                      ttGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGCACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTc             >  4:139036/1‑149 (MQ=255)
                                                                                                                                        gCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTc             <  3:202197/147‑1 (MQ=255)
                                                                                                                                        gCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa            <  6:325582/148‑1 (MQ=255)
                                                                                                                                           gCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCt        >  8:25028/1‑149 (MQ=255)
                                                                                                                                           gCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTCGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCt        >  2:300626/1‑149 (MQ=255)
                                                                                                                                            cTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCt        >  5:316186/1‑148 (MQ=255)
                                                                                                                                               aGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAgg    >  6:172662/1‑149 (MQ=255)
                                                                                                                                                 cacaGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGGGc  <  4:287948/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGGGC  >  NC_002947/77358‑77650

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: