Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,350,8530AC41.8% 163.4 / 11.8 91S725A (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 (42/11);  new base C (0/38);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 5.85e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTT  >  NC_000913/3350704‑3351002
                                                                                                                                                       |                                                                                                                                                     
ttttGACACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGa                                                                                                                                                       >  1:140784/1‑150 (MQ=255)
   tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                                   >  2:394332/1‑151 (MQ=255)
    gCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                                   <  2:778904/150‑1 (MQ=255)
     ccacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgc                                                                                                                                                 >  2:587016/1‑151 (MQ=255)
      cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                                >  1:938956/1‑151 (MQ=255)
      cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                                <  1:314768/151‑1 (MQ=255)
      cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                                >  1:505080/1‑151 (MQ=255)
              cTTCCAGCACTTCTCCGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATCCCCCGCCGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGTTGTTGCAGATGGCGTACCCCCCCTGACCCCGCCGCACCtt                                                                                                                                         <  1:12133/150‑1 (MQ=255)
               tttCAGCACTTCTACGTCGGGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                                        >  1:649510/1‑150 (MQ=255)
                  cAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtt                                                                                                                                    >  2:591158/1‑151 (MQ=255)
                  cAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtt                                                                                                                                    >  3:512591/1‑151 (MQ=255)
                     cccTTCTCCGTCGTCACGCCCCTTTTCTTTCATCTCTTCAACCCATCCCCCGACGTTTTCTTCCCAGGCCGGAAGGTCAGGACACTAATTTTGCTGACCCGGCTGTGGCAGTTGGTGTACCCCCACAGCCCCCGCCGCACCTTTa                                                                                                                                       <  1:246180/143‑1 (MQ=255)
                     catttCTCCGTCGTCACCCCCCTTTTCTTCCATCTCTCCAACCCATCCCCCGACTTTATCTTCCAAGGCGGGAGGGTCAGGAGACGAAATTGCCTGCCCCATTTTTTGAAGAGGGCGAAACCCCACGCACCCCGCCGCACCTTTAATTTTa                                                                                                                                 <  4:297437/148‑1 (MQ=255)
                     cACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                                  >  2:978019/1‑150 (MQ=255)
                             cgtcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttc                                                                                                                         >  1:76435/1‑151 (MQ=255)
                                cgggACCCCCCTCTTTTTTCATCTCTTCAACCCATCCCCCGCCTTTATCTCCCCGGGCCGGAAGGTCAGGAGACTAAATTTGCTGCCCAGTTTTTTGAAGAGGGCGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTa                                                                                                                                 <  1:736517/137‑1 (MQ=255)
                                   gccccccctctTCTTTCATCTCTTCAATCCATCCCCCGCCTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGACCCAGCTTTTGCAGATGGCGAAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTca                                                                                                                     <  4:346051/142‑1 (MQ=255)
                                     cGCCACTCTTCTTTCATCTCTTAAATCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTAAATTTGCTGACCCAGCTTTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                                                       <  2:766811/145‑1 (MQ=255)
                                       cccctctTCTTTCATCTCTTCAACCCATCCCCCGACGTTATCTCCCAAGGCGGGAAGGTGAGGAGACTAAATTGGCTGACCAAGCTGTTGAAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAAt                                                                                                                <  2:404967/147‑1 (MQ=255)
                                        ccctctTCTTTCATCTCTTCAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCGGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                               <  1:1521771/148‑1 (MQ=255)
                                        ccctctTCTTCCATCCCTTCAACCCATCCCCCGATGTTATCTTCCCGGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGTTTTTGCAGATGGCGTACCCCCCCGGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                                              <  1:586933/149‑1 (MQ=255)
                                         aCTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                                                                                             >  1:1086914/1‑151 (MQ=255)
                                           tcttctTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt                                                                                                           >  1:423051/1‑151 (MQ=255)
                                            cttcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                          >  2:1560776/1‑151 (MQ=255)
                                             ttttttCACCTCTTCAATCCATTCCCCGACTTTATCTCCCCAGGCGGGAGGGTGAGGAGCCTGATTTTCCTGACCCAGCTGTTGAAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                          <  2:299076/147‑1 (MQ=255)
                                              tcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                        >  1:1702106/1‑151 (MQ=255)
                                                tttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                       >  2:1063980/1‑150 (MQ=255)
                                                tttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                       >  1:802857/1‑150 (MQ=255)
                                                ttcattctctcCAACCCCTCCCCCGAATTTTTCTTCCCGGGGGGGAGGGACAGGAGATTGAATTGCCGGACCCGTGTTTTGCAGGGGGCCCACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                          <  1:897257/142‑1 (MQ=255)
                                                 ttCATCTCTTCAACCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                        <  2:717469/148‑1 (MQ=255)
                                                     tctcTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCAGCACCTTTAATTTTATGTCCTTCCTCAACAATGCATTTTTTATCCt                                                                                                  >  1:67142/1‑150 (MQ=255)
                                                      ctctTCAATCAATCCCCCGACTTTATCTTCCCGGGCGGGAAGGTCAGGAGAATGATTTTGCTGCCCTAGTTTTTGCAGAGGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                                                     <  1:492005/146‑1 (MQ=255)
                                                        ctacAACCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTAAATTTGCTGACCCAGCTTTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                                 <  4:48330/145‑1 (MQ=255)
                                                        ccacaacccctcccccGACGTTTTTTCCCCAGCCGGGAGGGCAAGAACAGTAAATTTGCTCACCAAGTTTTTGGAGAGGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                                 <  1:1266247/134‑1 (MQ=255)
                                                         ttCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                             >  2:1619601/1‑151 (MQ=255)
                                                            aaTCCATTCCCCCACGTTATCTCCCCAGGCCGGAAGGTAAGGAAACTAAATTGCCTCCCCATGTTTTTGCGGGTGGCTAACCCCCTCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                                            <  1:658659/149‑1 (MQ=255)
                                                            aaTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                          >  2:967090/1‑151 (MQ=255)
                                                               ccATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                          >  4:16383/1‑148 (MQ=255)
                                                                ccccccccctttttATCTCCCCAGGCGGGAAGCTCAGGAGACTGATTTTCCAAACCCAGTTGTTGATGAGGGAGAACCCCCATTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                                      <  1:581636/139‑1 (MQ=255)
                                                                 atcccccGACGTTATCTTCCCGGCCGGGAAGGTAAGGAGACTGATTTTGCTGACCCAGTTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                                      <  2:215129/145‑1 (MQ=255)
                                                                 atccACCGATGTTATCTCCCAAGGCCGAAGGATAAGATGATGAAATTGGCTGACCCAGCTTTTGAAGATGGCTTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCa                                                                                         <  1:428699/144‑1 (MQ=255)
                                                                    cACCGACTTTTTCTCCCCAGGCCGAAGGGCCAGGAGACTAATTTTGTTGACCCAGCTTTTGAAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGa                                                                                  <  1:115400/151‑1 (MQ=255)
                                                                       cGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                               >  2:1296092/1‑151 (MQ=255)
                                                                        gACGTTATCTCCCCGGGCGGGAAGGTCAGAAGCCTGATTGTGCTCCCCCAGTGTTTGAAGAGGGCGTACCCCCTCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                              <  1:1605813/151‑1 (MQ=255)
                                                                         actttATCTCCCGAGGCTGGAGGCCAGGGACACTGATTTTGCTGCCCCCGCTGTGGCGGAGGGGTTACCCCCTCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACt                                                                                <  2:109326/145‑1 (MQ=255)
                                                                         aCGTTATCTCCCCAGGCCGGAAGGTAAGGAGACTAATTTGGCTGACCCAGCTGTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                              <  2:270793/150‑1 (MQ=255)
                                                                            ttATCTCCCCAGGCGGGAAGGTCAGGAGACTGAATTTGCTGCCCCTTTTGTTGAAGTGGGCTAAACCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCa                                                                             <  1:561509/148‑1 (MQ=255)
                                                                             tATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                         >  1:737210/1‑151 (MQ=255)
                                                                                cTTCCCATTCCGGAAGGACAGCAGACGCAATTTGCTGACACAGCTGTTGCAGATGGCGTAACCCCGCTGACCCCCCCGCACCTTTAATTTTATGTCCTTCCTCATCAATGCCTTTTTTATCCTGCGCAGTCAGATTCCACTCGTGCAc                                                                         >  1:312885/1‑148 (MQ=255)
                                                                                   cccGGGCGGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                                                   <  2:1529468/151‑1 (MQ=255)
                                                                                     cAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                                 >  2:1283233/1‑151 (MQ=255)
                                                                                        gCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGc                                                              >  2:739315/1‑151 (MQ=255)
                                                                                         ccGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGcc                                                             >  2:256185/1‑151 (MQ=255)
                                                                                          cGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                            >  2:167014/1‑151 (MQ=255)
                                                                                           ggAGGGTCAGGGGACTAAATTTGCTGCCCCTGCTGTTGCAGTTGCCTTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                            <  1:102991/150‑1 (MQ=255)
                                                                                           ggAGGGTCAGGAGACTGAATTTGCTGACCCAGTTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                            <  1:1520339/150‑1 (MQ=255)
                                                                                           ggAAGGTCGGGAGCCTGATTTTGCGCCCCCGGCTGTTGCGGAGGGCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAg                                                               <  1:354922/147‑1 (MQ=255)
                                                                                           ggAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCg                                                           >  2:724567/1‑151 (MQ=255)
                                                                                             aaGGTCAGGAGACTAAATTTGCTGCCCCAGCTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCgg                                                          <  3:197716/150‑1 (MQ=255)
                                                                                             aaGGCCAGGAGACAAATTTTGCTGCCCCAGCTGTTGCAGAGGGCTTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                         <  1:776180/151‑1 (MQ=255)
                                                                                               ggTCGGGAGAGGGATTTTGCTAACCCCGCTTTTGCAGGGGTCGGACCCCCCGGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGca                                                        <  2:1271681/150‑1 (MQ=255)
                                                                                               ggTCAGGAGACTGATTTGGCTGACCCAGTTTTTGCAGATGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcat                                                       <  1:1263094/151‑1 (MQ=255)
                                                                                                     gaagaCTAATTTTGCTCACCCAGTTTTTGAAGAGGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatc                                                      <  1:1000783/145‑1 (MQ=255)
                                                                                                       agaCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttc                                                >  3:435692/1‑151 (MQ=255)
                                                                                                       agaCTAATTTTGCTGACCTGGTTTTTGCAGATGGCATACCCCCCCTACCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat                                                    <  1:617316/147‑1 (MQ=255)
                                                                                                            gAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa                                            >  2:383458/1‑150 (MQ=255)
                                                                                                             aaTTTGCTCACCCAGTTTTTGCAGGTGGCTCACCCCCCCCCACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaaa                                           <  4:480764/150‑1 (MQ=255)
                                                                                                               ttggCTGCCCCAGGTTTTGCAGAGGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCa                                             <  1:148976/143‑1 (MQ=255)
                                                                                                                ttGCTGACCCAGCTTTTTCAGAGGGCGTAACCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCGCGCTGACATAGCCCGGCATCATCTTCTCAAaca                                         <  2:667508/149‑1 (MQ=255)
                                                                                                                 tGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                                      >  2:931121/1‑151 (MQ=255)
                                                                                                                    tGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt                                    >  4:9397/1‑150 (MQ=255)
                                                                                                                      gccccgtcTGTTGCAGAGGGCATACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                                      <  1:1657232/140‑1 (MQ=255)
                                                                                                                      ccccAGCTGTTGCAGATGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAc                                 <  3:57859/150‑1 (MQ=255)
                                                                                                                       cccAGTTGTTGCAGATGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTa                                   <  4:496588/148‑1 (MQ=255)
                                                                                                                              tGTTGCAGATGGCTAACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGt                             <  1:189925/148‑1 (MQ=255)
                                                                                                                              tGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg                          >  1:1063129/1‑151 (MQ=255)
                                                                                                                              tGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTCCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg                          >  2:418558/1‑151 (MQ=255)
                                                                                                                                ttGCAGGTGACGAACCCCCCCGCACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCGTCTCAAACGCCGCTAACCCGTCgg                          <  1:1412392/149‑1 (MQ=255)
                                                                                                                                ttGAAGATGGCTTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGt                         <  2:1489612/150‑1 (MQ=255)
                                                                                                                                 tgaggcgggCGTACCCCCCCCGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg                        <  1:427038/143‑1 (MQ=255)
                                                                                                                                 tGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg                        >  1:665066/1‑150 (MQ=255)
                                                                                                                                    aGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATc                     >  2:548854/1‑150 (MQ=255)
                                                                                                                                     gATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATc                     >  3:309372/1‑149 (MQ=255)
                                                                                                                                       tGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCt                 >  1:281682/1‑151 (MQ=255)
                                                                                                                                         gctaacccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTc               <  1:433763/145‑1 (MQ=255)
                                                                                                                                         gCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTc               >  1:1052336/1‑151 (MQ=255)
                                                                                                                                         gCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGGCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTc               >  2:1569102/1‑151 (MQ=255)
                                                                                                                                          cGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg              <  2:823564/151‑1 (MQ=255)
                                                                                                                                               ccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCt         <  2:161764/151‑1 (MQ=255)
                                                                                                                                                  cACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAc       >  2:1599481/1‑150 (MQ=255)
                                                                                                                                                   aCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACa      >  2:1124674/1‑150 (MQ=255)
                                                                                                                                                      gACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGtt  >  1:1661698/1‑151 (MQ=255)
                                                                                                                                                       aCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGtt  >  1:510848/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTT  >  NC_000913/3350704‑3351002

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