Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 853,762 G→A 9.6% R62H (CGT→CAT) ‡ prmC → release factor‑(glutamine‑N5) methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947853,7620GA 9.6% 169.2 / 5.3 62R62H (CGT→CAT) ‡prmCrelease factor‑(glutamine‑N5) methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base A (4/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

GCCTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTCGC  >  NC_002947/853616‑853902
                                                                                                                                                  |                                                                                                                                                     
gCCTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTg                                                                                                                                                     >  7:95357/1‑149 (MQ=255)
  cTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTggc                                                                                                                                                   >  6:313149/1‑149 (MQ=255)
       tCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCACGgcggcgc                                                                                                                                              <  7:281513/149‑2 (MQ=255)
         gCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGc                                                                                                                                            >  1:123974/1‑149 (MQ=255)
           ccACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGc                                                                                                                                            >  1:159601/1‑147 (MQ=255)
           ccACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCACGgcggcgcttca                                                                                                                                          <  6:156254/149‑6 (MQ=255)
            cACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGa                                                                                                                                              <  1:293429/144‑1 (MQ=255)
            cACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGa                                                                                                                                              >  2:293429/1‑144 (MQ=255)
                    gcCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGAGGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCtat                                                                                                                                 <  3:204373/149‑1 (MQ=255)
                     ccTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCtata                                                                                                                                >  1:258385/1‑149 (MQ=255)
                      cTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTatat                                                                                                                               >  8:20715/1‑149 (MQ=255)
                        ggATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACGTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATcc                                                                                                                             <  1:119666/149‑1 (MQ=255)
                        ggATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATcc                                                                                                                             >  4:120927/1‑149 (MQ=255)
                             ccGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTgg                                                                                                                         <  5:147797/148‑1 (MQ=255)
                              cGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGc                                                                                                                       >  3:246265/1‑149 (MQ=255)
                                     ctgctgGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAg                                                                                                                  <  8:95357/147‑1 (MQ=255)
                                          ggctgCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCt                                                                                                           <  8:218621/149‑1 (MQ=255)
                                          ggctgCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCt                                                                                                           >  5:8761/1‑149 (MQ=255)
                                               cGGCCATCGGTAAATCGCGCTGCTTCCTGCATACCTGGCCCGCGCGTATGGTCAGCAGCGAAGATGCCGAGTCCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGaa                                                                                                       <  7:258893/148‑1 (MQ=255)
                                                  ccATCGGTAAATGCCGCAGCTACCTTCAAACCTGGCCCTAGCGTATTGTCCGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATc                                                                                                   <  4:4499/149‑1 (MQ=255)
                                                    aTCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGa                                                                                                 >  7:171956/1‑149 (MQ=255)
                                                     tCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGAc                                                                                                <  6:277227/149‑1 (MQ=255)
                                                         tAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTgg                                                                                            <  2:258385/149‑1 (MQ=255)
                                                             tCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGt                                                                                    >  4:153119/1‑149 (MQ=255)
                                                                 gcAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTcgcg                                                                                >  1:44963/1‑149 (MQ=255)
                                                                  cAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCgcgc                                                                               <  8:171956/149‑1 (MQ=255)
                                                                    gCTACCTGAATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGcgcc‑‑‑‑‑g                                                                        <  2:159601/149‑1 (MQ=255)
                                                                            cATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATAccc                                                                 <  3:120927/148‑1 (MQ=255)
                                                                                 cTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATcc                                                           <  4:246265/149‑1 (MQ=255)
                                                                                 cTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGCCTTCTGGAAGATCGACCTCGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATcc                                                           >  5:298722/1‑149 (MQ=255)
                                                                                      ccGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATAcc                                                                  >  6:50712/1‑137 (MQ=255)
                                                                                      ccGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATAcc                                                                  <  5:50712/137‑1 (MQ=255)
                                                                                      ccGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCgg                                                      >  6:98222/1‑149 (MQ=255)
                                                                                        gAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGcc                                                    >  2:149895/1‑149 (MQ=255)
                                                                                             tATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATa                                               <  3:223597/149‑1 (MQ=255)
                                                                                               ttGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATAcc                                             <  6:199488/149‑1 (MQ=255)
                                                                                                  tCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAg                                          >  5:313736/1‑149 (MQ=255)
                                                                                                     gcagcGAAGATGCCGAGACCTTTGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATAACCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGcgcc‑‑‑‑‑gc                                                                       <  1:266881/117‑1 (MQ=255)
                                                                                                     gcagcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATAACCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGcgcc‑‑‑‑‑gc                                                                       >  2:266881/1‑117 (MQ=255)
                                                                                                      cagcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAgctgc                                      >  1:170747/1‑149 (MQ=255)
                                                                                                       agcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGAGGCCTATATCCTTGGTCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑CGTCGCGCC‑‑‑‑‑GCATACACTGATCCCGCGGCCGGATACCGAgctgct                                     >  4:117923/1‑149 (MQ=255)
                                                                                                          gAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGt                                  <  5:313149/149‑1 (MQ=255)
                                                                                                           aaGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTc                                 <  2:44963/149‑1 (MQ=255)
                                                                                                            aGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCg                                <  6:313736/149‑1 (MQ=255)
                                                                                                                       ccTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGGAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCtt                     <  7:93087/149‑1 (MQ=255)
                                                                                                                        cTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCgg                                                  >  3:127093/1‑119 (MQ=255)
                                                                                                                        cTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCgg                                                  <  4:127093/119‑1 (MQ=255)
                                                                                                                        cTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATAc                                              >  5:94951/1‑123 (MQ=255)
                                                                                                                        cTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATAc                                              <  6:94951/123‑1 (MQ=255)
                                                                                                                        cTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCtt                     >  7:77418/1‑148 (MQ=255)
                                                                                                                           tGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGaa                  <  1:149895/148‑1 (MQ=255)
                                                                                                                           tGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAAc                 >  2:297804/1‑149 (MQ=255)
                                                                                                                              cGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGc              <  5:98222/149‑1 (MQ=255)
                                                                                                                                   aCCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCGATATCCTTGGGCTGCAGGGCGTCTGGAAGCTCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAAc           >  5:60148/1‑147 (MQ=255)
                                                                                                                                    ccTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCg        >  6:262111/1‑149 (MQ=255)
                                                                                                                                      tggctcgccgcCACGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCt      >  1:56526/6‑149 (MQ=255)
                                                                                                                                      tGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAAAAGGGGGGCCCCGGGAGGGATACACAGATCCCGCGGCCGGATACCGAGCTGATGGTCGAGACGGCGCTTGTAcag               >  1:265331/1‑147 (MQ=255)
                                                                                                                                       gAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCtt     >  1:270225/1‑149 (MQ=255)
                                                                                                                                        gctcgccgcCACGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTc    >  8:79474/4‑149 (MQ=255)
                                                                                                                                        aaGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTc    <  1:297804/149‑1 (MQ=255)
                                                                                                                                         ctcgccgcCACGGCGGCGAGCCAGTGGCCTTTATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTcg   >  8:192298/3‑149 (MQ=255)
                                                                                                                                          tcgccgcCACGGCGGCGAGCCAGTGGCCTATATCCATGGCCTGCACGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATCACCTGATCCCGCGGCCGGATACCGAGCTCCTGGTCGAGACGGCGCTTCAACTGCAACCCGCTTcgc  >  6:63765/2‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCCTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGA‑‑‑‑GGTCGCGCC‑‑‑‑‑GCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTCGC  >  NC_002947/853616‑853902

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