Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,350,8530AC47.9% 145.0 / 22.8 96S725A (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 (41/9);  new base C (0/46);  total (41/55)
Fisher's exact test for biased strand distribution p-value = 1.27e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTT  >  NC_000913/3350707‑3351002
                                                                                                                                                    |                                                                                                                                                     
tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGAcccc                                                                                                                                                   >  2:1283782/1‑151 (MQ=255)
  ccacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCccgc                                                                                                                                                 >  3:1139759/1‑151 (MQ=255)
   cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCccgcc                                                                                                                                                >  2:897622/1‑151 (MQ=255)
   cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCccgcc                                                                                                                                                <  2:559861/151‑1 (MQ=255)
   cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCccgcc                                                                                                                                                >  3:241958/1‑151 (MQ=255)
         cGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCt                                                                                                                                          >  2:1551657/1‑151 (MQ=255)
          gCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCtt                                                                                                                                         >  3:69660/1‑151 (MQ=255)
             ttCAGCACTTCTCCGTCGTGCCCCCCCTCTTCTTCCATCTCTTCAACCCATCCCCCGACGTTATCTCCCCGGGCGGGAGGGTCAGGAGATTGATTTTGCTGCCCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACGGCCCCCGCCGCACCttt                                                                                                                                        <  3:943121/149‑1 (MQ=255)
               cAGCACTTCTCCGTCGTGCCGCCACTCTTCTTTCATCTCTTAAACCCATCCCCCGCCGTTATCTCCCCGGGCGGGAAGGTCAGGAGACTAAATTTGCGGCCCCAGTTGTTGCAGATGGCG‑AAC‑CCCCACTGACCCCGCCGCACCTTTAAt                                                                                                                                     <  3:744200/150‑1 (MQ=255)
               cAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAAtt                                                                                                                                    >  3:1401750/1‑151 (MQ=255)
                aGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAAttt                                                                                                                                   >  2:1097195/1‑151 (MQ=255)
                  cACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                                  >  2:760895/1‑150 (MQ=255)
                     ttCTCCGTCGTGACGCCACTCTTCTTTCATCTCTTAAATCCATTCCCCGACTTTATCTCCCCAGGCGGGAGGGTCGGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCG‑TAC‑CCCCCCTGCCCCCGCCGCACCTTTAAtttt                                                                                                                                  <  2:599027/147‑1 (MQ=255)
                        tccgtcgtGCCCCCACTCTTCTTTCTTCTCTTCAATCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTTTTGCAGATGGCT‑TAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGt                                                                                                                              <  2:11854/146‑1 (MQ=255)
                              gtgccccccctctTCTTTCATCTCTTCAATCCATCCCCCGACGTTATCTCCCCAGGCCGGAGGGTCAGGAGACTGATTTTGCTGACCCAGTTGTTGCAGATGGCT‑TAC‑CCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTca                                                                                                                     <  3:97329/142‑1 (MQ=255)
                               tgccccccctctTTTTCCATCTCTTCAACCCATCCCCCGACGTTATCTCCCCAGGCGGGAAGGTCAGGAGACTGATTTGGCTGCCCCAGCTGTTGCAGAGGGCT‑TAC‑CCCCACGGCCCCCGCCGCACCTTTAATTTTATGtc                                                                                                                             <  2:1154122/134‑1 (MQ=255)
                                gACGCCCCTCTTCTTTAATCTCTTCAATCCTTTCACCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGTTGTTGCAGATGGCG‑TAA‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                                      <  2:398977/148‑1 (MQ=255)
                                 aCGCCACTTTTTTTTCATCTCTTCAACCCATTCCCCGACGTTATCTCCCCGGGCGGGAAGGTCGGGAGACTGATTTTGCTGCCCCAGCTGTTGAAGATGGCG‑TAC‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaac                                                                                                                   <  2:544711/150‑1 (MQ=255)
                                  cGCCCCTCTTCTCCCATCTCTTAAACCCATACACCGACGTTACCTCCCCGGGCCGAAAGGTCAGAAGCCAGAATTTCCTGACCTGGCTGTTGCAGA‑GGCC‑TAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaac                                                                                                                   <  2:747238/149‑1 (MQ=255)
                                    ccccttttcttTCATCTCTTCAACCCTTCCCCCGCCTTTCTCTCCCCGGGCGGGAAGGTCAGGAGACTAATTTTGCGGCCCCAGTTTTTGCAGAGGGCG‑AAC‑CCCCACGGCCCCCGCCGCACCTTTAATTTTATGtcc                                                                                                                            <  3:523172/132‑1 (MQ=255)
                                    cccctctTCTTTCATCTCTTCAATCCATTCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCG‑TAC‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                               <  2:993222/148‑1 (MQ=255)
                                    cccctctTCTTTCATCTCTTCAACCCATTCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCG‑TAC‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                               <  2:551991/148‑1 (MQ=255)
                                    ccACTTTTCTTTCACCTCTAAAATCCTTCCCCCGATTTTACCTCCCCAGGCGGGAGGGTAAGGCGACAAAATTTGTGGCCCCACTTGTTGAAGTTGGCG‑AAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaa                                                                                                                    <  2:951187/146‑1 (MQ=255)
                                     cACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                                              >  3:426952/1‑151 (MQ=255)
                                     cACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                                              >  2:1067709/1‑151 (MQ=255)
                                       ctctTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCt                                                                                                            >  2:1337032/1‑151 (MQ=255)
                                        tcttctTTCATCTCTTCAATCCATACCCCGACGTTATCTCCCCAGGCCGGAGGGTAGGGAGACTGAATTTGCTGCCCCAGCTTTTGCAGATGGCG‑TAC‑CCCCGCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaac                                                                                                                   <  2:904132/143‑1 (MQ=255)
                                        tcttctCTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt                                                                                                           >  2:1226636/1‑151 (MQ=255)
                                           tcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                        >  2:159541/1‑151 (MQ=255)
                                            cttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                       >  2:884791/1‑151 (MQ=255)
                                            cttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                       >  2:1470217/1‑151 (MQ=255)
                                             tttCATCTCTTCACCCCATCCCCCGACTTTATCTCCCCGGGCGGGAAGGTCAGGAGACTAATTTTGCTGCCCCATTTTTTGCAGTGGGCT‑AAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                          <  3:189426/147‑1 (MQ=255)
                                              ttCATCTCTTCAACCCATCCCCCGACGTTATCTCCCCAGGCGGGAAGGTCAGGAGACTAAATTTGCTGCCCCAGCTGTTGCAGATGGCG‑TAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                       <  2:985300/149‑1 (MQ=255)
                                              ttCATCTCTCCAACCCATCCCCCGACTTTATCTCCCCGGGCGGGAGGGTCAGGAGACTGAATTTGCTGCCCCAGTTTTTGCAGATGGCG‑AAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                       <  3:905571/149‑1 (MQ=255)
                                                     ctTCAACCCATCCCCCGACTTTCTCTCCCCAGCCCGGAAGGTCAGGAGACTAATTTTGCTGCCCCATTTTTTACAGATGGCG‑TAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                                 <  2:404171/148‑1 (MQ=255)
                                                     ctTCAACCCATCCCCCGACTTTATCTCCCCAGCCGGGAAGGTCAGGAGATTAATTTTGCTCCCCCAGTTTTTGAAGGTGGCT‑TAC‑CCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcg                                                                                               <  3:899344/150‑1 (MQ=255)
                                                         accccATCCCCCGCCTTTATCTCCCCGGGCCGGAAGGTCGGGAGCCTGATTTTGCTGCCCCAGTTGTTGAAGATGGCT‑TAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                             <  3:1278778/145‑1 (MQ=255)
                                                         accccATCCCCCGACGTTATCTCCCCGGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCG‑TCC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                                            <  3:105043/146‑1 (MQ=255)
                                                         aaTCCATCCCCCGACTTTATCTCCCCAGGCCGGAGGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAC‑CCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcg                                                                                               <  2:70490/146‑1 (MQ=255)
                                                           ccctaccACCGTCTTTATCTCCCCGGGCGGGAAGTTAGGGAGACTAAATTCGCTGCCCCAGTTGTTGCAGGGGGCG‑TAC‑CCCCACTACCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAg                                                                                        <  3:995708/145‑1 (MQ=255)
                                                            ccATTCCCCGACGTTATCTCCCCGGGCGGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGTTGTTGCAGATGGCG‑TAC‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGa                                                                                       <  2:1457658/151‑1 (MQ=255)
                                                            ccATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                          >  3:1404170/1‑148 (MQ=255)
                                                            ccATCCCCCGCCTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTGTTGCGGATGGCG‑TAC‑CCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCa                                                                                         <  2:746913/149‑1 (MQ=255)
                                                             cATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCC‑TTTTTATCCTGCGCCGTCAGAtt                                                                                     >  3:916398/1‑151 (MQ=255)
                                                               tcccccGACTTTATCTCCCCGGGCGGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTTTTGCAGATGGCG‑AAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                                      <  2:1204708/145‑1 (MQ=255)
                                                                 ccccGTCGTTACCTCCCCTGCCCGGAAGGCCAGGAAACTTATTTTGCTGCCCCTGATGTTGAATATGGCG‑TAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                                    <  3:1575530/147‑1 (MQ=255)
                                                                 cACCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGAGGGCG‑TAC‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                                    <  2:518691/149‑1 (MQ=255)
                                                                  aCCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGAc                                                                                 >  2:1612354/1‑151 (MQ=255)
                                                                    cGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                               >  3:235732/1‑151 (MQ=255)
                                                                     gACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                              >  2:275750/1‑151 (MQ=255)
                                                                      actttATCTCCCCGGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCG‑TAA‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCa                                                                             <  2:914533/148‑1 (MQ=255)
                                                                       cGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                              >  2:1020744/1‑149 (MQ=255)
                                                                            tCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGc                                                                       >  3:1008916/1‑151 (MQ=255)
                                                                              ttCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg                                                                     >  3:1524732/1‑151 (MQ=255)
                                                                               tCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg                                                                     >  2:1005102/1‑150 (MQ=255)
                                                                               ccccAGGCGGGAAGGTCAGGAGATTGATTTTGCTCCCCCATTTTTTGAAGTTGCCG‑ACC‑CCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                         <  3:1003135/145‑1 (MQ=255)
                                                                                cccGGGCCGGAAGGTCAGGAGACTAATTTTGCTGCCCCAGTTGTTGCAGATGGCA‑TAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                         <  3:1336903/145‑1 (MQ=255)
                                                                                cccAGGCCGGAGGGTCAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGATGGCG‑TAC‑CCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt                                                                      <  2:845820/148‑1 (MQ=255)
                                                                                  cAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                                 >  3:853956/1‑151 (MQ=255)
                                                                                      ccGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGATGGCG‑TAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAg                                                               <  2:1293753/149‑1 (MQ=255)
                                                                                      ccGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGcc                                                             >  3:670181/1‑151 (MQ=255)
                                                                                          aaGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                         >  2:1442377/1‑151 (MQ=255)
                                                                                           aGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                         >  2:621090/1‑150 (MQ=255)
                                                                                            ggttagatgacaaaaTTTGCTGACCGAGCTATGGCAGA‑GGCGTTAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcat                                                       <  2:812174/138‑1 (MQ=255)
                                                                                            ggTCAGGAGATTGAATTTGCTGCCCCAGCTGTTGCAGAGGGCG‑TAC‑CCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcat                                                       <  3:74801/151‑1 (MQ=255)
                                                                                             gTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatc                                                      >  3:1139712/1‑151 (MQ=255)
                                                                                             gTCAGAAGACTGATTTTGCGGCCCCAGTTGTTGCAGATGGCG‑TAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatc                                                      <  2:1285454/151‑1 (MQ=255)
                                                                                                     aCTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTtctc                                              >  2:1029392/1‑151 (MQ=255)
                                                                                                     aCTAAATTTGCTGCCCCACCTGTTGCGGGTGGCA‑TAA‑CCCCACGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttc                                                <  1:8792/149‑1 (MQ=255)
                                                                                                      cTGATTTTGCTGACCCAGCTGTTGCAGATGGCG‑TAC‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct                                               <  2:812507/149‑1 (MQ=255)
                                                                                                           tctgttcccccAGTTGTTTCTGATGCCA‑TAC‑CCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAACGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACCTAGCCCGGCATCAtct                                                  <  3:657662/133‑1 (MQ=255)
                                                                                                            ttcctcccccAGCTTTTGCAGATGGCG‑TAC‑CCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAacac                                        <  2:1178141/143‑1 (MQ=255)
                                                                                                             tGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACACTGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                                      >  3:966884/1‑151 (MQ=255)
                                                                                                                   cccAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAcc                                >  2:930241/1‑151 (MQ=255)
                                                                                                                      aGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGt                             >  2:303058/1‑151 (MQ=255)
                                                                                                                         tGTTGCAGACGGCC‑CAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg                           <  3:586369/150‑1 (MQ=255)
                                                                                                                          gTTGCAGATGGCG‑TAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTc                            <  3:575738/148‑1 (MQ=255)
                                                                                                                           ttGCAGATGGCG‑TAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGt                         <  2:1473257/150‑1 (MQ=255)
                                                                                                                            tGCAGAGGGCA‑TAC‑CCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg                        <  3:1294346/150‑1 (MQ=255)
                                                                                                                                gATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATc                     >  2:433841/1‑149 (MQ=255)
                                                                                                                                   ggCG‑AAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa                    <  2:1229124/147‑1 (MQ=255)
                                                                                                                                   ggCA‑TAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa                    <  3:303587/147‑1 (MQ=255)
                                                                                                                                       cgtccccccACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTc               <  2:27306/145‑1 (MQ=255)
                                                                                                                                       cgaacccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCt                 <  3:197104/143‑1 (MQ=255)
                                                                                                                                          acccccACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTc           <  2:427807/149‑1 (MQ=255)
                                                                                                                                            ccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg              <  3:271693/146‑1 (MQ=255)
                                                                                                                                            ccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCt         <  3:323969/151‑1 (MQ=255)
                                                                                                                                            ccccACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCt         >  2:608737/1‑151 (MQ=255)
                                                                                                                                             ccccCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTa        <  3:1135230/151‑1 (MQ=255)
                                                                                                                                             cccACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCt         >  2:419233/1‑150 (MQ=255)
                                                                                                                                              ccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAc       <  2:1493577/148‑1 (MQ=255)
                                                                                                                                               cACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATGGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAc       >  3:1520224/1‑150 (MQ=255)
                                                                                                                                               cACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAc       >  3:1542916/1‑150 (MQ=255)
                                                                                                                                                ccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACaa     <  2:1171021/150‑1 (MQ=255)
                                                                                                                                                   gACCCCGCCGCACCTTTAATTTTCTGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGtt  >  2:1562353/1‑151 (MQ=255)
                                                                                                                                                   gACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGtt  >  3:696803/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCG‑TAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTT  >  NC_000913/3350707‑3351002

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