Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,116,170 T→C 8.6% L303P (CTG→CCG)  PP_2732 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,116,1700TC 8.6% 164.5 / 5.1 58L303P (CTG→CCG) PP_2732hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/22);  new base C (3/2);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GAAGGCGAGCCGTACCTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATGGCGCAAGCCAT  >  NC_002947/3116022‑3116314
                                                                                                                                                    |                                                                                                                                                
gAAGGCTAGCGTTACCTGTTTCATGCGCGCATCAGTCCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCTACGACGTGGACAAAGCCTGGGGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCt                                                                                                                                                  <  8:16020/149‑1 (MQ=255)
gAAGGCGAGCCGTACCTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGAGGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCt                                                                                                                                                  <  1:187021/149‑1 (MQ=255)
gAAGGCGAGCCGTACCTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGA‑GGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCt                                                                                                                                                  >  2:291996/1‑148 (MQ=255)
    gCGAGCCGTACCTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCCGATc                                                                                                                                              <  4:6594/149‑1 (MQ=255)
      gAGCCGTACCTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCgg                                                                                                                                            <  4:49224/149‑1 (MQ=255)
                          gcgcATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTAcgtgcgtg                                                                                                                        >  8:243936/1‑149 (MQ=255)
                          gcgcATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTAcgtgcgtg                                                                                                                        >  1:337597/1‑149 (MQ=255)
                                cAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATc                                                                                                                   >  3:363246/1‑148 (MQ=255)
                                  gTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGCCGTCTGTGCCACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTAt                                                                                                                <  4:185001/149‑1 (MQ=255)
                                   tGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTAtt                                                                                                               >  2:43391/1‑149 (MQ=255)
                                       gcgcTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCa                                                                                                           <  2:337597/149‑1 (MQ=255)
                                          cTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAgcg                                                                                                        >  2:236623/1‑149 (MQ=255)
                                             aaCATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCAt                                                                                                     >  4:82533/1‑149 (MQ=255)
                                               cATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGc                                                                                                   >  3:84302/1‑149 (MQ=255)
                                                    gCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCa                                                                                              >  8:164903/1‑149 (MQ=255)
                                                      cTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCata                                                                                            >  3:259445/1‑149 (MQ=255)
                                                       ttCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCatat                                                                                           >  2:143530/1‑149 (MQ=255)
                                                          tGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGcc                                                                                        >  2:366365/1‑149 (MQ=255)
                                                                 gCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTc                                                                                  >  5:113513/1‑148 (MQ=255)
                                                                 gCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTc                                                                                  >  5:254969/1‑148 (MQ=255)
                                                                 gCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCa                                                                                 <  7:328951/149‑1 (MQ=255)
                                                                   gtcgtcTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGc                                                                                               >  1:46373/1‑133 (MQ=255)
                                                                   gtcgtcTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGc                                                                                               <  2:46373/133‑1 (MQ=255)
                                                                      gtcTGTGCGACGACGTCGACAAAGCCTGGCGCGA‑GGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCAcc                                                                           <  1:291996/149‑1 (MQ=255)
                                                                        cTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCt                                                                          <  3:82533/149‑1 (MQ=255)
                                                                         tgtgCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTAcgtgcgt                                                                                                                         >  2:352027/1‑101 (MQ=255)
                                                                         tgtgCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTg                                                                         <  3:78291/149‑1 (MQ=255)
                                                                           tgCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCt                                                                          <  1:167083/146‑1 (MQ=255)
                                                                           tgCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGaa                                                                       <  1:206381/149‑1 (MQ=255)
                                                                             cgacgaCGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAAc                                                                      <  6:113513/148‑1 (MQ=255)
                                                                              gacgacGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAGGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACaa                                                                    >  1:161254/1‑149 (MQ=255)
                                                                              gacgacGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACaa                                                                    >  3:324946/1‑149 (MQ=255)
                                                                                    ggcGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCg                                                              <  6:125556/147‑1 (MQ=255)
                                                                                     tCGACAAAGCCTGGCGGGAGGGGCGGGTGCCGCTCAATGCAGCAGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGACAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCgg                                                             <  3:248551/149‑1 (MQ=255)
                                                                                       gACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTGCGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGgc                                                           >  2:197284/1‑149 (MQ=255)
                                                                                       gACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGgc                                                           >  4:331888/1‑149 (MQ=255)
                                                                                          aaaGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAGGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCt                                                        <  2:161254/149‑1 (MQ=255)
                                                                                               cTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGc                                                      >  5:212922/1‑146 (MQ=255)
                                                                                               cTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGc                                                      >  7:301213/1‑146 (MQ=255)
                                                                                                  gcgcGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCa                                                                             >  6:322080/1‑120 (MQ=255)
                                                                                                  gcgcGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGaa                                                 <  1:366365/148‑1 (MQ=255)
                                                                                                  gcgcGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAAt                                                >  2:304961/1‑149 (MQ=255)
                                                                                                     cgAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCt                                             <  4:337587/149‑1 (MQ=255)
                                                                                                      gAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTa                                            <  7:164903/149‑1 (MQ=255)
                                                                                                       aGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTAc                                           <  6:212922/149‑1 (MQ=255)
                                                                                                          ggCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTgg                                        >  8:355409/1‑149 (MQ=255)
                                                                                                                  ccGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATAGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGgcg                                >  7:299708/1‑149 (MQ=255)
                                                                                                                          tGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCa                        >  1:333824/1‑149 (MQ=255)
                                                                                                                           gCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCAt                       >  5:67584/1‑149 (MQ=255)
                                                                                                                               ccGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATgcgc                   >  1:210057/1‑149 (MQ=255)
                                                                                                                                cGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCt                  <  1:236623/149‑1 (MQ=255)
                                                                                                                                 gAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTg                 <  4:259445/149‑1 (MQ=255)
                                                                                                                                  aaGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTg                 >  4:248627/1‑148 (MQ=255)
                                                                                                                                     gCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCGGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATg             <  7:243936/149‑1 (MQ=255)
                                                                                                                                      cTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATgg            <  6:254969/149‑1 (MQ=255)
                                                                                                                                      cTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATgg            >  7:32765/1‑149 (MQ=255)
                                                                                                                                          tcacGCCAGCCGATCAGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCa               >  2:285326/4‑142 (MQ=255)
                                                                                                                                          tcacGCCAGCCGATCAGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCa               <  1:285326/139‑1 (MQ=255)
                                                                                                                                          aTTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATGgcgc         >  5:336414/1‑148 (MQ=255)
                                                                                                                                            acGCCAGCCGATCAGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATGGCGCAAg      >  5:83804/2‑149 (MQ=255)
                                                                                                                                                cAGCCGATCAGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATGGCGCAAGCCAt  >  8:68635/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GAAGGCGAGCCGTACCTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATGGCGCAAGCCAT  >  NC_002947/3116022‑3116314

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