Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,348,812:1 +A 5.6% intergenic (+44/+10) yfaE → / ← inaA ferredoxin‑like diferric‑tyrosyl radical cofactor maintenance protein YfaE/putative lipopolysaccharide kinase InaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,348,8121.A 5.6% 297.5 / 11.6 72intergenic (+44/+10)yfaE/inaAferredoxin‑like diferric‑tyrosyl radical cofactor maintenance protein YfaE/putative lipopolysaccharide kinase InaA
Reads supporting (aligned to +/- strand):  ref base . (38/30);  new base A (2/2);  total (40/32)
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.60e-01

TGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTTTCTCCAGGTCAA  >  NC_000913/2348663‑2348958
                                                                                                                                                      |                                                                                                                                                  
tGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                                                     <  1:493595/150‑1 (MQ=255)
tGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                                                     <  2:570018/150‑1 (MQ=255)
  cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCt                                                                                                                                                  >  1:733761/1‑150 (MQ=255)
   aGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTa                                                                                                                                                 <  2:181342/150‑1 (MQ=255)
        aaGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGGCCAGATAAGACGCATAAGCGTTGTATCTGGCACTACAAcc                                                                                                                                            <  1:518663/150‑2 (MQ=37)
          gTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGAt                                                                                                                                          >  2:737589/1‑150 (MQ=255)
               cTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACAt                                                                                                                                     >  2:1205725/1‑150 (MQ=255)
                 ggATTGCCGAACCGTTAGCCTTTATTCAGCTGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATgg                                                                                                                                   <  2:186162/150‑1 (MQ=255)
                 ggATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATgg                                                                                                                                   >  1:1185845/1‑150 (MQ=255)
                 ggATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATgg                                                                                                                                   <  2:424956/150‑1 (MQ=255)
                 ggATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATgg                                                                                                                                   <  1:737589/150‑1 (MQ=255)
                  gATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGc                                                                                                                                  >  1:699370/1‑150 (MQ=255)
                    ttGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgt                                                                                                                                <  2:733761/150‑1 (MQ=255)
                      gCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtag                                                                                                                              <  2:1116591/150‑1 (MQ=255)
                       ccGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtagt                                                                                                                             >  2:1163525/1‑150 (MQ=255)
                         gAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTaa                                                                                                                           >  1:282960/1‑150 (MQ=255)
                          aaCCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAAt                                                                                                                          <  2:1187355/150‑1 (MQ=255)
                              gTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGc                                                                                                                      >  1:1373615/1‑150 (MQ=255)
                                   ccTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCa                                                                                                                 >  1:281117/1‑150 (MQ=255)
                                       tATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTg                                                                                                             >  2:46126/1‑150 (MQ=255)
                                        aTTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGc                                                                                                            <  1:1384599/150‑1 (MQ=255)
                                           cAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTcc                                                                                                         <  2:1090004/150‑1 (MQ=255)
                                            aGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTccc                                                                                                        <  1:779780/150‑1 (MQ=255)
                                             gCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCCCCCa                                                                                                       >  2:433401/1‑150 (MQ=255)
                                              ccGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAg                                                                                                      <  2:1373615/150‑1 (MQ=255)
                                               cGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGt                                                                                                     >  2:1109600/1‑150 (MQ=255)
                                                  gggAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCg                                                                                                  >  1:1415540/1‑150 (MQ=255)
                                                   ggAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGc                                                                                                 >  1:134068/1‑150 (MQ=255)
                                                     aaaTTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCtt                                                                                               >  2:978942/1‑150 (MQ=255)
                                                       aTTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTg                                                                                             <  1:287346/150‑1 (MQ=255)
                                                        ttttGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                                            >  2:941878/1‑150 (MQ=255)
                                                          ttGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgggg                                                                                          >  1:716367/1‑150 (MQ=255)
                                                               cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTgg                                                                                     <  2:281117/150‑1 (MQ=255)
                                                               cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGACAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTgg                                                                                     >  2:1400807/1‑150 (MQ=255)
                                                                  ttgttgCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTc                                                                                  <  1:1205725/150‑1 (MQ=255)
                                                                    gttgCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGTCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCa                                                                                <  1:674690/150‑1 (MQ=255)
                                                                      tgCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGa                                                                              >  2:50713/1‑150 (MQ=255)
                                                                       gCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGat                                                                             >  2:998459/1‑150 (MQ=255)
                                                                         cGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGatat                                                                           >  1:868583/1‑150 (MQ=255)
                                                                         cGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGatat                                                                           <  2:1415540/150‑1 (MQ=255)
                                                                          gggCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATAtt                                                                          >  1:849669/1‑150 (MQ=255)
                                                                                aaGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGTGGATTGGCTCCAGATATTTTTCCa                                                                    <  2:703420/150‑1 (MQ=255)
                                                                                aaGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGACAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCa                                                                    <  1:1400807/150‑1 (MQ=255)
                                                                                    cGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTg                                                                >  1:434635/1‑150 (MQ=255)
                                                                                       tatTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCg                                                             >  1:428553/1‑150 (MQ=255)
                                                                                         ttGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGa                                                           >  2:210624/1‑150 (MQ=255)
                                                                                            aaaTCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGt                                                        <  1:998459/150‑1 (MQ=255)
                                                                                               tCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGTTCCAGATATTTTTCCAGCTGGCGGAAGTCAt                                                     <  2:1185845/150‑1 (MQ=255)
                                                                                               tCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCAt                                                     <  2:282960/150‑1 (MQ=255)
                                                                                               tCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCAt                                                     <  1:46126/150‑1 (MQ=255)
                                                                                                cGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGCTTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATg                                                    >  2:878359/1‑150 (MQ=255)
                                                                                                 gagaTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCCCCCAGTCTGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGa                                                   <  1:433401/150‑1 (MQ=255)
                                                                                                  agaTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGAt                                                  <  1:394197/150‑1 (MQ=255)
                                                                                                     tgtgAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTaa                                               >  1:335416/1‑150 (MQ=255)
                                                                                                      gtgAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAAt                                              >  2:1329524/1‑150 (MQ=255)
                                                                                                       tgAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACATGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAAta                                             >  2:130771/1‑149 (MQ=255)
                                                                                                        gAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCg                                            >  2:1094415/1‑150 (MQ=255)
                                                                                                           ttGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCtt                                         <  1:50713/150‑1 (MQ=255)
                                                                                                              gTTGTAGGCCAGATAAGACGCATAAGCGTTGTATCTGGCACTACAACGATTACTTGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCACATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTAt                                      >  2:829859/1‑150 (MQ=25)
                                                                                                                 gTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATcgc                                   >  1:666571/1‑150 (MQ=255)
                                                                                                                     ggCCAGATAAGACGCATAAGCGTTGTATCTGGCACTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCt                                                                                   <  1:146834/98‑1 (MQ=21)
                                                                                                                     ggCCAGATAAGACGCATAAGCGTTGTATCTGGCACTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCt                                                                                   >  2:146834/2‑99 (MQ=21)
                                                                                                                      gCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGc                              >  1:619713/1‑150 (MQ=255)
                                                                                                                      gCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGc                              >  2:775346/1‑150 (MQ=255)
                                                                                                                         aGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAg                           <  2:716367/150‑1 (MQ=255)
                                                                                                                              cAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAgcggcg                      >  1:119379/1‑150 (MQ=255)
                                                                                                                               aaCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGAAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGTGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGATTTATAGCGACGCAAGCGgaga                     >  2:596236/1‑150 (MQ=255)
                                                                                                                                   cTTATGCGTCTTATCTGGCCTACAACGATACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACt                >  1:1393593/1‑150 (MQ=255)
                                                                                                                                      aTGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACttt              >  1:68350/1‑150 (MQ=255)
                                                                                                                                        gCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTTtc            <  1:775346/150‑1 (MQ=255)
                                                                                                                                        gCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTTtc            <  1:1163525/150‑1 (MQ=255)
                                                                                                                                           tCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTTTCTcc         <  1:1337467/150‑1 (MQ=255)
                                                                                                                                              tATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTTTCTCCAgg      <  1:1329524/150‑1 (MQ=255)
                                                                                                                                                  tGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTTTCTCCAGGTCaa  >  1:637907/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTTTCTCCAGGTCAA  >  NC_000913/2348663‑2348958

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