Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 559,579 G→T 8.7% G43V (GGT→GTT)  rplO → 50S ribosomal protein L15

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947559,5790GT 8.7% 170.9 / 4.4 57G43V (GGT→GTT) rplO50S ribosomal protein L15
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base T (2/3);  total (28/29)
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.47e-01

GGCTAATCGATGAAACTCAATGATCTGAGTCCAGCGCCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCCAAAGTGGAAGGCGACGTGATCTCCGTG  >  NC_002947/559443‑559727
                                                                                                                                         |                                                                                                                                                      
ggCTAATCGATGAAACTCAATGATCTGAGTCCAGCGCCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAACGTCAGACTTCCCGTTCGGGTGGTTCGATCg                                                                                                                                             >  2:144663/1‑148 (MQ=255)
     aTCGATGAAACTCAATGATCTGAGTCCAGCGCCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTcc                                                                                                                                        <  4:249072/147‑1 (MQ=255)
      tCGATGAAACTCAATGATCTGAGTCCAGCGCCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTgg                                                                                                                                     >  7:228908/1‑149 (MQ=255)
                        cTGAGTCCAGCGCCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACa                                                                                                                  >  8:57625/1‑149 (MQ=255)
                               cAGCGCCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTg                                                                                                           >  5:277973/1‑149 (MQ=255)
                                 gcgcCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCa                                                                                                         >  1:235954/1‑149 (MQ=255)
                                 gcgcCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCa                                                                                                         <  3:138575/149‑1 (MQ=255)
                                  cgcCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCAc                                                                                                        <  6:277973/149‑1 (MQ=255)
                                   gcCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGACGCTGcccc                                                                                                       >  4:239009/1‑149 (MQ=255)
                                         ttCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGcc                                                                                                 >  7:130899/1‑149 (MQ=255)
                                           cccGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGTGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCt                                                                                                <  7:127750/148‑1 (MQ=255)
                                           cccGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCt                                                                                                <  6:119472/148‑1 (MQ=255)
                                                  cgAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAg                                                                                         <  8:228908/148‑1 (MQ=255)
                                                     gaAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGATCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCg                                                                                     <  5:318938/149‑1 (MQ=255)
                                                       aGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGc                                                                                   <  8:130899/149‑1 (MQ=255)
                                                          aTCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTc                                                                                <  2:235954/149‑1 (MQ=255)
                                                              tCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGcc                                                                                                               <  7:227247/114‑1 (MQ=255)
                                                              tCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGcc                                                                                                               >  8:227247/1‑114 (MQ=255)
                                                              tCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGttt                                                                            >  3:76348/1‑149 (MQ=255)
                                                                   gTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGGTCGGCTTCGTTTCCCt                                                                        >  2:220967/1‑148 (MQ=255)
                                                                    tCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGTTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTg                                                                       >  5:325034/1‑148 (MQ=255)
                                                                    tCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGTGGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTg                                                                       >  5:231375/1‑148 (MQ=255)
                                                                      gTGGTATCGGTAGCGGTCTGGTTAAGACCGGCGGCCGTGGCCAAAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGaaa                                                                    <  6:282665/149‑1 (MQ=255)
                                                                        ggTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGc                                                                  >  5:49692/1‑149 (MQ=255)
                                                                              ggTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATgg                                                             <  2:100642/148‑1 (MQ=255)
                                                                              ggTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGa                                                            >  1:110534/1‑149 (MQ=255)
                                                                                  gCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAAc                                                                                                                   <  7:280643/90‑1 (MQ=255)
                                                                                  gCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAAc                                                                                                                   >  8:280643/1‑90 (MQ=255)
                                                                                  gCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGCTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACcgc                                                        <  4:86892/149‑1 (MQ=255)
                                                                                   cGGTCTGGGTAAGACCGGCGGCCGTGGCCACAACAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTACGTTTCCCTGAAAGCCATGGACcgc                                                        >  1:201045/1‑149 (MQ=255)
                                                                                    ggTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGTTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCgcgc                                                      <  6:325034/149‑1 (MQ=255)
                                                                                    ggTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGTTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCgcgc                                                      <  6:231375/149‑1 (MQ=255)
                                                                                    ggTCTGGGTAAGACCAGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGATCGATCACTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCgcgc                                                      <  8:199408/149‑1 (MQ=255)
                                                                                     gTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCt                                                     >  5:106492/1‑149 (MQ=255)
                                                                                     gTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCt                                                     <  2:201045/149‑1 (MQ=255)
                                                                                          ggTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGaa                                                  <  6:130550/147‑1 (MQ=255)
                                                                                           gTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCTGTGGCTACGAAGGCGGTCAACAGCCGGTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCGTGAAAGCCATGGACCGCGCTGAAGt                                                >  1:92312/1‑149 (MQ=255)
                                                                                            tAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGc                                              >  5:329411/1‑149 (MQ=255)
                                                                                                aCCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCt                                          >  5:156389/1‑149 (MQ=255)
                                                                                                  cggcggCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGTTGGTTCGATCGCTCCTTGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGAcc                                                          <  1:323799/131‑1 (MQ=255)
                                                                                                  cggcggCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGAcc                                                          >  2:323799/1‑131 (MQ=255)
                                                                                                  cggcggCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGt                                        <  7:57625/149‑1 (MQ=255)
                                                                                                    gcggcCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTcc                                      >  8:236404/1‑149 (MQ=255)
                                                                                                        ccGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACAGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGc                                  >  7:20398/1‑149 (MQ=255)
                                                                                                          gTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGTTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTg                                <  3:286121/149‑1 (MQ=255)
                                                                                                             gCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGc                              <  6:49692/148‑1 (MQ=255)
                                                                                                              ccACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCCa                            <  4:76348/149‑1 (MQ=255)
                                                                                                               cacaAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGcc                             >  1:245623/1‑148 (MQ=255)
                                                                                                                   aaGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCCAAAGTg                       <  2:92312/149‑1 (MQ=255)
                                                                                                                        tCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCt                                          >  2:47649/1‑126 (MQ=255)
                                                                                                                        tCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCt                                          <  1:47649/126‑1 (MQ=255)
                                                                                                                        tCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCCAAAGTGGAAg                   >  8:89024/1‑149 (MQ=255)
                                                                                                                        tCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCCAAAGTGGAAg                   <  6:106492/149‑1 (MQ=255)
                                                                                                                               tCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCCAAAGTGGAAGGCGACGt            >  7:305183/1‑149 (MQ=255)
                                                                                                                                cccGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAACTGGCCAAAGTGGAAGGCGACg             <  2:110534/147‑1 (MQ=255)
                                                                                                                                 ccGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCGAAAGTGGAAGGCGACGTGa          >  3:182510/1‑149 (MQ=255)
                                                                                                                                  cGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCCAAAGTGGAAGGCGACGTGAt         <  1:144663/149‑1 (MQ=255)
                                                                                                                                         gtggtTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCCAAAGTGGAAGGCGACGTGATCTCCGTg  <  7:236404/149‑1 (MQ=255)
                                                                                                                                         gtggtTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCCAAAGTGGAAGGCGACGTGATCTCCGTg  <  6:156389/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                      
GGCTAATCGATGAAACTCAATGATCTGAGTCCAGCGCCGGGTTCCCGTCGCGAGAAGCATCGTCCAGGTCGTGGTATCGGTAGCGGTCTGGGTAAGACCGGCGGCCGTGGCCACAAAGGTCAGACTTCCCGTTCGGGTGGTTCGATCGCTCCTGGCTTCGAAGGCGGTCAACAGCCGCTGCACCGTCGCCTGCCGAAGTTCGGCTTCGTTTCCCTGAAAGCCATGGACCGCGCTGAAGTGCGTCTGTCCGAGCTGGCCAAAGTGGAAGGCGACGTGATCTCCGTG  >  NC_002947/559443‑559727

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: