Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 612,487 T→C 6.0% K3R (AAG→AGG)  fepA ← ferrienterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913612,4870TC 6.0% 202.2 / 2.5 66K3R (AAG→AGG) fepAferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (38/24);  new base C (2/2);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGTCCGC  >  NC_000913/612340‑612625
                                                                                                                                                   |                                                                                                                                          
gtagAAACGCCAGTCTCATGTAAGTTATGATCGGCGGCTGTAACGACAATAGTATCTTAATGTGAAACAGGAGTATCGGTCTGCTCTTGTGCCTGCGCTACCCAATAAATCCCAAGATTGACCAACAAGGCCAGGGAATGAATcttctt                                                                                                                                           <  1:141058/146‑1 (MQ=255)
    aaaCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTc                                                                                                                                       <  3:324827/149‑1 (MQ=255)
       cGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCAtt                                                                                                                                    <  5:131035/149‑1 (MQ=255)
             gcgcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTa                                                                                                                              <  7:393449/149‑1 (MQ=255)
                    tAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTg                                                                                                                        >  5:347349/1‑148 (MQ=255)
                       gTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAtt                                                                                                                    <  3:107640/149‑1 (MQ=255)
                           tGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGt                                                                                                                                         <  3:269215/124‑1 (MQ=255)
                           tGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGt                                                                                                                                         >  4:269215/1‑124 (MQ=255)
                            gCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGc                                                                                                               >  5:348260/1‑149 (MQ=255)
                             cTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAttttt                                                                                                                 <  3:251973/146‑1 (MQ=255)
                              tCGGCGGCGGTAACGACAATAGTATCGACATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTAACCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTAATTGTTTTATTCCTGCATTTTTGACa                                                                                                             >  8:402770/1‑149 (MQ=255)
                                    gcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAAtt                                                                                                       <  2:123485/149‑1 (MQ=255)
                                         aaCGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTAt                                                                                                                             <  8:241626/122‑1 (MQ=255)
                                         aaCGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTAt                                                                                                                             >  7:241626/1‑122 (MQ=255)
                                            gACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGt                                                                                               >  6:158105/1‑149 (MQ=255)
                                                 tAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTGATTGTTTTATTACTGCATTTTTGCCACGAATTGCAACTGTCggg                                                                                           >  7:38144/1‑148 (MQ=255)
                                                 tAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCggg                                                                                           >  7:211848/1‑148 (MQ=255)
                                                 tAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGACACGAATTGCAACTGTCggg                                                                                           >  7:76450/1‑148 (MQ=255)
                                                  aGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCa                                                                                         >  3:9194/1‑149 (MQ=255)
                                                   gTAGCGGCATGTGCAACAGGCGTATCGGGCGTCTCTTGTGCCTGCGCTCCCCCATACATCCCCAGAGTGCCCAACAAGGCCCGGGAATGAATCTTCTTGTTCATTGTGTTATTCCTGCCTGTTTGCCACGCATTGCAACTGTCGGGcct                                                                                        <  1:43539/149‑3 (MQ=255)
                                                              tgAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTAt                                                                                                                             <  6:341424/101‑1 (MQ=255)
                                                              tgAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTAt                                                                                                                             >  5:341424/1‑101 (MQ=255)
                                                              tgAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAttttt                                                                                                                 >  3:58083/1‑113 (MQ=255)
                                                              tgAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAttttt                                                                                                                 <  4:58083/113‑1 (MQ=255)
                                                               gAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTg                                                                                                      >  4:255185/1‑123 (MQ=255)
                                                               gAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTg                                                                                                      <  3:255185/123‑1 (MQ=255)
                                                                aaaCATGAGTATCTCTCGTCTCTTGTGCCTTCGCTACCCCATAAATCCCAAGATTGAACAAAAATCCCAGGGAATTAATCTTATTGTTCATTGTTTTATTCCTTCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCaa                                                                            <  5:110851/148‑1 (MQ=255)
                                                                  aCAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAacac                                                                         <  6:348260/149‑1 (MQ=255)
                                                                        gTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                                   >  7:223297/1‑149 (MQ=255)
                                                                         tATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTTTCGGGCATGGTCGTCATCAACACGACGCAt                                                                  >  7:160493/1‑149 (MQ=255)
                                                                              gtcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGc                                                             >  6:158440/1‑149 (MQ=255)
                                                                               tcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                                      <  6:302532/139‑1 (MQ=255)
                                                                               tcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                                      >  5:302532/1‑139 (MQ=255)
                                                                               tcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGc                                                             >  4:204073/1‑148 (MQ=255)
                                                                                  gCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                         >  7:18894/1‑149 (MQ=255)
                                                                                      ttGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcgcg                                                      >  6:299754/1‑148 (MQ=255)
                                                                                      ttGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCCTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcgcg                                                      >  7:36275/1‑148 (MQ=255)
                                                                                      ttGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCCTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcgcg                                                      >  5:374825/1‑148 (MQ=255)
                                                                                        gtgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaa                                                    >  8:346360/1‑148 (MQ=255)
                                                                                         tgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaaaa                                                  >  2:377791/1‑149 (MQ=255)
                                                                                            cTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCt                                               >  1:291994/1‑149 (MQ=255)
                                                                                             tGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCtt                                              >  8:245250/1‑149 (MQ=255)
                                                                                              gcgcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                             >  4:15709/1‑149 (MQ=255)
                                                                                               cgcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                             <  5:299754/148‑1 (MQ=255)
                                                                                               cgcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCtt                                              <  5:158105/147‑1 (MQ=255)
                                                                                                 cTACCCAATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTTTTTTATTCCTTCATTTTTGCCACTAATTGCATCTGTCGGTCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGAt                                          <  6:44157/149‑1 (MQ=255)
                                                                                                    ccccATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATc                                         <  4:364530/147‑1 (MQ=255)
                                                                                                    ccccATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCt                                       >  7:193733/1‑149 (MQ=255)
                                                                                                       cATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGa                                     >  5:61217/1‑148 (MQ=255)
                                                                                                        aTAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTATTATTCCTCCATTTTTGCCACGAATTGc                                                                                                     >  1:48283/1‑83 (MQ=255)
                                                                                                               cccAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACg                             <  1:377791/148‑1 (MQ=255)
                                                                                                                ccAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGACATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGa                                     >  8:184810/1‑139 (MQ=255)
                                                                                                                ccAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGACATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGa                                     <  7:184810/139‑1 (MQ=255)
                                                                                                                ccAGATTGAACAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCGTCATTTTTGCCACGAATTGCAACTGTCGGGCATTGTCGTCATCAACAAGAAGCATCCCGCTACCGCTAAAACCTTGGATCCTGAAAGACACgaa                           >  6:220800/1‑147 (MQ=255)
                                                                                                                          caacaaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTgg                 >  1:183009/1‑149 (MQ=255)
                                                                                                                          caacaaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTATGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTgg                 >  1:85670/1‑149 (MQ=255)
                                                                                                                            acaaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCTAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGtt               <  2:291994/149‑1 (MQ=255)
                                                                                                                            acaaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGt                >  5:381667/1‑148 (MQ=255)
                                                                                                                                   cAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTtaattaat        <  3:204073/149‑1 (MQ=255)
                                                                                                                                    aGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATg       >  2:224692/1‑149 (MQ=255)
                                                                                                                                    aGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATg       >  7:296632/1‑149 (MQ=255)
                                                                                                                                    aGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATg       >  4:309301/1‑149 (MQ=255)
                                                                                                                                      ggAATGAATCTTCCTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGTc     <  6:374825/149‑1 (MQ=255)
                                                                                                                                      ggAATGAATCTTCCTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGTc     <  8:36275/149‑1 (MQ=255)
                                                                                                                                         atgaatCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                             >  3:21374/1‑106 (MQ=255)
                                                                                                                                         atgaatCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                             <  4:21374/106‑1 (MQ=255)
                                                                                                                                         atgaatCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGTCcgc  >  5:259690/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
GTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGTCCGC  >  NC_000913/612340‑612625

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: