Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,059,500 Δ1 bp 8.1% coding (1334/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,5000C. 8.1% 227.2 / 14.3 74coding (1334/1581 nt)PP_4458opine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (35/32);  new base . (3/3);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

GCCCAGACCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAG  >  NC_002947/5059353‑5059641
                                                                                                                                                   |                                                                                                                                              
gCCCAGACCTGTTTCAGACCAAACAAAGGGATGGCGCAAACCACGTCGGGAATCTGAACCTGCGCGTTTTTCCACAAATCGAGTTGAGTATGCTGATCAGCCTCGCCCCTGGCTTTCCGGATGGCATCAGCGGGTACTGACCAATGTcc                                                                                                                                               >  6:115593/1‑149 (MQ=255)
 cccAGACCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTcc                                                                                                                                               >  2:160200/1‑148 (MQ=255)
 cccAGACCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTcc                                                                                                                                               >  7:198440/1‑148 (MQ=255)
   cAGACCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGa                                                                                                                                             >  1:56314/1‑148 (MQ=255)
      aCCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAAtt                                                                                                                                         >  4:127678/1‑149 (MQ=255)
       ccttgttCAGTCCAAACAAAGGGATGGCGCATACGTCGGGGTGAATCTGAACCTGCGCCTTTTTCGACTACTCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTCCTGAGCAATGTCCGAAATTc                                                                                                                                        <  1:118009/144‑1 (MQ=255)
         tGTTTCAGACCAAACAAAGGGAGGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGAAGGCATCATCAGCTAGTGAGCAATGTCCGAAATTCGa                                                                                                                                     >  7:227561/1‑149 (MQ=255)
           tttCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCTGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGAc                                                                                                                                    <  8:176636/148‑1 (MQ=255)
           cctCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGa                                                                                                                                     <  8:227561/145‑1 (MQ=255)
                gACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTg                                                                                                                                >  5:166375/1‑147 (MQ=255)
                 aCCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCGAATTTCGGACATTGctc                                                                                                                             <  3:94185/149‑3 (MQ=255)
                   caaacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCg                                                                                                                           >  8:282530/1‑149 (MQ=255)
                    aaacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCgg                                                                                                                          <  2:56314/149‑1 (MQ=255)
                     aacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCggg                                                                                                                         >  7:303033/1‑149 (MQ=255)
                       caaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTg                                                                                                                       <  6:135459/149‑1 (MQ=255)
                        aaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTgc                                                                                                                      >  7:29259/1‑149 (MQ=255)
                        aaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTgc                                                                                                                      >  1:348426/1‑149 (MQ=255)
                                 gcgcATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATgg                                                                                                              <  3:344529/148‑1 (MQ=255)
                                 gcgcATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATgg                                                                                                              <  5:293432/148‑1 (MQ=255)
                                             tcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATc                                                                                                 <  8:198440/149‑1 (MQ=255)
                                                   aTCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAATTTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTaaaa                                                                                           <  5:160591/149‑1 (MQ=255)
                                                         aCCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGt                                                                                     <  6:166375/149‑1 (MQ=255)
                                                           cTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCa                                                                                   <  8:29259/149‑1 (MQ=255)
                                                             gcgcCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGc                                                                                 >  7:88562/1‑149 (MQ=255)
                                                              cgcCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGc                                                                                 <  6:354133/148‑1 (MQ=255)
                                                              cgcCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGAAGAAGCCGAATCGTAAAACTCTGTCAGc                                                                                 >  3:100618/1‑148 (MQ=255)
                                                                         acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAg                                                                                                                             >  4:223746/1‑93 (MQ=255)
                                                                         acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAg                                                                                                                             <  3:223746/93‑1 (MQ=255)
                                                                         acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAAt                                                                                                  >  5:64384/1‑120 (MQ=255)
                                                                         acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgcg                                                                     >  4:41575/1‑149 (MQ=255)
                                                                         acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgcg                                                                     >  1:181235/1‑149 (MQ=255)
                                                                         acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgcg                                                                     >  5:290001/1‑149 (MQ=255)
                                                                         acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCCGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgca                                                                     >  2:53586/1‑148 (MQ=255)
                                                                         acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTACGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgcg                                                                     >  5:163207/1‑149 (MQ=255)
                                                                              tCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCgg                                                                <  7:282530/149‑1 (MQ=255)
                                                                                  gTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGatat                                                            >  7:217585/1‑149 (MQ=255)
                                                                                   ttGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATc                                                           <  8:88562/149‑1 (MQ=255)
                                                                                      aGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                        <  1:160200/149‑1 (MQ=255)
                                                                                      aGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                        <  3:127678/149‑1 (MQ=255)
                                                                                        tCTGCGGATCAGCCTCGCCCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                        >  6:273896/1‑147 (MQ=255)
                                                                                           gCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAgcagca                                                   >  4:120988/1‑149 (MQ=255)
                                                                                              gATCATCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCg                                                <  5:115593/149‑1 (MQ=255)
                                                                                                   gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                        <  1:76950/136‑1 (MQ=255)
                                                                                                   gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                        >  2:76950/1‑136 (MQ=255)
                                                                                                    ccTCGACCCAGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAac                                          >  7:121180/1‑149 (MQ=255)
                                                                                                     cTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAaca                                         >  3:67501/1‑149 (MQ=255)
                                                                                                      tCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCg                                                                 <  6:13092/124‑1 (MQ=255)
                                                                                                      tCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCg                                                                 >  5:13092/1‑124 (MQ=255)
                                                                                                           ccTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGc                                                                                 >  3:331338/1‑103 (MQ=255)
                                                                                                           ccTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGc                                                                                 <  4:331338/103‑1 (MQ=255)
                                                                                                               gCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGct                               >  7:279919/1‑149 (MQ=255)
                                                                                                                 tttCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGctct                             >  5:274904/1‑149 (MQ=255)
                                                                                                                   tCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTg                            <  4:67501/148‑1 (MQ=255)
                                                                                                                     cGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCa                                                                               >  4:41433/1‑95 (MQ=255)
                                                                                                                     cGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCa                                                                               <  3:41433/95‑1 (MQ=255)
                                                                                                                      ggATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGACGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCggtctggtct                        >  4:248816/1‑149 (MQ=255)
                                                                                                                         tGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGa                     >  5:78042/1‑149 (MQ=255)
                                                                                                                           gCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGAc                    <  5:273896/148‑1 (MQ=255)
                                                                                                                            catcatCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACgg                  >  6:318695/1‑149 (MQ=255)
                                                                                                                             atcatcAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGc                 <  8:279919/149‑1 (MQ=255)
                                                                                                                                atcaGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCttt              <  6:290001/149‑1 (MQ=255)
                                                                                                                                 tcaGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCtttt             <  2:348426/149‑1 (MQ=255)
                                                                                                                                   aGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGTGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTgg           <  6:163207/149‑1 (MQ=255)
                                                                                                                                     acccgctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGcc         >  2:190694/8‑149 (MQ=255)
                                                                                                                                      tACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCt        >  8:198330/1‑149 (MQ=255)
                                                                                                                                        cTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCt        <  5:318695/147‑1 (MQ=255)
                                                                                                                                        cTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGATTTTGGCCTgg      >  4:213252/1‑149 (MQ=255)
                                                                                                                                         tGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTgg      >  5:166016/1‑148 (MQ=255)
                                                                                                                                          ctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTg       <  8:277919/144‑1 (MQ=255)
                                                                                                                                          ctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTg       >  7:277919/3‑146 (MQ=255)
                                                                                                                                           aGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCGGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGt     >  3:83842/1‑147 (MQ=255)
                                                                                                                                           aGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAActct                                                                                       >  8:115468/1‑65 (MQ=255)
                                                                                                                                           aGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAActct                                                                                       <  7:115468/65‑1 (MQ=255)
                                                                                                                                           aGCAATGGCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGt     <  6:78042/147‑1 (MQ=255)
                                                                                                                                            gCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAg  <  3:120988/149‑1 (MQ=255)
                                                                                                                                             cAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGccc                                                                  >  6:52833/1‑84 (MQ=38)
                                                                                                                                             cAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGccc                                                                  <  5:52833/84‑1 (MQ=38)
                                                                                                                                                   |                                                                                                                                              
GCCCAGACCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAG  >  NC_002947/5059353‑5059641

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: