Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,059,508:1 +G 5.2% coding (1326/1581 nt) PP_4458 ← opine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,059,5081.G 5.2% 320.5 / 9.0 77coding (1326/1581 nt)PP_4458opine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (36/37);  new base G (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

CTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGTGATCCAC  >  NC_002947/5059361‑5059654
                                                                                                                                                     |                                                                                                                                                  
cTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGACTCGACCCTGGCTTTCCGGATGGCCTCATCAGCTAATGAGCAATGTACGCAATTc                                                                                                                                                     >  8:338379/1‑148 (MQ=255)
   tttCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACa                                                                                                                                                <  4:109057/149‑1 (MQ=255)
        gACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTg                                                                                                                                             >  4:92986/1‑147 (MQ=255)
          ccAAACAAAGGGATGGCGCATTCGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGTATCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCGAATTTCGGACATTGctca                                                                                                                                         <  1:118171/149‑4 (MQ=255)
          ccAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGc                                                                                                                                         >  4:237868/1‑149 (MQ=255)
          ccAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGACTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGc                                                                                                                                         >  3:2381/1‑149 (MQ=255)
            aaacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCgg                                                                                                                                       <  7:257125/149‑1 (MQ=255)
                         gcgcATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGc                                                                                                                          <  6:275942/149‑1 (MQ=255)
                         gcgcATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGACTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATgg                                                                                                                           <  4:2381/148‑1 (MQ=255)
                          cgcATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCa                                                                                                                         <  2:100648/149‑1 (MQ=255)
                          cgcATACGTCGTCGTGAAGCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCa                                                                                                                         <  4:40033/149‑1 (MQ=255)
                           gcATACGTCGTCCTGAATCTCAACCTGCGCCTTTTTACACAAATCGAGTTGAGTATCACAGATCAGCCACGACCCTCTCTTTCCGCATGGCATCAAAAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCTCCGATGGCa                                                                                                                         >  1:107792/1‑149 (MQ=255)
                                  tcgtcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGa                                                                                                                 <  3:159669/149‑1 (MQ=255)
                                   cgtcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGGCTGCCGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGaa                                                                                                                <  3:237098/149‑1 (MQ=255)
                                    gtcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAAt                                                                                                               >  6:67654/1‑149 (MQ=255)
                                    gtcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAAt                                                                                                               >  1:234837/1‑149 (MQ=255)
                                    gtcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCAGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGaa                                                                                                                <  3:328969/148‑1 (MQ=255)
                                    gtcgtGAATCTGAACCTCCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAAt                                                                                                               <  3:274387/149‑1 (MQ=255)
                                     tcgtGTATCTTTACCTGCGCGTTTTTCCTCAAATCTAGTTGTGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGATTTTCGTCTTTGCAGCTGGTGCGGCTATGGCAGAAGCCGAATc                                                                                                              <  2:107792/149‑1 (MQ=255)
                                               gAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAActct                                                                                                    >  3:211682/1‑149 (MQ=255)
                                               gAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAActct                                                                                                    >  7:298709/1‑149 (MQ=255)
                                               gAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAActct                                                                                                    >  2:382109/1‑149 (MQ=255)
                                                 aCCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGt                                                                                                  <  3:237868/149‑1 (MQ=255)
                                                   cTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCa                                                                                                <  7:287682/149‑1 (MQ=255)
                                                      cgcCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGc                                                                                              <  4:211682/148‑1 (MQ=255)
                                                       gcCCTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCa                                                                                            <  3:12149/149‑1 (MQ=255)
                                                         cTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCag                                                                                           >  2:386261/1‑148 (MQ=255)
                                                            tttCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCagtagt                                                                                       >  7:8898/1‑149 (MQ=255)
                                                              tCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCag                                                                                           <  5:104643/143‑1 (MQ=255)
                                                              tCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCag                                                                                           >  6:104643/1‑143 (MQ=255)
                                                                 acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgcg                                                                                  >  3:114673/1‑149 (MQ=255)
                                                                 acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGAATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTATTCGgcg                                                                                  >  1:149074/1‑149 (MQ=255)
                                                                 acaAATAGAGTTGAGTATGCGGATCAGCCTCGACCCTGGCTTTGCGGATGGCATCATCAGATACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgcg                                                                                  >  2:253978/1‑149 (MQ=255)
                                                                  caAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgcgc                                                                                 <  2:234837/149‑1 (MQ=255)
                                                                           ttGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATc                                                                        <  3:185354/149‑1 (MQ=255)
                                                                              aGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                                     <  3:198664/149‑1 (MQ=255)
                                                                              aGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                                     <  8:298709/149‑1 (MQ=255)
                                                                               gTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGagc                                                                    <  3:92986/149‑1 (MQ=255)
                                                                               gTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGagc                                                                    <  1:382109/149‑1 (MQ=255)
                                                                                tCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                                     >  2:250092/1‑147 (MQ=255)
                                                                                  tGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGagcag                                                                  >  6:317070/1‑148 (MQ=255)
                                                                                      ggATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCAcc                                                              <  1:250092/149‑1 (MQ=255)
                                                                                       gATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCg                                                             >  5:187883/1‑149 (MQ=255)
                                                                                         tCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGt                                                            >  5:257639/1‑148 (MQ=255)
                                                                                            gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTaaaa                                                        >  3:237050/1‑149 (MQ=255)
                                                                                            gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTaaaa                                                        <  7:338379/149‑1 (MQ=255)
                                                                                            gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTaaaa                                                        >  7:282668/1‑149 (MQ=255)
                                                                                              cTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAaca                                                      >  1:307989/1‑149 (MQ=255)
                                                                                                     cTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATcgc                                               >  4:257362/1‑149 (MQ=255)
                                                                                                      tGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCgg                                                                             >  7:365643/1‑118 (MQ=255)
                                                                                                      tGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCgg                                                                             <  8:365643/118‑1 (MQ=255)
                                                                                                       ggCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCgcgc                                             >  4:347329/1‑149 (MQ=255)
                                                                                                       ggCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCgcgc                                             >  4:353010/1‑149 (MQ=255)
                                                                                                        gCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGct                                            <  2:307989/149‑1 (MQ=255)
                                                                                                          tttCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGctct                                          >  2:9621/1‑149 (MQ=255)
                                                                                                                gATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCtt                                    >  1:211197/1‑149 (MQ=255)
                                                                                                                gATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGGAAAACACGATCGCGCTCTGGTCtt                                    >  6:319939/1‑149 (MQ=255)
                                                                                                                      atcatcAGCTACTGAGCAATGTCCGAATTTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGc                              >  1:16702/1‑149 (MQ=255)
                                                                                                                      atcatcAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCgg                                                                             <  1:20274/102‑1 (MQ=255)
                                                                                                                      atcatcAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCgg                                                                             >  2:20274/1‑102 (MQ=255)
                                                                                                                         atcaGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCtt                            <  2:149074/148‑1 (MQ=255)
                                                                                                                            aGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTg                         <  5:317070/148‑1 (MQ=255)
                                                                                                                            aGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTg                         <  5:67654/148‑1 (MQ=255)
                                                                                                                             gCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGc                       <  4:237050/149‑1 (MQ=255)
                                                                                                                             gCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGc                       <  4:114673/149‑1 (MQ=255)
                                                                                                                               tACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCt                     >  1:369813/1‑149 (MQ=255)
                                                                                                                                 cTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTg                    <  2:369813/148‑1 (MQ=255)
                                                                                                                                 cTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTg                    <  6:257639/148‑1 (MQ=255)
                                                                                                                                  gctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCa                                                                                            >  2:113652/4‑75 (MQ=25)
                                                                                                                                  gctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCa                                                                                            <  1:113652/72‑1 (MQ=25)
                                                                                                                                   gAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGt                  >  7:287121/1‑148 (MQ=255)
                                                                                                                                    tgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGt                  >  7:119135/2‑147 (MQ=255)
                                                                                                                                     gCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGGAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTaa                <  5:319939/148‑1 (MQ=255)
                                                                                                                                         tGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTcgc           <  3:257362/149‑1 (MQ=255)
                                                                                                                                              gAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGAAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTcgc           >  1:184119/1‑144 (MQ=255)
                                                                                                                                              gAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCAGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGTGAt      >  3:89565/1‑149 (MQ=255)
                                                                                                                                                  ttCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGTGATCCAc  >  8:73908/1‑149 (MQ=255)
                                                                                                                                                  ttCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGACAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGTGATCCAc  <  8:8898/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGTGATCCAC  >  NC_002947/5059361‑5059654

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: