Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,373,382 C→T A285T (GCA→ACA)  yjeJ ← protein YjeJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,373,3820CT98.0% 192.1 / ‑6.6 51A285T (GCA→ACA) yjeJprotein YjeJ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/23);  total (28/23)
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.99e-01

GCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAATT  >  NC_000913/4373152‑4373606
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
gCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGc                                                                                                                                                                                                                          >  1:548427/1‑239 (MQ=255)
          aTGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTg                                                                                                                                                                                                                <  1:199155/239‑1 (MQ=255)
          aTGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTg                                                                                                                                                                                                                <  1:232695/239‑1 (MQ=255)
          aTGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGa                                                                                                                                                                                                               >  1:22948/1‑240 (MQ=255)
           tgtgATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGaa                                                                                                                                                                                                              <  1:488789/240‑1 (MQ=255)
           tgtgATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGaa                                                                                                                                                                                                              <  1:318263/240‑1 (MQ=255)
            gtgATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGaa                                                                                                                                                                                                              <  2:162094/239‑1 (MQ=255)
              gATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAacac                                                                                                                                                                                                           <  1:402388/240‑1 (MQ=255)
                tCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGt                                                                                                                                                                                                         >  1:683989/1‑240 (MQ=255)
                  cAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGt                                                                                                                                                                                                         <  1:602481/238‑1 (MQ=255)
                    gagGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAg                                                                                                                                                                                                      >  2:418694/1‑239 (MQ=255)
                          gggTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGc                                                                                                                                                                                                >  1:669055/1‑239 (MQ=255)
                           ggTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCa                                                                                                                                                                                               >  2:347540/1‑239 (MQ=255)
                                  aaCGTATGATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGt                                                                                                                                                                                       <  1:612465/240‑1 (MQ=255)
                                        taATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCAt                                                                                                                                                                                 <  1:369051/240‑1 (MQ=255)
                                        taATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCAt                                                                                                                                                                                 <  1:616402/240‑1 (MQ=255)
                                             gCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGt                                                                                                                                                                                                         >  2:502880/1‑211 (MQ=255)
                                             gCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGt                                                                                                                                                                                                         <  1:502889/211‑1 (MQ=255)
                                                   tGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATTACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGt                                                                                                                                                                      >  1:212929/1‑240 (MQ=255)
                                                              tATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTggg                                                                                                                                                           <  2:255862/240‑1 (MQ=255)
                                                                    gCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTAt                                                                                                                                                     >  1:290811/1‑240 (MQ=255)
                                                                    gCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTAt                                                                                                                                                     >  1:426998/1‑240 (MQ=255)
                                                                     ccTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTAtt                                                                                                                                                    <  2:135421/240‑1 (MQ=255)
                                                                      cTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCAt                                                                                                                                                                                 <  2:530633/210‑1 (MQ=255)
                                                                       tACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTc                                                                                                                                                   >  1:489098/1‑239 (MQ=255)
                                                                        aCATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTcg                                                                                                                                                  <  1:459267/239‑1 (MQ=255)
                                                                              ggCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAGGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTAc                                                                                                                                           >  2:651763/1‑240 (MQ=255)
                                                                              ggCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTAc                                                                                                                                           >  2:481512/1‑240 (MQ=255)
                                                                              ggCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTAc                                                                                                                                           >  1:77840/1‑240 (MQ=255)
                                                                              ggCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTAc                                                                                                                                           <  1:347545/240‑1 (MQ=255)
                                                                                   tGCAACTGATTGAATGCTTGGGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATGTCGCCTGTGGTGCAGTTTTGCtg                                                                                                                                                                                                                        >  1:412507/1‑158 (MQ=255)
                                                                                          gATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGCTTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGc                                                                                                                                                                  >  2:215548/1‑205 (MQ=255)
                                                                                          gATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACAt                                                                                                                                <  2:489089/239‑1 (MQ=255)
                                                                                                      gagaATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGa                                                                                                                                                                                                               >  1:67890/1‑148 (MQ=255)
                                                                                                                      gATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGtttt                                                                                                   >  1:221811/1‑240 (MQ=255)
                                                                                                                                ttACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTg                                                                                         >  1:41972/1‑240 (MQ=255)
                                                                                                                                     ccgcATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCcg                                                                                                           <  2:490884/217‑1 (MQ=255)
                                                                                                                                      cgcATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGcaca                                                                                                               >  2:379031/1‑212 (MQ=255)
                                                                                                                                              gCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGc                                                                           >  1:447011/1‑240 (MQ=255)
                                                                                                                                                  gACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCAtt                                                                        >  1:196288/1‑239 (MQ=255)
                                                                                                                                                            cACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCg                                                                      <  2:476436/231‑1 (MQ=255)
                                                                                                                                                                               ggCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTAcc                                           <  2:202602/239‑1 (MQ=255)
                                                                                                                                                                               ggCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTAcc                                           <  1:476868/239‑1 (MQ=255)
                                                                                                                                                                                                  cGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGttt                       >  1:586088/1‑240 (MQ=255)
                                                                                                                                                                                                   gCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGtttt                      >  1:484728/1‑240 (MQ=255)
                                                                                                                                                                                                             aGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGtttt                      <  2:50363/230‑1 (MQ=255)
                                                                                                                                                                                                             aGATGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTgc             >  1:122728/1‑239 (MQ=255)
                                                                                                                                                                                                               aTGAGATTATTTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGAGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGcgcc          >  1:206042/1‑240 (MQ=255)
                                                                                                                                                                                                                      tattTCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAAt   <  2:206039/240‑1 (MQ=255)
                                                                                                                                                                                                                        tttCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAAtt  >  1:349025/1‑239 (MQ=255)
                                                                                                                                                                                                                         ttCGCCTGTGGTGTAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACAACCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGcgcc          >  2:213597/1‑230 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
GCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAATT  >  NC_000913/4373152‑4373606

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