Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,665,7020CA54.4% 27.9 / 103.9 79A37S (GCG→TCG) gadXacid resistance regulon transcriptional activator; autoactivator
Reads supporting (aligned to +/- strand):  ref base C (20/16);  new base A (0/43);  total (20/59)
Fisher's exact test for biased strand distribution p-value = 2.76e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-07

GTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTAC  >  NC_000913/3665553‑3665840
                                                                                                                                                           |                                                                                                                                              
gTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCa                                                                                                                                               >  1:1009744/1‑151 (MQ=255)
      aaCGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAcc                                                                                                                                        <  1:542632/151‑1 (MQ=255)
               tGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCg                                                                                                                            <  2:182516/151‑1 (MQ=255)
                 aagaaACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGc                                                                                                                           <  2:782318/150‑1 (MQ=255)
                       cTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaa                                                                                                                    >  1:37762/1‑151 (MQ=255)
                           tttATGGGAAGGAAGATTTAATGCCCCCTCCTTGAGCATGGGAAAAAAAAATAAAAGCGTGCCCCTTAACACAAATAAAAAATTTTCAACACATTTATCTACTCAAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAg                                                                                                                <  2:562492/151‑1 (MQ=255)
                                         gATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                                                                                  >  2:390446/1‑151 (MQ=255)
                                           tttAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTaa                                                                                                >  2:73100/1‑151 (MQ=255)
                                           tttAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTaa                                                                                                >  1:92645/1‑151 (MQ=255)
                                             taacccctcctcctTGAGCATCGGAAAAAATAATGAAGGCGTTCCCCTTGACACAAAAAAAAAATTTTAACCACACTTATCTACTAAAATTTGCCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                              <  1:366060/146‑1 (MQ=255)
                                             taacccccccccctTGACCATGGGGAAAAAAAATGAAGGCGGTTCCCTTGACACGAAAAAAAAATTTACAACCCACTTATCTACTCGAATTTGGCTTCCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                              <  2:85666/140‑1 (MQ=255)
                                             tAATGCCTCCTCCGTGAGCATCGGGAAAAATAATAAACGCGTGTCCCTTGACACAAAAAAAAAATTTTCAACACCCTTATCTTCTCAATTTTGGTTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                              <  1:425347/151‑1 (MQ=255)
                                             tAATCCCCCCTCCTTGAGCAGCGGGAAAAATAAGAAAAGCGTCTCCTTTGACACGAATAAAAAATTTTCAACACACTTACCTACAAAAATTGGGCTTGCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                              <  1:416985/151‑1 (MQ=255)
                                             tAATCCCCCCCCCTTGAACATCGGGAAAAATAATGAAAGCGTTCCCCTTGACCAAAATACAAAATTTACAACACCCTTATTTCCCAAATTTTGTTTTCCACCCAAAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                              <  1:28022/151‑1 (MQ=255)
                                                 ccctcctcctGGAGCATCGGGAAAAATAATGAAAGTTTGTCCTTTGACAAAAATACAAAATTTCCAACCCACTTTCCTCCTGAATTTTGCTTTCCACCAAAAAAACCAGGTCCACCGCCATTAAAATAGCGATATTCACCATTAACCATgg                                                                                          <  1:567244/150‑1 (MQ=255)
                                                   cccctcctgaagatacGGGAAAAAAAATAAAGGGGTTCCCTTTGCCACAAATAAAAATTTTCCACCCCCCTTTTCTCCTCAATTTTGGTTTCCACCCCAAAAAACAGGTTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                                                                        <  2:615292/136‑1 (MQ=255)
                                                    tcctcctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACCCACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGCTATTCACCATTAACCATGGTga                                                                                       >  2:207218/1‑151 (MQ=255)
                                                        cctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat                                                                                   >  1:60419/1‑151 (MQ=255)
                                                        cctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat                                                                                   >  1:999054/1‑151 (MQ=255)
                                                        cctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat                                                                                   >  1:15385/1‑151 (MQ=255)
                                                           tGAGCATCGGGAAAAAAAATAAAAGCGTGTCCCTTGACACAAATAAAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCACCCGAAAAAACAAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatata                                                                                <  2:530027/151‑1 (MQ=255)
                                                               tcatgggAAAAATTATTAAAGGCTTCCCCTTGACCAGAATGCCAAAATTTCCAACACATTTTTCTGTTGAATTTGGCTTCACCCCGCCAAAAACGTGGTCTCCGCCATTAAAATAGCGATATGCTCCATTAGCCATGGTGAGAATATAtt                                                                              <  1:144336/146‑1 (MQ=255)
                                                               cagggggAAAAATAATGAAAGCGTTTCCCTTGACACAAAAAAAAAATTTTCAACCCCCTTATTTACTAGATTTTGGTTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                              <  2:223265/145‑1 (MQ=255)
                                                               cagcGGGAAAAATAATAAAGGCGTGTCCCTTGCAAAGAATAAAAAATTTACAACACCCTTATCTACTCGAATTTGGCTTGCACCCGAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                              <  1:622815/146‑1 (MQ=255)
                                                               cATCGGGAAAAATAATGAAAACGTGGCCCTGGACAAAAAAAAAAAATTAACAACCCTATTAAAAACTCGAATTTGGCTTACATCCAAAAAACCCAGTTTTCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAttt                                                                             <  2:926202/150‑1 (MQ=255)
                                                               cATCGGGAAAAAAAATGAAAGCGTCCCCCTTGACACGAAAAAAAAATTTTCAACACACTTATCTCCTCAAATTTGGCTTGCACCCAAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                              <  2:71648/149‑1 (MQ=255)
                                                                  cGGGAAAAAAAATAAAGGCTTGTCCCTTGACACAAATACAAAATTTCCAACCCATTTATCTCCTCGATTTTTGCTTGCATCCCAAAAAACCAGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                         <  2:521527/151‑1 (MQ=255)
                                                                  cGGGAAAAAAAAAAAAAGCGTTTCCCTTGACACAAAAAAAAATTTTTCAACCCCCTTATCTACCAAATTTTGGTTTCCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                         <  2:947097/151‑1 (MQ=255)
                                                                   gggAAAAATAATGAAAGCGTGTCCCTTAACACGAAAAAAAAATTTACAACACATTTATCAACTGGATTTTGGTTTTCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                         <  2:1032754/150‑1 (MQ=255)
                                                                   gggAAAAATAATAAAATCGTTCCCTTTGACACAAATAAAAATTTTTCACCACACTTTTCTCCTCGAATTTGTTTTGCATCCACAAAAACCAGGTCTCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                         <  2:118205/150‑1 (MQ=255)
                                                                   gggAAAAAAAATAAAGGCGTGCCCTTTGCCAAAAATAAAAATTTTTAAACACATTTTTCCACGAAAATTGGGCTTCCATCCGAAAAAACCGGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                         <  2:1010103/150‑1 (MQ=255)
                                                                   ggaaaaaaaaaTGAAAGCGTCTCCCTTAAAACAAAAAAAAAATTTTCACCCCATTTATCTCCCCAAATTTGGTTTGCACCCGAAAAAACAAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                         <  2:1053841/145‑1 (MQ=255)
                                                                     gaaaaaTAATAAAGGCGTTTCCTTTGACAAAAATACAAAATTTCAAACACCTTTATCTACTCGATTTTGGCTTCCACCCGAAAAAACCAGTTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCt                                                                       <  2:545355/147‑1 (MQ=255)
                                                                          ataatGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATa                                                                 >  2:1022608/1‑151 (MQ=255)
                                                                                      ttccctttGACACAAAAAAAAAATTTCCACCCCTATTACCTACCCAATTTTGGCTTGCACCCCCAAAAACCGGGTCTCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                      <  2:163671/144‑1 (MQ=255)
                                                                                      ttccccTTGACAAAAAAAAAAATTTTACACCACACTTATCTCCTCAATTTTGGTTTGCACCCGAAAAAACCGGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                      <  2:654755/147‑1 (MQ=255)
                                                                                      tgtCCTTTGACACAAATACAAAATTTACAACACACTTATCTACTAGATTTTGCCTTGCATCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                      <  1:67383/150‑1 (MQ=255)
                                                                                      tgtCCTTTGACAAAAAAAAAAAATTTCCACCACACTTATCTACACAATTTTGGCTTGCACCCAAAAAAACCGGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                      <  1:205086/150‑1 (MQ=255)
                                                                                      tcccccTTGACACAAAAAAAAAATTTCCACCACCCTTATCTACCCGATTTTGGCTTGCCCCCAAAAAAACCGGGCCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                      <  2:37768/147‑1 (MQ=255)
                                                                                       gtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACaa                                                    >  1:907726/1‑151 (MQ=255)
                                                                                         cccTTAACAGAAATCAAAAATTTACAACCCATTTAGTTACTCAAATTTGGCTTTAACCCAAAAAAACCAGGTATCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTATGTCTTGCATACGCAATGAGACAATTc                                                 <  2:486796/151‑1 (MQ=255)
                                                                                           cTTGACACAACAAAAAATTTTTCAACTCTTTTCTCTACTAGAATTTGGTTTCCACCCGAAAAACCGAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTcc                                                <  2:966686/151‑1 (MQ=255)
                                                                                                cacGAATAAAAAATTTCCAACACACTTATCTACTCAAATTTGGCTTGCACCCAAAAAACCAAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                           <  1:858047/151‑1 (MQ=255)
                                                                                                 cacaaaaaaaaaaTTTTCACCACATTTATCTACTCAAATTGGGTTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                           <  2:311430/142‑1 (MQ=255)
                                                                                                  cacaaacaaaaaTTTTCCAACACTTTTATCTACTCGAATTGGGCTTCCCTCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCCCCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                           <  2:11779/143‑1 (MQ=255)
                                                                                                  caaaTAAAAATTTTTCAACACACTTATCTACTCAATTTTGCCTTGCATCCACAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAg                                         <  1:81292/149‑1 (MQ=255)
                                                                                                  caaaTAAAAATTTTTAAACCCATTTATCTACTCAAATTTGGCTTGCACCCAAAAAACCAGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAg                                         <  2:842965/149‑1 (MQ=255)
                                                                                                  caaaTAAAAAATTTTCAACACCCTTATCTACTCAATTTTGGTTTCCATCCGCAAAAACCAGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAg                                         <  1:171051/149‑1 (MQ=255)
                                                                                                   aaaTAAAAAATTAACAACACATTTACCATCTGAATTTTCTGTTGCATCCGAAAAAACCGTAACACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATg                                            <  1:962528/146‑1 (MQ=255)
                                                                                                             aTTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTc                                                                                                                                >  2:105587/1‑54 (MQ=255)
                                                                                                             aTTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTc                                                                                                                                <  1:105575/54‑1 (MQ=255)
                                                                                                                ttccAACCCACTTTTCTACTCAAATTTGTTTTGCACCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGt                             <  1:390311/148‑1 (MQ=255)
                                                                                                                ttCAACCACACTTATCTCCCCAATTTTGTCTTCCACCCAAAAAAACCAGGCAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGt                             <  1:291696/151‑1 (MQ=255)
                                                                                                                 tCAACACATTTATCTACTCAATTTTGGTTTGCCCCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTg                           <  2:764993/151‑1 (MQ=255)
                                                                                                                   tcaccccacTTATCTACTCAATTTTGGTTTGCAGCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTg                           <  1:919477/145‑1 (MQ=255)
                                                                                                                     acacacTTATCTACTCAATTTTGGCTTCCCTCCGCAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACt                        <  2:148623/151‑1 (MQ=255)
                                                                                                                        caccTTATCTACCCGATTTGGGTTTGCACCCGAAAAACCCAGTCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTa                      <  2:736624/148‑1 (MQ=255)
                                                                                                                        acacTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTa                      >  2:862166/1‑151 (MQ=255)
                                                                                                                        acacTTATCTACCCAATTTTGGCTTCCACCCGAAAAAACCGGGTCGCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTa                      <  2:336960/151‑1 (MQ=255)
                                                                                                                        acacTTATCTACCCAATTTGGTTTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTa                      <  2:999410/151‑1 (MQ=255)
                                                                                                                         ccacttacctcaccaaTTTGGGCTGGCACCCGCAAAAACAAGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTa                      <  1:954763/137‑1 (MQ=255)
                                                                                                                         cacTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTaa                     >  2:601034/1‑151 (MQ=255)
                                                                                                                         cacTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTaa                     >  1:96593/1‑151 (MQ=255)
                                                                                                                         cacTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTaa                     >  1:451454/1‑151 (MQ=255)
                                                                                                                           cTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAata                   >  2:65460/1‑151 (MQ=255)
                                                                                                                           cTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAata                   >  2:1042342/1‑151 (MQ=255)
                                                                                                                           cTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAata                   >  2:1002689/1‑151 (MQ=255)
                                                                                                                            ttttctcCTCAATTTTGGCTTGCACCCAAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAatat                  <  2:246652/148‑1 (MQ=255)
                                                                                                                            ttATCTATTAAATTTTGGTTTTCACCCAAAAAAACAGGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAatat                  <  2:493429/151‑1 (MQ=255)
                                                                                                                            ttATCTACTAAATTTTGGCTTACACCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAatat                  <  1:357783/151‑1 (MQ=255)
                                                                                                                                 tACTCAATTTTGGCTTGCATCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCAACATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACAt             <  2:852542/151‑1 (MQ=255)
                                                                                                                                 tACTCAATTTGGGTTTCCATCCACAAAACCCAGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACAt             <  1:56256/151‑1 (MQ=255)
                                                                                                                                 aacgcaaaTTTGTCTTGCACCCGAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACAt             <  1:239906/145‑1 (MQ=255)
                                                                                                                                     caaaTTTGCCTTGCACCCGAAAAACCGAGGTAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAAc         <  1:853820/149‑1 (MQ=255)
                                                                                                                                       attttGGTTTGCATCCAAAAAAACCAGTCAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACat       <  2:278044/149‑1 (MQ=255)
                                                                                                                                       aaTTTGTCTTGCTCCCAAAAAACCCAGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACat       <  1:226426/151‑1 (MQ=255)
                                                                                                                                       aaTTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACat       >  2:439029/1‑151 (MQ=255)
                                                                                                                                              gCTTGCTCCCCAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTAc  <  2:15388/150‑1 (MQ=255)
                                                                                                                                               atttGCATCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTAc  <  1:686519/148‑1 (MQ=255)
                                                                                                                                                           |                                                                                                                                              
GTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTAC  >  NC_000913/3665553‑3665840

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