Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,803,922 G→T 6.4% L245L (CTG→CTT)  rbsC → D‑ribose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,803,9220GT 6.4% 200.6 / 3.4 62L245L (CTG→CTT) rbsCD‑ribose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base T (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TGTTCGGCCGTTACCTGATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCACCTTCTTTGGCGTACTGATCATTTCCGTATTGGCCGCCGGG  >  NC_002947/2803777‑2804069
                                                                                                                                                  |                                                                                                                                                   
tgtTCGGCCGTTACCTGATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGaa                                                                                                                                                  <  6:80671/149‑1 (MQ=255)
     ggCCGTTACCTGATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAgccgc                                                                                                                                             >  6:251235/1‑149 (MQ=255)
      gCCGTTACCTGATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGGATATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTCTTCCAGATCTCGCGGCTGGAAgccgc                                                                                                                                             >  1:243899/1‑149 (MQ=255)
      gCCGTTACCTGATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCTTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAgccgcc                                                                                                                                            >  7:221784/1‑149 (MQ=255)
       ccGTTACCTGATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGccgccg                                                                                                                                           >  7:67603/1‑149 (MQ=255)
        cGTTACCTGATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTCATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGa                                                                                                                                          >  3:42204/1‑149 (MQ=255)
               tGATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAAt                                                                                                                                   >  7:106898/1‑149 (MQ=255)
                gATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTg                                                                                                                                                     <  8:292190/130‑1 (MQ=255)
                gATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTg                                                                                                                                                     >  7:292190/1‑130 (MQ=255)
                 aTCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAAt                                                                                                                                   <  2:236559/147‑1 (MQ=255)
                   cGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCg                                                                                                                               >  1:127784/1‑149 (MQ=255)
                    gttATGGGCACCAACGAAGAGGCCGGGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCGTGATGGGGGTGCTCGCGGGCCTGGCCGAACTGTTCCATATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCgg                                                                                                                              <  2:243899/147‑1 (MQ=255)
                    ggaaTCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TAGCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCgg                                                                                                                              <  3:85653/146‑1 (MQ=255)
                    ggTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCgg                                                                                                                              <  4:125230/149‑1 (MQ=255)
                    agTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGAGCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCg                                                                                                                               <  5:176955/147‑1 (MQ=255)
                        tCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGctct                                                                                                                          <  3:95586/149‑1 (MQ=255)
                                   gAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCt                                                                                                               >  3:45613/1‑149 (MQ=255)
                                     agagGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGc                                                                                                             <  5:85569/149‑1 (MQ=255)
                                        ggCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGACCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGCCCCCAATGCCGGCTCTGGCCTGGAGCTGCAgg                                                                                                          >  7:283754/1‑149 (MQ=255)
                                          ccGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTc                                                                                                        <  2:103541/149‑1 (MQ=255)
                                                gcCTGGCCGG‑TTTCGATCCGCGCCTTTAAAATGTGCTGGTGTTCGCCGTGAGGGGCCTGCGCGCCGGCCTGGCCGCAGTGTTCCTTATCTGGGGGGTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTgcc                                                                                                  <  8:106898/149‑1 (MQ=255)
                                                     gCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCcgt                                                                                             >  1:223480/1‑149 (MQ=255)
                                                      ccGG‑GATCGATCCGCGCCCTTCCCACGTGAGGGGGTTCGCCCTGATGGGCCGGCTCGCCGGCCTGGCCGCACTGTGCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCcgtc                                                                                            <  5:22683/149‑1 (MQ=255)
                                                               gATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGccgccgatcac                                                                                         >  3:146452/1‑139 (MQ=255)
                                                                aTCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGTTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATcggc                                                                                   <  2:108326/149‑1 (MQ=255)
                                                                 tCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATcggcg                                                                                  <  5:107671/149‑1 (MQ=255)
                                                                  ccGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATcggcgg                                                                                 >  7:223074/1‑149 (MQ=255)
                                                                  ccGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATcggcgg                                                                                 <  4:45613/149‑1 (MQ=255)
                                                                   cgcgCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATcggcgg                                                                                 >  8:113010/1‑148 (MQ=255)
                                                                         cTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGc                                                                          <  6:150167/149‑1 (MQ=255)
                                                                          ttACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTTCCAgccgcga                                                                                                                                           <  8:127646/83‑3 (MQ=255)
                                                                          ttACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTTCCAgccgcga                                                                                                                                           >  7:123328/1‑81 (MQ=255)
                                                                          ttACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTTCCAgccgcga                                                                                                                                           >  7:127646/1‑81 (MQ=255)
                                                                          ttACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTTCCAgccgcga                                                                                                                                           <  8:123328/83‑3 (MQ=255)
                                                                            aCAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCt                                                                        <  4:42204/148‑1 (MQ=255)
                                                                            aCAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGATCATTGCCGCCGTCGTGATCGGCGGCACCAGCCt                                                                        >  6:296005/1‑148 (MQ=255)
                                                                              aaaGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGAt                                                                     >  6:265438/1‑149 (MQ=255)
                                                                              aaaGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGAt                                                                     >  2:262871/1‑149 (MQ=255)
                                                                               aaGTGGTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGACGTGATCGGCGGCACCAGCCTGAtg                                                                    >  3:168984/1‑149 (MQ=255)
                                                                                aGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCcgtc                                                                                            >  6:138778/1‑124 (MQ=255)
                                                                                aGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCcgtc                                                                                            <  5:138778/124‑1 (MQ=255)
                                                                                     tGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGAtgggtggg                                                              >  5:262132/1‑149 (MQ=255)
                                                                                     gggTGTGCCCCCTGATGGGCCTGCGCGCCGGCCGGGCCGCACTGTGCCAGATCTCGCGGCTGGAAGCCGCCGACCCCTATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGAtgggtgg                                                               <  6:18591/147‑1 (MQ=255)
                                                                                       gtgtTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCTATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGAtgggtgg                                                               >  5:18591/1‑146 (MQ=255)
                                                                                        tgtTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGt                                                           <  8:67603/149‑1 (MQ=255)
                                                                                        tgtTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGt                                                           <  5:251235/149‑1 (MQ=255)
                                                                                          ttCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTgg                                                         >  7:86966/1‑149 (MQ=255)
                                                                                          ttCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGATCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTgg                                                         <  5:296005/149‑1 (MQ=255)
                                                                                              cccTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCg                                                     >  4:7210/1‑149 (MQ=255)
                                                                                               ccTGATGGGCCTGCTCGCCGGCCGGGCCGCACTGTTCCAGATGTCGCGGCTGGATGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCTCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCgg                                                    <  6:285400/149‑1 (MQ=255)
                                                                                                cTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCggc                                                   <  2:223480/149‑2 (MQ=255)
                                                                                                      ggCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAg                                             <  6:262132/149‑1 (MQ=255)
                                                                                                      ggCCTGCTCGCCGCCCGGCCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGtcatca                                              <  2:238725/148‑1 (MQ=255)
                                                                                                        ccTGCTCGCCGGCCTGGCCGCACTGTTCCAGAACTCGCGGCGGGACGCCGACGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCGCGGTCATCAGCa                                           <  8:177203/149‑1 (MQ=255)
                                                                                                         cTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGACGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCAc                                          >  2:3911/1‑149 (MQ=255)
                                                                                                          tGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGAGCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCAGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCAcc                                         >  3:122425/1‑149 (MQ=255)
                                                                                                                 cGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCACCTTCTTTg                                  <  1:262871/149‑1 (MQ=255)
                                                                                                                       ggCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGACGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCACCTTCTTTGGCGTAc                            <  1:3911/149‑1 (MQ=255)
                                                                                                                            cACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCACCTTCTTTGGCGTACTGAt                        >  3:179591/1‑148 (MQ=255)
                                                                                                                                     aGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCACCTTCTTTGGCGTACTGATCATTTCCGTa              >  5:254040/1‑149 (MQ=255)
                                                                                                                                      gATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCACCTTCTTTGGCGTACTGATCATTTCCGTAt             >  1:136461/1‑149 (MQ=255)
                                                                                                                                         ctcGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCACCTTCTTTGGCGTACTGATCATTTCCGTATTgg          >  3:249121/1‑149 (MQ=255)
                                                                                                                                                 tGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCACCTTCTTTGGCGTACTGATCATTTCCGTATTGGCCGCCggg  >  4:255658/1‑149 (MQ=255)
                                                                                                                                                 tGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCACCTTCTTTGGCGTACTGATCATTTCCGTATTGGCCGCCggg  >  4:260039/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGTTCGGCCGTTACCTGATCGGTATCGGCACCAACGAAGAGGCCGTGCGCCTGGCCGG‑TATCGATCCGCGCCCTTACAAAGTGCTGGTGTTCGCCCTGATGGGCCTGCTCGCCGGCCTGGCCGCACTGTTCCAGATCTCGCGGCTGGAAGCCGCCGACCCCAATGCCGGCTCTGGCCTGGAGCTGCAGGTCATTGCCGCCGTCGTGATCGGCGGCACCAGCCTGATGGGTGGGCGTGGCTCGGTCATCAGCACCTTCTTTGGCGTACTGATCATTTCCGTATTGGCCGCCGGG  >  NC_002947/2803777‑2804069

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