Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,350,8530AC41.5% 176.3 / 13.3 95S725A (TCA→GCA) arcBaerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base A (47/8);  new base C (0/39);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 1.50e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTT  >  NC_000913/3350704‑3351002
                                                                                                                                                       |                                                                                                                                                     
ttttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAc                                                                                                                                                      >  3:19394/1‑151 (MQ=255)
   tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                                   >  3:1320689/1‑151 (MQ=255)
   tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                                   >  1:459535/1‑151 (MQ=255)
      cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                                <  1:1257553/151‑1 (MQ=255)
      cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                                <  1:73937/151‑1 (MQ=255)
         ccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCcgccgc                                                                                                                                              >  2:236200/1‑150 (MQ=255)
           acGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCAcc                                                                                                                                           >  1:1211416/1‑151 (MQ=255)
           acGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCAcc                                                                                                                                           >  3:120402/1‑151 (MQ=255)
                ttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTa                                                                                                                                       >  1:330016/1‑150 (MQ=255)
                      aCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTAt                                                                                                                                >  1:1322933/1‑151 (MQ=255)
                      aCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTAt                                                                                                                                >  3:331776/1‑151 (MQ=255)
                        ttCTCCGTCGTGCCCCCCCCCTTCTTTCTTCTCTTAAACCCTTCCCCCGATGTTTTTTCCCCAGGCCGGAAGGTCAGGAAACTTAATTTGCTCCCCTTGTTTTTGAAGTGGGGGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTAt                                                                                                                                <  4:77337/149‑1 (MQ=255)
                              gtcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcc                                                                                                                        >  3:616430/1‑151 (MQ=255)
                               tcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                                                       >  3:1075356/1‑151 (MQ=255)
                               tcgtGACCCCCCTCTTCTTTCATCTCTCCAACCCATCCACCGCTGTTATCTCCCCAGGCCGGAAGGTCAGGGAACTGATTTTTCTGACCCAGCTGTCGGAGTTGGCGAAACCCCTCGGTCCCCGCCGCACCTTTAATTTTATGtccttcc                                                                                                                        <  3:1317533/150‑1 (MQ=255)
                                  tGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTca                                                                                                                     >  1:797392/1‑150 (MQ=255)
                                  tGACCCCCCCCTTCTTCCACCTCTAAAACCCATCCCCCGACTTTTTCTCCCCGGCCCGGAGGTGAAGGGCACTAATTTTGCTCACCCAGCTTTTGCAGAGGGCGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGtccttc                                                                                                                         <  1:93368/146‑1 (MQ=255)
                                     cccccctctTTTTTCATCTCTTCAACCCATCCCCCGCCTTTATCTCCCCGGGCGGGAAGGTCAGGAGACTGAATTTGCTGACCCAGTTGTTGCAGAGGGCGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca                                                                                                                  <  3:510288/145‑1 (MQ=255)
                                     cccccctctTCTTCCATCTCTTAAACCCATCCCCCGACTTTATCTCCCCAGGCGGGAAGGTCGGGAGATTAATTTTCCTGCCCTATTTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaac                                                                                                                   <  4:301104/144‑1 (MQ=255)
                                      cccccttttcttcCATCTCTTCAACCCATCCCCCGACGTTATCTCCCAGGGCGGGAGGGTCAGGAGACTGAATTTGCTGACCCAGTTGTTGCAGAGGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAAt                                                                                                                <  3:522050/144‑1 (MQ=255)
                                       cccctctTTTTTCATCTCTTCAATCCATCCCCCGCCGTTATCTTCCCAGGCCGGAGGGTCAGGAGACTAATTTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                               <  3:187218/148‑1 (MQ=255)
                                       ccACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                               >  1:151567/1‑151 (MQ=255)
                                            cttcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                          >  4:358009/1‑151 (MQ=255)
                                              tcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                        >  4:111663/1‑151 (MQ=255)
                                                tttCATCTCTTCAATCCATTCCCCGACGTTATCTTCCCAGGCCGGAGGGTCAGGAGACTGAATTTGCTCACCCAGCTGTTGCAGAGGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt                                                                                                           <  2:547233/146‑1 (MQ=255)
                                                tttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                       >  4:104228/1‑150 (MQ=255)
                                                tttCATCTCTTCAACCCATTCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                          <  3:568499/147‑1 (MQ=255)
                                                tttCATCTCTACATGCCTTCCCCCGACTTTATCTTCCCGGCCCGGAGGGTTGGGAGGCTGATTGGGCTGACCCAGTTGTTGCAGAGGGGGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                        <  1:499069/149‑1 (MQ=255)
                                                 ttctcctgttCAACCCATCCCCCATCTTTATCTTCCAAGCCAGGGAGGTCGGGAAACAGATTTTTCTGCCCCAGTTGTTGCAGATGGCGAAACCCCCTTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTa                                                                                                      <  3:1378747/142‑1 (MQ=255)
                                                   cATCTCTTCAACCCATCCACCGACTTTATCTCCCCAGGCCGGAAGGCCAGGAGACTAATTTGGCTACCCCATTTGTTGCAGTGGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                        <  1:904760/146‑1 (MQ=255)
                                                      ctctCCAATCCATCCCCCGCCGTTATCTCCCCAGGCGGGAAGGTCAGGAGATTGAATTTGCTGACCCGGCTGTTGCAGATGGCGTACCCCCACTAACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATc                                                                                                    <  2:540767/147‑1 (MQ=255)
                                                      ctctCCAACCCATCCCCCGACTTTATCTCCCCAGGCGGGAGGGCCAGGAGACTGATTTGGCTCCCCCATCTTTTGAAGGGGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATcc                                                                                                   <  1:166487/148‑1 (MQ=255)
                                                      ctcctcAACCCATCCGCCCACTTTTTCTCCCCAGGCCGGAAGGTCAGAAGAATGATTTTGCTAACCCTGCTGTTCCAAATGGCAAACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt                                                                                                           <  3:1280970/140‑1 (MQ=255)
                                                        ctTCAACCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                                 <  3:1202578/148‑1 (MQ=255)
                                                         ttCAATCCATCCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTCACCCAGTTGTTGCAGATGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                                <  4:540816/148‑1 (MQ=255)
                                                          tCAATCCATTCACCGACGTTATGTTCCAAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcg                                                                                               >  4:245116/1‑148 (MQ=255)
                                                           cAACCCATCCCCCGACGTTTTCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                                            <  2:93061/150‑1 (MQ=255)
                                                            aacccATCCACCGACGTTATCTCCCCAGGCGGGAAGGCCAGGAGACTAATTTTGCTGCCCCGGCTGTTGCAGAGGGCTTACCCCCACTGCCCCCGCCGCACCTTAAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                                            <  3:505359/146‑1 (MQ=255)
                                                            aaTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGt                                                                                           >  1:168163/1‑150 (MQ=255)
                                                              ccaatcccccGACGTTATCTTCCCAGGCCGGAGGGTCAGGAGACTAAATTTGCTGCCCCAGTTGTTGCAGAGGGCGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCa                                                                                         <  2:515250/142‑1 (MQ=255)
                                                               ccATTCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGACCCAGTTGTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCa                                                                                         <  3:330052/149‑1 (MQ=255)
                                                                 atcccccGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAg                                                                                        <  3:452526/143‑1 (MQ=255)
                                                                  ttCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                                    >  1:849713/1‑151 (MQ=255)
                                                                    ccccgcctttATCTCCCCGGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTGTTGAAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                   <  1:778167/142‑1 (MQ=255)
                                                                    ccccGACTTTATCTCCCCGGGCCGGAAGGTCGGGAGACTGAATTTGCTGCCCCATTTTTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                                            <  1:570402/139‑1 (MQ=255)
                                                                       cGACTTTTTCTTCCCAGGCCGGGAGGGCAGGGGCTTGATTTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                               <  1:159855/151‑1 (MQ=255)
                                                                       cGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                               >  4:493184/1‑151 (MQ=255)
                                                                          cGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                              >  1:1344145/1‑149 (MQ=255)
                                                                           gTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                                           >  4:244342/1‑151 (MQ=255)
                                                                            ttATCTCCCCAGGCCGGAGGGTCGGGAGACTGATTTTGCTGCCCCAGTTTTTGCAGATGCCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAg                                                                            <  4:311120/149‑1 (MQ=255)
                                                                             tATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                         >  1:1352844/1‑151 (MQ=255)
                                                                             tATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                         >  1:622524/1‑151 (MQ=255)
                                                                               tCTTCCCAGGCCGGAAGGTCGGGAGACTAAATTTGCTCACCCATCTTTTGCAGATGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCa                                                                          <  1:705765/148‑1 (MQ=255)
                                                                               acTTCCCGAGGGCAGAGGGTCAGGAGACTTATTGTGCTGCCCGAGCTTTTGCAGACGGCTTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAACGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                         <  1:477399/149‑1 (MQ=255)
                                                                                cTTCCCAGGCCGGAAGGTCAGGAGACTAAATTTGCTCCCCCATTTTTTGCAGGTGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                         <  1:1027873/148‑1 (MQ=255)
                                                                                 tccccAGGCCGAAAGGTCAGGAGACTAATTTGCCTGCCCCAGTTGTTGCAGAGGGCTTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCa                                                                          <  1:835531/144‑1 (MQ=255)
                                                                                      cacgCCGGAAGGGCGGGAGATTGATTCTGCTGACCCACCTGTTGCAGATGCCATCCCCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGa                                                                    <  1:788033/145‑1 (MQ=255)
                                                                                       aGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                                 >  4:208486/1‑150 (MQ=255)
                                                                                       aGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                                 >  4:234492/1‑150 (MQ=255)
                                                                                        ggCCGGAGGGTCAGGAGACTGAATTTCCTGCCCCAGTTGTTGAAGATGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATa                                                                <  1:389558/150‑1 (MQ=255)
                                                                                        ggCCGGAAGGCCGGGAGACTAAATTTGCTCACCCAGCTTTTGCAGAGGGCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                                 <  1:704557/149‑1 (MQ=255)
                                                                                           cGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                            >  2:46979/1‑151 (MQ=255)
                                                                                           cGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                            >  3:349060/1‑151 (MQ=255)
                                                                                               aGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                         >  3:1262274/1‑150 (MQ=255)
                                                                                                ggccAGGAGACTGAATTTGCTGACCCAGCTGTGGCAGTGGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcat                                                       <  1:1011415/148‑1 (MQ=255)
                                                                                                  tCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatc                                                      >  1:1354504/1‑150 (MQ=255)
                                                                                                    aGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCatcatc                                                   >  3:1141993/1‑151 (MQ=255)
                                                                                                     ggactcgcattttCCTGCCCCAGCTGTTGCAGATGCCTAACCCCCCCTTCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat                                                    <  1:677700/139‑1 (MQ=255)
                                                                                                     ggAGCCTAATTTTGCGGCCCCATTTGTTGCAGAGGGCGTACCCCCCCCAACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtct                                                  <  4:380437/151‑1 (MQ=255)
                                                                                                     ggAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtct                                                  >  3:1363829/1‑151 (MQ=255)
                                                                                                      gagaCTGATTTTCCTGCCCCATCTTTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt                                                 <  4:401302/151‑1 (MQ=255)
                                                                                                       agaCTGATTTTCCGCCCCTAGCTGTTGCAGTTGGCTTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttc                                                <  1:858768/151‑1 (MQ=255)
                                                                                                          cTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTtctc                                              >  4:421863/1‑150 (MQ=255)
                                                                                                            gAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa                                            >  1:317902/1‑150 (MQ=255)
                                                                                                             aaTTTGCTGACCCAGCTGTTACAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTGTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCAAGAACGCTGACATAGCACGGCATCATCTTCTCAAac                                          >  4:13340/1‑151 (MQ=255)
                                                                                                               tttGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAaca                                         <  3:1220010/150‑1 (MQ=255)
                                                                                                               attGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAaca                                         >  1:928896/2‑150 (MQ=255)
                                                                                                                     tGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt                                    >  3:1413929/1‑150 (MQ=255)
                                                                                                                      gACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTa                                   >  1:1387840/1‑150 (MQ=255)
                                                                                                                       aCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTaa                                  <  3:667367/150‑1 (MQ=255)
                                                                                                                        cccAGTTTTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt                                    <  3:832629/147‑1 (MQ=255)
                                                                                                                        cccAGCTGTTGCAGAGGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt                                    <  3:624368/147‑1 (MQ=255)
                                                                                                                         ccAGTTTTTGAAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAccc                               <  2:100346/151‑1 (MQ=255)
                                                                                                                         ccAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAccc                               >  3:529021/1‑151 (MQ=255)
                                                                                                                            gCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTc                            >  1:564977/1‑151 (MQ=255)
                                                                                                                                 tGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg                        >  1:1314598/1‑150 (MQ=255)
                                                                                                                                   cAGATGGCGTAACCCCACTAACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGAt                      <  2:432433/150‑1 (MQ=255)
                                                                                                                                          cGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGAGATCAGCTTCg              >  4:380451/1‑151 (MQ=255)
                                                                                                                                            tAACCCCACTCACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACCATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATCGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCgg             >  4:506055/1‑150 (MQ=255)
                                                                                                                                             acccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGt            <  2:437047/148‑1 (MQ=255)
                                                                                                                                               ccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCgg             <  1:616345/147‑1 (MQ=255)
                                                                                                                                                cccACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCt         >  3:149567/1‑150 (MQ=255)
                                                                                                                                                 ccACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAc       >  4:442923/1‑151 (MQ=255)
                                                                                                                                                      gACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACa      >  4:29870/1‑147 (MQ=255)
                                                                                                                                                      gACCCCACCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGtt  >  1:769163/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTT  >  NC_000913/3350704‑3351002

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