Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,900,610 C→T W537* (TGG→TGA bglH ← carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,900,6100CT97.9% 142.7 / NA 47W537* (TGG→TGAbglHcarbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/25);  minor base G (0/1);  total (21/26)
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.67e-01

GAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCT  >  NC_000913/3900370‑3900853
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
gaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTc                                                                                                                                                                                                                                             >  2:475287/1‑249 (MQ=255)
 aaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTTAGCAGAAATACCGCTGGCAATATCCATCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCAGCAGATTTCAGc                                                                                                                                                                                                                                          <  2:233896/251‑1 (MQ=255)
  accggaacGCCCACCACTTCATAGACATTTTTTGCTCCCGGGGCCAAGGAGCTAAAGGTCACGCTGTTAGCGGCATTGACCTGAGTCACCTACTGTTTCCTCGGGCTAGTTGGTGCTGTGCTGGCTGGCATATCCGCTGCTATGGCCCACTCAGCAGAAATTCCGCTGTCAATAGCCATCGTTAAACCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTATTTTCCATTCCCATCAGATTTCaccc                                                                                                                                                                                                                                         <  2:453817/251‑4 (MQ=255)
       aacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGc                                                                                                                                                                                                                                    >  1:434660/1‑251 (MQ=255)
                               ttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCtt                                                                                                                                                                                                             <  1:97881/250‑1 (MQ=255)
                                         ggggCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCa                                                                                                                                                                                                  <  1:557800/251‑1 (MQ=255)
                                            gCAAAGTAGCGAAAGGTCAAGCTGTTATCGGCCTTGACCTGCGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATCTCCGCTGCTATGGCCAACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCaa                                                                                                                                                                                                 >  2:305303/1‑249 (MQ=255)
                                                        aaGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATcc                                                                                                                                                                                   >  2:275984/1‑251 (MQ=255)
                                                             cACGCTGGTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTc                                                                                                                                                                              <  1:464707/251‑1 (MQ=255)
                                                                  tGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTCTATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGtttt                                                                                                                                                                         >  2:454249/1‑251 (MQ=255)
                                                                   gTTATCGGCATTGCCCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGACGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGtttt                                                                                                                                                                         <  2:9683/250‑1 (MQ=255)
                                                                    ttATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTcc                                                                                                                                                                       >  1:161050/1‑251 (MQ=255)
                                                                      aTCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAg                                                                                                                                                                     <  1:535054/251‑1 (MQ=255)
                                                                      aTCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTATTCCa                                                                                                                                                                      <  2:297696/250‑1 (MQ=255)
                                                                               gcccTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGat                                                                                                                                                            <  2:387124/249‑1 (MQ=255)
                                                                                aCCTGAGTCACATACTGTTTCACCGGGATAGTTGGAGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAATCAGCAATTTTTATATTCATCGTAATATTACTTTCATTAAGTTCAATTACCTTCATATTTCTGCCTGGGCACCGCCAGCAAAATGGTCGTCTTTCTTGTCTTCAAAGGTGCAGCCATCCAGTTCGTTTTCCAGTGCTTTGat                                                                                                                                                            >  2:187985/1‑250 (MQ=255)
                                                                                 ccTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGatat                                                                                                                                                          >  1:318543/1‑251 (MQ=255)
                                                                                        cacaTACTGTTTCGCCGGGAGAGGGGGTGCCGGGCTGGCTGGGATATCCGGTGCTATGGCCCACTGAGCAGAAAGACCGCTGGCAGCAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTACATGAAGTTCCATTACCGTCAGATTTCAGCCTGGGCACCGACAGCAAAGGGGTCGTCTTTATTTTCTGCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGc                                                                                                                                                   <  1:32514/251‑1 (MQ=255)
                                                                                        cacaTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTGTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTGCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATGGTCTTCAAAGGTGAAGCCAGCCAGTGCGTTTTCCAGGTCTTTGATAGACGTGGc                                                                                                                                                   <  2:371073/251‑1 (MQ=255)
                                                                                             aCTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGa                                                                                                                                              <  1:356680/251‑1 (MQ=255)
                                                                                                   tCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGACATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATAc                                                                                                                                                        >  2:489874/1‑235 (MQ=255)
                                                                                                   tCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCAGGGCACCGACAGCAAACTGATCGGCTTGATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGa                                                                                                                                         <  1:571671/250‑1 (MQ=255)
                                                                                                   tCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTGTATATTCATCGTAATATTCCTTTCATTAAGGTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCGTTGATATAc                                                                                                                                                        <  1:489874/235‑1 (MQ=255)
                                                                                                        gggATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTcc                                                                                                                                   >  2:219588/1‑251 (MQ=255)
                                                                                                                      cGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCa                                                                                                                      >  2:86005/1‑250 (MQ=255)
                                                                                                                           tggctggcATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATAc                                                                                                                >  1:295656/1‑251 (MQ=255)
                                                                                                                           tggctggcATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGAAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCCACATAc                                                                                                                >  2:99810/1‑251 (MQ=255)
                                                                                                                            ggctggcATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTGCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATa                                                                                                                 <  1:149621/249‑1 (MQ=255)
                                                                                                                                     atCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTGAAGCAGCAATTTTTATATTCATCGTAATATTACTTTCATTAAGTTCAATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAAAGAATTTACGGCAGCGGGGTCAACATAATGGAATTGAc                                                                                                      >  1:568548/1‑251 (MQ=255)
                                                                                                                                       ccGCTGCTATGGCCCACGGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATGTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGATATTGACt                                                                                                     <  2:550145/250‑1 (MQ=255)
                                                                                                                                         gctgctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTa                                                                                                  <  1:280433/251‑1 (MQ=255)
                                                                                                                                           tgctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTaa                                                                                                 >  1:243754/1‑250 (MQ=255)
                                                                                                                                            gctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAGTTTTATATTCAGCGTAATATTCCTTTCATTAAGTTCCATTGCCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTTGTTTTCCAGGGTTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGGCACCATACTGGTATTGACTTTaag                                                                                                <  2:111189/250‑2 (MQ=255)
                                                                                                                                                        tGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGc                                                                                   >  1:38642/1‑251 (MQ=255)
                                                                                                                                                                         ggCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCa                                                                    <  1:573856/249‑1 (MQ=255)
                                                                                                                                                                               aGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCAtt                                                            <  1:537012/251‑1 (MQ=255)
                                                                                                                                                                                         aaaGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTGCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAagag                                                                                    <  1:521028/217‑1 (MQ=255)
                                                                                                                                                                                         aaaGCAGAAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACAGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAagag                                                                                    >  2:521028/1‑217 (MQ=255)
                                                                                                                                                                                           agcagcAATTTTTAGATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACt                                                <  2:434660/251‑1 (MQ=255)
                                                                                                                                                                                            gcagcaATTTTTCTATTCATCGTAATATTCCTTTCATTAAGTTACATTACCATCAGATTTCATCCTGGGCACCGACAGCAACCTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGGGGTCAACATACTGGAATTGAATTTAAAGGAATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACtt                                               >  2:313389/1‑250 (MQ=255)
                                                                                                                                                                                               gcaATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTc                                            <  1:532978/251‑1 (MQ=255)
                                                                                                                                                                                                    ttttATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTCTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATc                                       >  1:311721/1‑251 (MQ=255)
                                                                                                                                                                                                    ttttATATTCATCGTAATATTCCTTTAATTAAGTCCCATTACCATCAGATTTCAGCATGGGCACCGAAAGCAACCTGGCCGTCTTTATTGTCTTAACAGGTGATGCCATACAGTTCGTTTTCAAGGGCTTTGCTATAAGTGGCGAAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATAAACTTTATTACTATTCGCATc                                       >  2:558722/1‑251 (MQ=255)
                                                                                                                                                                                                          atTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGt                                  <  1:305303/250‑1 (MQ=255)
                                                                                                                                                                                                                                    gTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCaa       <  1:441573/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCaa       <  1:275984/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aTTACCATCAGATTTCAGCCTGGGCACCGCCAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGAGACGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGCGGTCAACATACTGGTATTGACTTTAAAGGCATGGAAGAGCGTGGTTTTATAACCAGAGTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCCGTTCAACACCt  >  2:441/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
GAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCT  >  NC_000913/3900370‑3900853

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