Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,971,929 A→C 5.1% S255A (TCA→GCA)  tar ← methyl‑accepting chemotaxis protein Tar

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,971,9290AC 5.1% 250.2 / 4.9 79S255A (TCA→GCA) tarmethyl‑accepting chemotaxis protein Tar
Reads supporting (aligned to +/- strand):  ref base A (38/37);  new base C (2/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

AGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCA  >  NC_000913/1971783‑1972078
                                                                                                                                                  |                                                                                                                                                     
aGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAAc                                                                                                                                                    <  2:594359/150‑1 (MQ=255)
  tttCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGc                                                                                                                                                  <  1:320768/150‑1 (MQ=255)
        ccAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgcg                                                                                                                                            >  2:157743/1‑150 (MQ=255)
          aGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGcc                                                                                                                                          >  2:454663/1‑150 (MQ=255)
           gcgcGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGCACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCa                                                                                                                                         <  2:179705/150‑1 (MQ=255)
           gcgcGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCa                                                                                                                                         >  2:19987/1‑150 (MQ=255)
           gcgcGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCa                                                                                                                                         <  1:303749/150‑1 (MQ=255)
            cgcgGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAg                                                                                                                                        <  2:405206/150‑1 (MQ=255)
             gcgGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAgg                                                                                                                                       <  2:799611/150‑1 (MQ=255)
                   gCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGccc                                                                                                                                 >  1:861293/1‑150 (MQ=255)
                    ccTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCa                                                                                                                                >  2:710037/1‑150 (MQ=255)
                    ccTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGGGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCa                                                                                                                                <  1:564341/150‑1 (MQ=255)
                     ctgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCAGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCAt                                                                                                                               >  2:115037/1‑150 (MQ=255)
                         tgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctctttgac                                                                                                                                          >  1:1046250/1‑130 (MQ=255)
                         tgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctctttgac                                                                                                                                          <  2:1046250/135‑6 (MQ=255)
                         tgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCa                                                                                                                           <  2:76631/150‑1 (MQ=255)
                           ttCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGGCGCCCATTTCACt                                                                                                                         <  1:1096583/150‑1 (MQ=255)
                            tCAGTTCGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTg                                                                                                                        >  1:867880/1‑150 (MQ=255)
                                tACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGcc                                                                                                                    <  2:621043/150‑1 (MQ=255)
                                 aCGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGccc                                                                                                                   >  2:933480/1‑150 (MQ=255)
                                  cGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCg                                                                                                                  >  2:438410/1‑150 (MQ=255)
                                   gggAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGt                                                                                                                 >  1:440886/1‑150 (MQ=255)
                                    ggaggaAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTc                                                                                                                >  1:59190/1‑150 (MQ=255)
                                     gaggaAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCa                                                                                                               >  2:1000591/1‑150 (MQ=255)
                                       ggaAAGATCGGTGGTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                                             <  2:1126370/150‑1 (MQ=255)
                                         aaaGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATgg                                                                                                           <  1:454663/150‑1 (MQ=255)
                                           aGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTc                                                                                                         <  2:1243903/150‑1 (MQ=255)
                                            gATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCa                                                                                                        >  2:1168520/1‑150 (MQ=255)
                                            gATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCa                                                                                                        >  1:780980/1‑150 (MQ=255)
                                               cGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAggg                                                                                                     <  2:346911/150‑1 (MQ=255)
                                               cGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAggg                                                                                                     <  1:596190/150‑1 (MQ=255)
                                               cGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAggg                                                                                                     <  2:1192422/150‑1 (MQ=255)
                                                ggTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGt                                                                                                    >  2:963733/1‑150 (MQ=255)
                                                     tGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCg                                                                                               >  2:854811/1‑150 (MQ=255)
                                                     tGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCg                                                                                               <  1:1277108/150‑1 (MQ=255)
                                                      gcccgcccCAATTACACGTGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGc                                                                                              >  2:554159/1‑150 (MQ=255)
                                                       cccGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGcc                                                                                             >  2:630788/1‑150 (MQ=255)
                                                       cccGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGcc                                                                                             >  1:829287/1‑150 (MQ=255)
                                                          gccgcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAgg                                                                                          >  1:38304/1‑150 (MQ=255)
                                                            cgcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGTGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCTCCCATTTCACTGCGCCCGTAAATGGTCAGGTTATTCGCCAGGtt                                                                                        <  2:549383/150‑1 (MQ=255)
                                                            cgcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGtt                                                                                        <  2:59190/150‑1 (MQ=255)
                                                             gcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGtt                                                                                        <  2:42237/149‑1 (MQ=255)
                                                             gcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGtt                                                                                        >  1:42237/1‑149 (MQ=255)
                                                             gcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTa                                                                                       >  1:1170288/1‑150 (MQ=255)
                                                               aaTTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTacc                                                                                     <  2:861293/150‑1 (MQ=255)
                                                                 tttCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTaccac                                                                                   >  1:70145/1‑150 (MQ=255)
                                                                  ttCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTaccacc                                                                                  >  1:804720/1‑150 (MQ=255)
                                                                     aCGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGc                                                                               >  1:299342/1‑150 (MQ=255)
                                                                        ggTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGAt                                                                            <  1:866760/150‑1 (MQ=255)
                                                                         gTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGAtt                                                                           <  1:710037/150‑1 (MQ=255)
                                                                           aCCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTc                                                                         >  1:510633/1‑150 (MQ=255)
                                                                            ccGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTcg                                                                        >  1:436908/1‑150 (MQ=255)
                                                                              ggCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTcgcg                                                                      <  2:804720/150‑1 (MQ=255)
                                                                              ggCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTcgcg                                                                      <  2:48237/150‑1 (MQ=255)
                                                                                  tAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAAtg                                                                  <  1:630788/150‑1 (MQ=255)
                                                                                     aTGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTcgcg                                                                      >  2:1266676/1‑143 (MQ=255)
                                                                                     aTGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTcgcg                                                                      <  1:1266676/143‑1 (MQ=255)
                                                                                        gCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCa                                                            <  1:438410/150‑1 (MQ=255)
                                                                                           tCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCaata                                                         >  1:1065277/1‑150 (MQ=255)
                                                                                               aaCCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAAt                                                                   >  1:51794/1‑136 (MQ=255)
                                                                                               aaCCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAAt                                                                   <  2:51794/136‑1 (MQ=255)
                                                                                                  cTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTg                                                  >  2:248620/1‑150 (MQ=255)
                                                                                                  cTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTg                                                  <  1:59021/150‑1 (MQ=255)
                                                                                                       cgGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGc                                             <  1:19987/150‑1 (MQ=255)
                                                                                                             tGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTa                                       <  1:963733/150‑1 (MQ=255)
                                                                                                                   aCGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGTCAGCGGAGTAAGCAAc                                 >  1:740823/1‑150 (MQ=255)
                                                                                                                   aCGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAAc                                 <  2:70145/150‑1 (MQ=255)
                                                                                                                            gtcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGa                        >  1:729460/1‑150 (MQ=255)
                                                                                                                                 agagCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGcc                   <  1:115037/150‑1 (MQ=255)
                                                                                                                                      gTTTCACATATGCAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTAcc              >  2:507684/1‑150 (MQ=255)
                                                                                                                                      gTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTAcc              >  2:1278618/1‑150 (MQ=255)
                                                                                                                                        tGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCAc            <  2:1170288/150‑1 (MQ=255)
                                                                                                                                         gCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACg           <  2:1065277/150‑1 (MQ=255)
                                                                                                                                         gCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACg           <  2:440886/150‑1 (MQ=255)
                                                                                                                                              tgtgAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGccacc      >  1:977878/1‑150 (MQ=255)
                                                                                                                                               gtgAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTcgcg                                                                      >  2:732474/1‑85 (MQ=255)
                                                                                                                                               gtgAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTcgcg                                                                      <  1:732474/85‑1 (MQ=255)
                                                                                                                                                tgAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACcag    >  2:467534/1‑150 (MQ=255)
                                                                                                                                                 gAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACcagc   >  1:1188110/1‑150 (MQ=255)
                                                                                                                                                  aaaCGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACcagca  <  2:867880/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCA  >  NC_000913/1971783‑1972078

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