Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,903,692	0	C	A	43.1%	85.0 / 44.2	58	intergenic (+318/+1055)	ygcE/queE	putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (6/27);  new base A (25/0);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 4.74e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTC  >  NC_000913/2903555‑2903835                                                                                                                                           |                                                                                                                                                     atatTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTT‑AAAAACAC‑CC‑Atttgac                                                                                                                                        >  1:770766/1‑146 (MQ=255)    ttAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑AAAAACAC‑CC‑ATTAGcc                                                                                                                                        >  1:125742/1‑147 (MQ=255)        cTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAAACAC‑TC‑CTTAAACTC‑aaaaa                                                                                                                                >  1:935892/1‑149 (MQ=255)           cacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACaaa                                                                                                                             <  2:413766/150‑1 (MQ=255)           cacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAg                                                                                                                            <  2:445027/151‑1 (MQ=255)               cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACaa                                                                                                                        <  1:663626/151‑1 (MQ=255)                     cTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAACAC‑CC‑ACTAACCCC‑CAaaaaaaaaaaaaaatcc                                                                                                                  >  1:445318/1‑146 (MQ=255)                         aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGa                                                                                                              <  2:295780/150‑1 (MQ=255)                         aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGAt                                                                                                             <  1:971843/151‑1 (MQ=255)                         aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAAACAC‑CC‑ATTAGCCCC‑CAaaaaaaaaaaaa                                                                                                                       >  1:471563/1‑142 (MQ=255)                                 aGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTT‑TAAAACAC‑CC‑ATTAACCTC‑CAAAAAAAAAAAAAAAT‑CCCGATttatt                                                                                                        >  1:519443/1‑148 (MQ=255)                                  gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACtt                                                                                                   <  2:617586/151‑1 (MQ=255)                                        aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTATCCaa                                                                                             <  2:838805/151‑1 (MQ=255)                                          gAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑AAAAACAC‑TC‑ATTTACCTCAAAAAAAAAAAAAAATTT‑CCAATTTTTTT‑ATTTAAcaa                                                                                              >  2:584853/1‑147 (MQ=255)                                          gAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATg                                                                                           <  2:448556/151‑1 (MQ=255)                                           aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGa                                                                                          <  2:855080/151‑1 (MQ=255)                                             tCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑AAAAACAC‑CC‑ATTAAACTCAAAAAAAAAAAAAAAATT‑CCAAATTTTTA‑ATTTAACaaaaaa                                                                                          >  1:773084/1‑147 (MQ=255)                                              cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAAt                                                                                        <  1:610180/150‑1 (MQ=255)                                              cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAAt                                                                                        <  1:682626/150‑1 (MQ=255)                                                 aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAACAC‑TC‑CTTAACCTC‑AAAACAAAAACAAAATT‑CCAGATTTATT‑ATTTTAACAAAGAaaaa                                                                                      >  1:118394/1‑149 (MQ=255)                                                   tttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGCTTTTT‑TAAAAAAA‑TC‑AATTACCTC‑AAAAAAAAAAAAAAAAT‑CCAAAATTTTT‑TTTTTAACAAAAAATAACaa                                                                                   >  2:765830/1‑150 (MQ=255)                                                      aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCt                                                                                <  2:9266/150‑1 (MQ=255)                                                             ttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTaaaaaa                                                                         <  1:203151/150‑1 (MQ=255)                                                              tCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCTGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAat                                                                        <  2:263744/150‑1 (MQ=255)                                                                   tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑TAAAACAC‑TC‑ATTAGCCTC‑CAAAAAAAAACAAAAAT‑CCCGATTTTTT‑TATTAAACAAAGAATAAAAAACTTaaaaaaaa                                                                       >  1:13377/1‑144 (MQ=255)                                                                    aaaGTATGAACTGGAAATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGCTTTGTAAAAAACAT‑CA‑ATTGGCCCC‑AAAA‑AAAAAAAAAAAT‑CCACAATTATT‑AATTAACAAAAGGAAAAAAAACTTAAAAAATATACttt                                                                 >  1:464785/1‑149 (MQ=255)                                                                        tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAAAAC‑AC‑ATTAACCTC‑AAAAAAAAAACAAAAAT‑CCCCAATTTTT‑TATTAAACAAAGAAAAAAAAACTTAAAAAAAATACTGTtttt                                                             >  2:933507/1‑151 (MQ=255)                                                                        tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑AAAAACAC‑CC‑ATTAACCTC‑AAAAAAAAAAAAAAAAT‑CCCGATTTATT‑TATTTAACAAAAAAAAAAAAACTTAAAAAAaataattttttt                                                             >  1:10250/1‑143 (MQ=255)                                                                        tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTT‑AAAAACAA‑AC‑CTTAGGCTC‑CAAAAAAAAAAAAAATA‑TCAGATTTTTTAAATTAAAAAATAAAAAACAAAGTAaaaaaaatacctttttt                                                              >  1:291303/1‑140 (MQ=255)                                                                          tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAAACAC‑CC‑ATTTACCTC‑CAAAAAAAAAAAAAAAT‑CCAGAATTTTT‑TCTTTACCAAAGAATAACAAACTTtaaaaaaaaaattttttttt                                                           >  1:386071/1‑137 (MQ=255)                                                                           gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAACAC‑AC‑ATTTTCCTC‑AAAACAAAAAAAAAAAT‑CCCGATTTTTT‑TCTTTACCAATGAAAAAAAAAATTTAAAAAAATAATGTtttttt                                                           >  2:940515/1‑148 (MQ=255)                                                                            aaCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaa                                                         <  2:756567/151‑1 (MQ=255)                                                                                ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaa                                                     <  1:286360/151‑1 (MQ=255)                                                                                ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGGTATGT‑GAAACCAC‑AC‑ATTTGCAAC‑TAAACAAAAAAAAACAT‑TCAGATTTTTT‑TCTTTAACAATTTAAAAAAAACCTAAAAAAGAAAATGTtttttt                                                           >  1:595849/1‑143 (MQ=255)                                                                                 gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑AAAAACAC‑CC‑AATAACCTC‑AAAAAAAAAACAAAAAT‑CCAGATTTTTT‑TATTAAACAAAGAATAAAAAAATTAAAAAAAATAATGTTTTTTTTAAGaa                                                     >  1:777521/1‑150 (MQ=255)                                                                                 gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTAT‑TAAAACACACC‑ATTA‑CCCC‑CAAACAAAAAAAAACAT‑CCCAATTTTTT‑TTTTTAACAACTAAAAAAAAACTTTAAAAAAATACATTTATTTTAAAGaa                                                     >  2:37856/1‑150 (MQ=255)                                                                                    atCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGTTTAAAAAATATACTGTTATTCTAAAGaaaaa                                                  <  1:132799/150‑1 (MQ=255)                                                                                     tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAACAC‑TC‑ATTAGCCTC‑AAAAAAAAAACAAAAAT‑CCAGATTTATT‑TTTTAACCAATGAATAAAAAAATTAAAAAAAATAATGTTTTTTTTAAGaaaaaa                                                 >  1:331213/1‑150 (MQ=255)                                                                                     tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAACAA‑TC‑CTTATCCTC‑AAAAAAAAAAAAAAATTCCCAAATTTTTA‑ATTTAACAAATAAAAAAAAACCTTAAAAAA‑ATAATTTTTTTTTTAaaaaaaaa                                                 >  2:119754/1‑150 (MQ=255)                                                                                     tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAAAAA‑TC‑ATTTTCCCC‑AAAACAAAAAAAAAACT‑TCAAATTTATT‑AATTAAACAAATAAAAAAAAAATTAAAAAAAATAATTTTTTTTTTAaaaaaaaa                                                 >  2:392246/1‑150 (MQ=255)                                                                                        ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAAtt                                               <  2:252695/149‑1 (MQ=255)                                                                                         tATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTT‑TAAAAAAC‑ACAATTAACC‑C‑AAAACAAAAAAAAAATT‑TCAAATTTTTT‑TATTTAACAAAGAATAAAAAACTTAAAAAAAAAAATTTTATTTTTAAAAAAAAAATaa                                             >  2:757560/1‑150 (MQ=255)                                                                                             cATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGT‑TAAAAAAC‑AC‑AATTACCTC‑AAAACAAAAAAAAAAAT‑CCACATTTTTT‑AATTAAACAATAAAAAAAAAAACTAAAAAAAATAATCTTATTTTAAAAAAAAAAATaaaaagc                                        >  2:526990/1‑146 (MQ=255)                                                                                              aTAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCa                                       <  2:255791/151‑1 (MQ=255)                                                                                                     cTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAAACAC‑TC‑CTTAACCCC‑AAAAAAAAAAAAAAATT‑TCAGATTTTTT‑TATTTAACAAAAAATAAAAAAATTTAAAAAAAAAATGTTATTTTAAAAAAAAAAATAAAAAGCATTaaat                                 >  1:109739/1‑148 (MQ=255)                                                                                                         tcatcaACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAat                             <  2:147771/150‑1 (MQ=255)                                                                                                           atcaACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATAtt                          <  1:685991/151‑1 (MQ=255)                                                                                                             caACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTa                         <  1:529174/150‑1 (MQ=255)                                                                                                                 aGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                    <  2:248826/151‑1 (MQ=255)                                                                                                                      ataGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaataat               <  2:519301/151‑1 (MQ=255)                                                                                                                         gagGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGa            <  2:434704/151‑1 (MQ=255)                                                                                                                         gagGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑AATTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaataat               >  2:579340/1‑148 (MQ=255)                                                                                                                           ggCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAtt          <  1:963388/151‑1 (MQ=255)                                                                                                                             cGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTa         >  1:476198/1‑150 (MQ=255)                                                                                                                              ggTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaaa       <  1:122803/151‑1 (MQ=255)                                                                                                                              ggTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaa        >  1:267485/1‑150 (MQ=255)                                                                                                                                ttATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAGAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAt     >  2:539341/1‑151 (MQ=255)                                                                                                                                   tGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTc  >  1:940984/1‑151 (MQ=255)                                                                                                                                           |                                                                                                                                                     ATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TC‑ATTAGCCTC‑AAAACAAAGACAAACTT‑CCAGATTTATT‑ACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTC  >  NC_000913/2903555‑2903835

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