Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,348,798 T→A 5.4% intergenic (+30/+24) 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,7980TA 5.4% 243.6 / 2.5 74intergenic (+30/+24)yfaE/inaAferredoxin‑like diferric‑tyrosyl radical cofactor maintenance protein YfaE/putative lipopolysaccharide kinase InaA
Reads supporting (aligned to +/- strand):  ref base T (40/30);  new base A (2/2);  total (42/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.86e-01

CGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTT  >  NC_000913/2348655‑2348946
                                                                                                                                               |                                                                                                                                                     
cGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCt                                                                                                                                                 >  2:1384599/1‑150 (MQ=255)
  ccTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTa                                                                                                                                               <  2:1311156/150‑1 (MQ=255)
   cTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTAt                                                                                                                                              >  2:394197/1‑150 (MQ=255)
        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:737589/150‑1 (MQ=255)
                         ggATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATgg                                                                                                                       <  2:424956/150‑1 (MQ=255)
                         ggATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATgg                                                                                                                       >  1:1185845/1‑150 (MQ=255)
                          gATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGc                                                                                                                      >  1:699370/1‑150 (MQ=255)
                            ttGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgt                                                                                                                    <  2:733761/150‑1 (MQ=255)
                              gCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCg                                                                                                                                                    <  2:1306829/117‑1 (MQ=255)
                              gCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCg                                                                                                                                                    >  1:1306829/1‑117 (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                                         <  1:46126/150‑1 (MQ=255)
                                                                                                       tCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCAt                                         <  2:282960/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)
                                                                                                                                               |                                                                                                                                                     
CGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTT  >  NC_000913/2348655‑2348946

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: