Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,224,429 C→A 6.0% intergenic (‑541/+280) PP_3703 ← / ← PP_3704 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,224,4290CA 6.0% 208.7 / 3.6 67intergenic (‑541/+280)PP_3703/PP_3704hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base A (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

TCAGCTTGCAGCTTACAGAAGACCTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCG  >  NC_002947/4224283‑4224575
                                                                                                                                                  |                                                                                                                                                    
tCAGCTTGCAGCTTACAGAAGCCCTCGAAATTCGCACTGGGTGAGCGAGGGAGAGTCCGTATGACAGATCGTCGGGGTATCACTGCAGCGGTAAATGTCGCTACGCCTGGCGGCGGGGAGCGGCCTCGCTCCCGTAGCAGAggcggc                                                                                                                                                      >  4:180631/1‑147 (MQ=255)
    cTTGCAGCTTACAGAAGACCTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGACAGATa                                                                                                                                                >  1:362267/1‑149 (MQ=255)
    cTTGCAGCTTACAGAAGACCTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGACAGATa                                                                                                                                                >  1:367846/1‑149 (MQ=255)
        cAGCTTACAGAAGACCTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTg                                                                                                                                            <  1:347118/149‑1 (MQ=255)
                agaagaCCTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATTGTGGCCCTATGAATACAAACCTAAAAGTTCCTAAGCCATGCCGCCGGGGACAACCTAGCTACCCTTACAAAGACGGCCAGAAACCTTATAAGCgg                                                                                                                                    >  1:284504/1‑149 (MQ=255)
                  aagaCCTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAg                                                                                                                                  <  8:172633/149‑1 (MQ=255)
                   agaCCTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGACAGATACGTGATAAGCGGAg                                                                                                                                  <  2:367846/148‑1 (MQ=255)
                   agaCCTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGACAGATACGTGATAAGCGGAg                                                                                                                                  <  2:362267/148‑1 (MQ=255)
                       cTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACg                                                                                                                                              >  5:44201/1‑132 (MQ=255)
                       cTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACg                                                                                                                                              <  6:44201/132‑1 (MQ=255)
                       cTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATg                                                                                                                             >  3:184330/1‑149 (MQ=255)
                          gAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCgg                                                                                                                          >  4:111865/1‑149 (MQ=255)
                          gAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCgg                                                                                                                          <  2:66090/149‑1 (MQ=255)
                                cacaCTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGaaaa                                                                                                                     <  6:36607/148‑1 (MQ=255)
                                   acTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAAcccc                                                                                                                 >  5:215066/1‑149 (MQ=255)
                                       gttgAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCcgc                                                                                                               <  3:308136/147‑1 (MQ=255)
                                           aaCGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGaaaa                                                                                                         <  8:82250/149‑1 (MQ=255)
                                              gATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATaa                                                                                                      <  3:111865/149‑1 (MQ=255)
                                              gATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATa                                                                                                       >  2:357449/1‑148 (MQ=255)
                                               aTGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATaaa                                                                                                     >  8:104219/1‑149 (MQ=255)
                                                tGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAAc                                                                                                    <  4:341739/149‑1 (MQ=255)
                                                     aTTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAgg                                                                                               <  2:284504/149‑1 (MQ=255)
                                                         cGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGc                                                                                           >  7:148637/1‑149 (MQ=255)
                                                             tGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑Ccag                                                                                     >  1:343809/1‑149 (MQ=255)
                                                               aCAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑Ccag                                                                                     >  7:217880/1‑147 (MQ=255)
                                                                    tcgtcgCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCa                                                                              <  4:184330/149‑1 (MQ=255)
                                                                         gCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCtgt                                                                         >  6:274134/1‑149 (MQ=255)
                                                                           ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCAGTAACAAAGACGGCCAGATAAGTGATAAGCGGAGCGATGCGGGAAAACCCCGAGGAAAATAAACATAGGTGGC‑‑CCAGCAGTGCACCtgtg                                                                        >  5:376390/1‑148 (MQ=255)
                                                                            cTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTa                                                                                                 >  7:76910/1‑124 (MQ=255)
                                                                            cTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTa                                                                                                 <  8:76910/124‑1 (MQ=255)
                                                                            cTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCa                                                                      <  7:99428/149‑1 (MQ=255)
                                                                                    ccAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGt                                                              >  5:380035/1‑149 (MQ=255)
                                                                                     cacaCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGtt                                                             <  4:149606/146‑1 (MQ=255)
                                                                                     cAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGTAAAATAAACCTAGGTGGC‑‑Ccagca                                                                                   <  4:155326/127‑1 (MQ=255)
                                                                                     cAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑Ccagca                                                                                   >  3:155326/1‑127 (MQ=255)
                                                                                      aaaCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTTGGTGGCAACCAGCAATGCACAGGTGCACTGCTggg                                                              <  1:56768/149‑2 (MQ=255)
                                                                                      aaaCGTAAATGTCCCTAAGCCTGGAGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTg                                                            >  5:341750/1‑149 (MQ=255)
                                                                                       aaCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGc                                                           <  6:341750/149‑1 (MQ=255)
                                                                                             aaTGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACgtg                                                     <  6:380035/149‑1 (MQ=255)
                                                                                                 tCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                 <  8:148637/149‑1 (MQ=255)
                                                                                                 tCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                 <  2:343809/149‑1 (MQ=255)
                                                                                                 tCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTa                                                   >  8:219496/1‑147 (MQ=255)
                                                                                                   ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAgg                                               >  1:12876/1‑149 (MQ=255)
                                                                                                   ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAgg                                               >  2:314578/1‑149 (MQ=255)
                                                                                                    cTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                              >  7:75532/1‑149 (MQ=255)
                                                                                                         ccTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                          >  6:23403/1‑148 (MQ=255)
                                                                                                         ccTGGCGGCGGGGAACAACCACGCTCACGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTCGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                          >  7:207682/1‑148 (MQ=255)
                                                                                                                  ggggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGa                                >  8:143083/1‑149 (MQ=255)
                                                                                                                  ggggAACAACCACGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTCGTTGCCACGTGTAAAGGCTTATGAGCTAAAGa                                >  4:269513/1‑149 (MQ=255)
                                                                                                                   gggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCGCCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGAc                               <  7:104219/149‑1 (MQ=255)
                                                                                                                     gAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAg                             <  6:215066/149‑1 (MQ=255)
                                                                                                                         aacCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACt                         >  8:160251/1‑149 (MQ=255)
                                                                                                                          cccTCGCTCCCGTAACAAAGACGGCCAGATTCGTGATAAGCGGAGCGATGCGGGAAAACCGCGCGGAAAATAAACCTAGGGGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGTCTg                        <  8:207682/148‑1 (MQ=255)
                                                                                                                           ccTCTCTCCATTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                       <  3:180631/149‑1 (MQ=255)
                                                                                                                           ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                       >  8:41941/1‑149 (MQ=255)
                                                                                                                             tCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                       <  1:220209/147‑1 (MQ=255)
                                                                                                                             tCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                       >  2:220209/1‑147 (MQ=255)
                                                                                                                              cGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTa                    >  6:294043/1‑149 (MQ=255)
                                                                                                                                    cGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCgcgc              <  1:314578/149‑1 (MQ=255)
                                                                                                                                     gTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGcc             >  4:146132/1‑149 (MQ=255)
                                                                                                                                        aCAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGa                                        >  6:165708/1‑119 (MQ=255)
                                                                                                                                        aCAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCa           <  7:143083/148‑1 (MQ=255)
                                                                                                                                         cAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATc         >  6:297122/1‑149 (MQ=255)
                                                                                                                                         cAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCActgctg                                                                <  4:104799/94‑1 (MQ=39)
                                                                                                                                         cAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCActgctg                                                                >  3:104799/1‑94 (MQ=39)
                                                                                                                                          aaaGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCt        <  5:23403/149‑1 (MQ=255)
                                                                                                                                           aaGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTACAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTg       >  4:241675/1‑149 (MQ=255)
                                                                                                                                                ggCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCGCATCTGTAGCg  >  4:203990/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCAGCTTGCAGCTTACAGAAGACCTCGAAATTCACACTTGTTGAACGATGGAGATTCCGTATGACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCG  >  NC_002947/4224283‑4224575

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: