Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,873 T→A 12.8% intergenic (+90/+62) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,8730TA12.8% 1726.7 / 56.9 639intergenic (+90/+62)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (153/376);  new base A (91/2);  total (258/381)
Fisher's exact test for biased strand distribution p-value = 9.11e-40
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 699 total aligned reads displayed.

TCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCT  >  NC_000913/3504724‑3505022
                                                                                                                                                     |                                                                                                                                                      
tcCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCct                                                                                                                                                        <  2:2596596/150‑1 (MQ=255)
   gAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCtc‑tc                                                                                                                                                    <  1:3606613/150‑1 (MQ=255)
       tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCCCCCGAAACCAC‑Accccc                                                                                                                                                >  1:14126388/1‑149 (MQ=39)
         ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCaaa                                                                                                                                              <  1:128595/150‑1 (MQ=255)
         ttctGCCGCGACAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCaaa                                                                                                                                              <  2:7463355/150‑1 (MQ=255)
          tctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAACCAC‑ACCCCaaac                                                                                                                                             >  1:13585060/1‑149 (MQ=39)
            tGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAggg                                                                                                                                           <  1:5324960/150‑1 (MQ=255)
             gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAgggg                                                                                                                                          <  2:12421358/150‑1 (MQ=255)
             gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCTAACCCTC‑ACCCCAAAgggg                                                                                                                                          >  2:15246251/1‑150 (MQ=255)
              ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGc                                                                                                                                         <  2:1081324/150‑1 (MQ=255)
               cgcgCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAACCACAACCCCAAAgggcg                                                                                                                                         >  1:658877/1‑148 (MQ=38)
                gcgcCAATCGCATTACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGa                                                                                                                                       <  1:8891835/150‑1 (MQ=255)
                 cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTACCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAg                                                                                                                                      <  1:3030589/150‑1 (MQ=255)
                 cgcCAAGCGCATCACCTTACCGATGGATAATAAAAGACCTTAACCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACAGCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCACAGGGGCGAg                                                                                                                                      <  2:4455665/150‑1 (MQ=38)
                  gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGGTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAgg                                                                                                                                     <  2:3358426/150‑1 (MQ=255)
                  gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAgg                                                                                                                                     <  1:7511633/150‑1 (MQ=255)
                  gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑GCCCCTTTGGGGAGAgg                                                                                                                                     <  2:13916195/150‑1 (MQ=39)
                    cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAgggg                                                                                                                                   <  2:654060/150‑1 (MQ=255)
                    cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAgggg                                                                                                                                   <  2:3741384/150‑1 (MQ=255)
                     aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGa                                                                                                                                  <  1:15049420/150‑1 (MQ=255)
                        cGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCg                                                                                                                               <  2:7767992/150‑1 (MQ=255)
                          cATCACCTTACCTCTTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑GCCCCTTTGGGGAGAgggttagggt                                                                                                                             <  2:9407100/150‑8 (MQ=38)
                            tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAACGGGGCGAGGGGACCGTTCa                                                                                                                           <  1:3423569/150‑1 (MQ=255)
                               ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                        <  1:728522/150‑1 (MQ=255)
                               ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                        <  2:10226083/150‑1 (MQ=255)
                                   aaCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGa                                                                                                                    <  2:13597825/150‑1 (MQ=255)
                                    aCTATTGATACTAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCACAGGGGCGAGGGGACCGTTCAGTACGGaa                                                                                                                   <  2:4913164/150‑1 (MQ=255)
                                    aCTAGTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGaa                                                                                                                   <  2:3925365/150‑1 (MQ=255)
                                      tagtGATAATAAAAGACATTAACCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAtt                                                                                                                 <  2:9727552/147‑1 (MQ=255)
                                      tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCAC‑ACCCCCAAGGGGCGCGGGGCCCGGTCAGTACCGaaat                                                                                                                 >  2:14701579/1‑148 (MQ=38)
                                        ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCAC‑TACCCCAAGGGGCGCGGGGCCCCGGGAGTACGGAAgttt                                                                                                               >  2:14552611/1‑149 (MQ=37)
                                         tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                              <  1:4923835/150‑1 (MQ=255)
                                         tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                              <  2:15575702/150‑1 (MQ=255)
                                         tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                              <  1:13171318/150‑1 (MQ=255)
                                           ataataAAAGACATTTACCGTTGGCCGGATAAGATGCGCCAGCTTCGCATGCGGCGATGGTGGCGCGTTGAATTTTTCATCCCGTACGTTCCCCTCTCGCTTAGCCTC‑TCCCCTATGGGGCGTGGGGACGGTTCAGTACGGTATTTTGCt                                                                                                            >  2:15480765/1‑150 (MQ=255)
                                           ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCACATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAATCCTC‑TCCCCaaa                                                                                                                                              >  1:2453769/1‑116 (MQ=255)
                                             aataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCAAAACCAC‑AACCCCAAGGGGCGGGGGGGCCCGTCACTACGGACATTATCTTg                                                                                                          >  1:9423715/1‑150 (MQ=25)
                                              ataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCACATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAAACCAC‑TACCCAAATGGGGGAGGGGGCCGTTAAGTACGGAATTTATCtttt                                                                                                         >  1:3847966/1‑148 (MQ=38)
                                                  aaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGa                                                                                                                                       <  1:966757/116‑1 (MQ=255)
                                                   aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                    <  2:5371527/150‑1 (MQ=255)
                                                   aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                    <  1:14039901/150‑1 (MQ=255)
                                                   aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTCCGGAATTTAGCTTGACGGTg                                                                                                    >  2:8937757/1‑150 (MQ=255)
                                                       aTTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCAC‑TCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTAt                                                                                                >  2:2461091/1‑150 (MQ=255)
                                                        ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATc                                                                                               >  2:13012155/1‑150 (MQ=255)
                                                        cgacccGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATc                                                                                               <  2:10391723/146‑1 (MQ=255)
                                                         tAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATcc                                                                                              <  2:3127903/150‑1 (MQ=255)
                                                          aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                             <  1:12537897/150‑1 (MQ=255)
                                                             cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                          <  2:11410285/150‑1 (MQ=255)
                                                                  gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                                     <  2:13373353/150‑1 (MQ=255)
                                                                    cGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATa                                                                                   <  2:5722798/150‑1 (MQ=255)
                                                                      gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACt                                                                                 <  1:14479547/150‑1 (MQ=255)
                                                                       aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                <  1:15307530/150‑1 (MQ=255)
                                                                       aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                <  2:13703757/150‑1 (MQ=255)
                                                                        tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                               <  2:9251980/150‑1 (MQ=255)
                                                                           gATGCGCCAGCATCGCATCCTGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCTAAGGGGCGAGTGGACCGTTCAGTACGGAATTGTGCTTGACGGTGCTATCCAGCATTAAATACTGGCaa                                                                            <  2:8175526/150‑1 (MQ=39)
                                                                           gATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCAACCGAAACCAC‑AACCCAAAGGGGCGAGGGGGCCCTTCAGTACGGAATTTAGCATGACGGTGCTATACAGCATTAAATAATGGcca                                                                            >  2:219849/1‑148 (MQ=255)
                                                                              gcgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAc                                                                         <  1:6135380/150‑1 (MQ=255)
                                                                              gcgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAc                                                                         <  2:14843430/150‑1 (MQ=255)
                                                                                gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATActc                                                                       <  1:7945492/150‑1 (MQ=255)
                                                                                    gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                   <  2:11838388/150‑1 (MQ=255)
                                                                                          catcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAACCAC‑ACCCCCAACGGGCGAGGGGACCGGTCAGTACGGCATTTAGATTGACGGTGCGATACAGCAATAAAAAATGGCAAAACACAAACCCaaaa                                                             >  1:11495570/1‑150 (MQ=255)
                                                                                           atcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGTCCGTTGAGTACGGAATTTAGCTTGTCGGTGCTATCCAGCATTAAATAATGGCAATAATCTCACGCAAAAt                                                            >  2:11030191/1‑150 (MQ=39)
                                                                                           atcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAAt                                                            <  2:11601467/150‑1 (MQ=255)
                                                                                               ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                        <  1:318594/150‑1 (MQ=255)
                                                                                                       tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGt                                                <  1:8458035/150‑1 (MQ=255)
                                                                                                        ggCGCGTTGAATTTTACCTCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtc                                               <  1:10026267/150‑1 (MQ=255)
                                                                                                             gTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCAC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTc                                          <  2:3515455/150‑1 (MQ=255)
                                                                                                              ttGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCa                                         <  2:5952212/150‑1 (MQ=255)
                                                                                                               tGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                                        <  1:2042230/150‑1 (MQ=255)
                                                                                                               tGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                                        <  1:8683877/150‑1 (MQ=255)
                                                                                                               tGAATTTTACATCCCGTACGGTCCCCTCACCCAAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                                        <  1:11090746/150‑1 (MQ=255)
                                                                                                                   ttttACATCCCGTACGTTCCCCTCACCCTAACCCAC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTAACGCAAAATCCGCCGCACCGTCTCTTCACCGcg                                    >  2:1861198/1‑149 (MQ=255)
                                                                                                                     ttACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGc                                  <  1:10451626/150‑1 (MQ=255)
                                                                                                                     ttACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGc                                  <  2:5607192/150‑1 (MQ=255)
                                                                                                                     ttACATCCCGTACGTTCCCCGCACCCTAACCCAC‑TACCCAAAGGGGCGGCGGGGCCGGTCAGTACGGAATATAGCTTGACGGTGCTATCCAGCAATAAATACTGGCAAAAAACTCACGCACAATCCCCCGCACCCTATCTTCACCGCTGc                                  >  2:8497241/1‑150 (MQ=255)
                                                                                                                          tCCCGTACGGTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTc                             <  1:14755758/150‑1 (MQ=255)
                                                                                                                           cccGTACGTTCCCCTCACCCTACCCCTC‑TCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCg                            <  1:4174480/150‑1 (MQ=255)
                                                                                                                           cccGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCg                            <  2:666584/150‑1 (MQ=255)
                                                                                                                           cccGTACGTTCCCCTAAACCTAAACCTA‑TACCCAAATGGGCGAGGGGACCGTTAAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAACGCTGCGATTCg                            >  2:2307394/1‑150 (MQ=255)
                                                                                                                              gTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGt                         <  1:6292520/150‑1 (MQ=255)
                                                                                                                              gTACGGTCCCCTCACCCTACCCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGt                         <  2:1504799/150‑1 (MQ=255)
                                                                                                                                 cGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                      <  1:4531541/150‑1 (MQ=255)
                                                                                                                                    tCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGc                   >  2:106156/1‑150 (MQ=255)
                                                                                                                                     ccccTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCg                  <  1:2623903/150‑1 (MQ=255)
                                                                                                                                      cccTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGa                 <  1:9232934/150‑1 (MQ=255)
                                                                                                                                       ccTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCtt                                           <  2:12668794/123‑1 (MQ=255)
                                                                                                                                       ccTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                >  1:1655095/1‑150 (MQ=255)
                                                                                                                                        cTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                                                        >  2:6255940/1‑61 (MQ=255)
                                                                                                                                        cTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt               >  2:8369282/1‑150 (MQ=255)
                                                                                                                                         tCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc              >  1:12118892/1‑150 (MQ=255)
                                                                                                                                           ccccTAACCCTC‑TCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc            <  1:10693386/149‑1 (MQ=255)
                                                                                                                                            cccgcccccTC‑TCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGCATTAACCCGAATCGca           <  2:14126388/144‑1 (MQ=255)
                                                                                                                                              cTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATg                    >  2:6139890/1‑139 (MQ=255)
                                                                                                                                               tAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACaa        <  1:11847519/150‑1 (MQ=255)
                                                                                                                                                  cccTC‑GCCCCTTTGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATg                    >  1:5850332/1‑135 (MQ=255)
                                                                                                                                                  cccTC‑GCCCCTTTGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt     >  2:6906736/1‑150 (MQ=255)
                                                                                                                                                  cccTC‑GCCCCTTTGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt     >  1:5403503/1‑150 (MQ=255)
                                                                                                                                                   ccTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt    <  1:13777326/150‑1 (MQ=255)
                                                                                                                                                    ctc‑tCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtctc                                             >  2:8762629/1‑108 (MQ=255)
                                                                                                                                                    ctc‑tCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACACTGTTc   >  2:4715741/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCT  >  NC_000913/3504724‑3505022

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑