Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 937,458 C→T S29S (AGC→AGT lolA → outer membrane lipoprotein carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913937,4580CT97.9% 151.2 / NA 48S29S (AGC→AGTlolAouter membrane lipoprotein carrier protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/17);  minor base G (1/0);  total (31/17)
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.86e-01

AATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACTGAAAATTCAGTATTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGCTACGGCAACCTGGCTGAAAGATGCCACCGGTAATACGCCGT  >  NC_000913/937216‑937702
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
aaTGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACTGAAAATTCAGTATTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATGTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGGTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGa                                                                                                                                                                                                                                              <  2:98355/251‑1 (MQ=255)
             ccGGATAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACTGAAAATTCAGTATTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCagcagc                                                                                                                                                                                                                                 >  2:133396/1‑251 (MQ=255)
              cgtaTAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACTGAAAATTCAGTATTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCt                                                                                                                                                                                                                                <  1:292333/248‑1 (MQ=255)
                                gacGCATCCGGCACTCTATCAACTGCAAATTCCGTATTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGCTCGGGAGCGACGTAGTTTGAGGAATACTGATGAAAACAATTGCAATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTcnccc                                                                                                                                                                                                              >  2:584878/3‑249 (MQ=255)
                                gTCGCATCCGGCACTCTATCAACTGAAAATTCAGTATTTTCTTCATTCCTCAAGCTGATTATTAGCCTGGAATAGAGCGCAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGACGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTcacac                                                                                                                                                                                                              >  1:36574/1‑251 (MQ=255)
                                     aTCCGGCACTCTATCAACTGAAAATTCAGTATTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACaaaa                                                                                                                                                                                                          >  1:337357/1‑250 (MQ=255)
                                          gccctctATCAACTGAACCTTCAGTATTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGACAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACACAAAGtgac                                                                                                                                                                                                    >  2:418402/4‑251 (MQ=255)
                                           cACTCTATCAACTGAAAATTCAGTATTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGtgact                                                                                                                                                                                                   >  2:90886/1‑251 (MQ=255)
                                                          aacttCAGTATTTTCTTCTTTCCACAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTAGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGCATGACGGTAgcggcgcg                                                                                                                                                                                    >  2:426803/4‑251 (MQ=255)
                                                                   aTTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAGGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAg                                                                                                                                                                           <  1:205588/251‑1 (MQ=255)
                                                                               ccTCAAGCTTATTATTAGCCTGGAATAGAGAGTAGAGGGAACGCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCGCCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCTGCGCGGCGGTGCAGGAAGGTCAggg                                                                                                                                                                <  2:429006/250‑1 (MQ=255)
                                                                               ccTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGc                                                                                                                                                               >  1:503468/1‑251 (MQ=255)
                                                                                         attattAGCCTGGAATAGAGAGTAGAGGGAAGAAAGGAAGGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGTCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCACGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTAACACAAAAAGTGACTGAAGGTAGCGGCGCGGCGGTGCAGGGAGGTCAGGGTGATcgctggg                                                                                                                                                     >  1:579952/1‑245 (MQ=255)
                                                                                            attaGCCTGGAATAGCGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTTCACGAAGGTCCGGGCGGTCTGTTGGTg                                                                                                                                                   >  1:554244/1‑250 (MQ=255)
                                                                                                 ccTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGGCAGGGCGATCTGTGGGTGAAACgg                                                                                                                                             >  1:201943/1‑250 (MQ=255)
                                                                                                      aaTAGAGAGTAGAGGGAACTCCCGATCCGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGCAGGACAAGGCGCTCTGTGGGGGAAAAGTCCaa                                                                                                                                         >  1:453446/1‑250 (MQ=255)
                                                                                                         agagagTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTa                                                                                                                                    <  2:106825/251‑1 (MQ=255)
                                                                                                            gagTAGAGGGAACTCCCGCTCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCGGGGAGGTAAGGGCGATCTGTGGGTGAAACGGCCAAACTTAtt                                                                                                                                  >  2:153391/1‑251 (MQ=255)
                                                                                                                 gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACt                                                                                                                             <  2:307100/251‑1 (MQ=255)
                                                                                                                    ggAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGc                                                                                                                          >  1:192140/1‑251 (MQ=255)
                                                                                                                        ttCCCGATCGGGAGTGCCGTAATTTGAGGAATAATGATGAAAAAACTTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGCTAACGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCCGGGCGATCTGTGGGTGAAACGTCCCAACTTATTCAACTAGCatat                                                                                                                      >  2:473886/2‑251 (MQ=255)
                                                                                                                           ccGATCGGGAGTGCCGTAATTTGAGCAATAATGATGACAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGCCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCCGGGCGATCTGTGGGTGAAACGCCCAAACTTATTCCACTGGCATATGac                                                                                                                   >  2:34039/1‑251 (MQ=255)
                                                                                                                                     gTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCtgat                                                                                                         <  1:42513/251‑1 (MQ=255)
                                                                                                                                           tAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAAc                                                                                                                              >  1:79063/1‑224 (MQ=255)
                                                                                                                                           tAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAAc                                                                                                                              <  2:79063/224‑1 (MQ=255)
                                                                                                                                                     aataatGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAATCTTATTCAACTGGCATATGACACAACCTTATGAAAGCATTCTGtttt                                                                                         <  2:33257/251‑1 (MQ=255)
                                                                                                                                                            atgaAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACgg                                                                                  <  1:659286/251‑1 (MQ=255)
                                                                                                                                                            atgaAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACgg                                                                                  >  1:490460/1‑251 (MQ=255)
                                                                                                                                                                          cATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGt                                                                    >  1:174808/1‑251 (MQ=255)
                                                                                                                                                                             cACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCt                                                                 <  2:428533/251‑1 (MQ=255)
                                                                                                                                                                             cACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCAGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGAGACTGACGGTAGCGGCGCGGCGGGGCCGGAAGGGCAGGGAGATCTGTGGGGGAAACGGCCAAAATTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAAAACTGTGGTTCt                                                                 >  1:400966/1‑251 (MQ=255)
                                                                                                                                                                              aCCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGGAAAACACTGTGGTTCt                                                                 >  2:447197/1‑250 (MQ=255)
                                                                                                                                                                               ccTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCtat                                                               >  1:136977/1‑251 (MQ=255)
                                                                                                                                                                                 tgtgCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCt                                                                                                                                                          >  2:166977/1‑158 (MQ=255)
                                                                                                                                                                                 tgtgCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCt                                                                                                                                                          <  1:166977/158‑1 (MQ=255)
                                                                                                                                                                                  gtgCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATaa                                                             <  2:174808/250‑1 (MQ=255)
                                                                                                                                                                                              cAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCcgttcgtt                                                   <  2:136977/248‑1 (MQ=255)
                                                                                                                                                                                              cAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAAAACTGTGGTTCTATAAACCGTTCGTTg                                                  >  2:352661/1‑249 (MQ=255)
                                                                                                                                                                                                 gCTTAGTAGCAAGCAGCGTTTGAGCTGATGCCGCAAGCGATCTGAAAAGGCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGGGCAGGAAAGGCAGGGCGATCTGTGGGTGAAACGGCCAAACTTATTCAACTGGAATATGACACAACCTGATAAAAGCCTTATGGGTTCTGAAGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAgga                                              >  2:449013/1‑248 (MQ=255)
                                                                                                                                                                                                     agtagCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGGGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGCt                                          >  1:162935/1‑250 (MQ=255)
                                                                                                                                                                                                           aagcaGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGGAGGAAAGTCAGGGCGGTCTGAGGGGGAAACGACCAAACTTATTCAAATGGCAGATGACAAAACCCGATGGAAGCATTCTGGTTTCAGACGGTAAAACACTGTTGGTCTATAACCCGTTCGTTGAGCAAGCCACGGCa                                    >  2:495601/1‑250 (MQ=255)
                                                                                                                                                                                                             gcagcGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGGGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGCTACGGCAAcc                                 >  2:297461/1‑251 (MQ=255)
                                                                                                                                                                                                                gcGTTTGGGCTGAGGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTGGTTGAGCAAGCTACGGCAACCTgg                              <  2:368187/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGCTACGGCAACCTGGc                             <  2:636856/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cGATCTGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGGCAGGGCGATCTGGGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGCTACGGCAACCTGGCTGAAAGATGCCACCGGTAAt         >  2:435423/1‑250 (MQ=255)
                                                                                                                                                                                                                                           tGAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAAAACTGTGGTTCTATAACCCGTTCGTTGAGCAAGCTACGGCAACCTGGCTGAAAGATGCCACCGGTAATACGCCg   >  1:561351/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGCAGGTCAGGGCGCTCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTTTGGTTTCTGACGGTAAAACACTGTGGTTATATAACCCGTTCGTTGAGCAAGACACGGCAACCTGGCTGAAAGATGCCACCGGTAATACGCCGt  >  2:267043/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gAAAAGTCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGCTACGGCAACCTGGCTGAAAGATGCCACCGGTAATACGc     <  1:153391/248‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
AATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACTGAAAATTCAGTATTTTCTTCTTTCCTCAAGCTGATTATTAGCCTGGAATAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACGGTAGCGGCGCGGCGGTGCAGGAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGCTACGGCAACCTGGCTGAAAGATGCCACCGGTAATACGCCGT  >  NC_000913/937216‑937702

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