Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,196,204	G→A	74.5%	intergenic (‑105/+128)	thiC ← / ← rsd	phosphomethylpyrimidine synthase/regulator of sigma D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,196,204	0	G	A	74.5%	85.8 / 29.3	62	intergenic (‑105/+128)	thiC/rsd	phosphomethylpyrimidine synthase/regulator of sigma D
Reads supporting (aligned to +/- strand):  ref base G (8/8);  new base A (22/22);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

TTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTT  >  NC_000913/4196067‑4196346                                                                                                                                          |                                                                                                                                                 ttgttCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGGTGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGt                                                                                                                                             >  1:85070/1‑143 (MQ=255)    ttCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTAtt                                                                                                                                          >  2:639665/1‑143 (MQ=255)       gcggcgGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGc                                                                                                                                       <  1:642718/143‑1 (MQ=255)           cggGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGa                                                                                                                                   >  1:26756/1‑143 (MQ=255)           cggGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGa                                                                                                                                   >  1:368185/1‑143 (MQ=255)            gggTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGAt                                                                                                                                  <  1:910080/143‑1 (MQ=255)                 aGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTGGCCTGAGCAGGt                                                                                                                             >  1:146699/1‑143 (MQ=255)                   ttttGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAACCCACAGAAGAAGCAAGAGACGAGCGCCGCAGGGGCGATGCGATAGCAGATTAACCCTTGTTCCCTTCACAGGTATTAGCCGGATCAGGGTc                                                                                                                           >  2:786293/1‑143 (MQ=255)                        ttgCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCgg                                                                                                                      >  2:565215/1‑143 (MQ=255)                                   gCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATaa                                                                                                           >  1:700155/1‑143 (MQ=255)                                    cTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATaaa                                                                                                          >  1:1105750/1‑143 (MQ=255)                                     tCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCACTTCGCAGGTATTAGCCTGATCAGGTTCCGCCGATCCCGAATaacc                                                                                                         >  2:343537/1‑141 (MQ=255)                                      cATTCCAAAAGGTTACGGACGTGGCTTGTCAGCCGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTGGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACg                                                                                                        <  2:905073/143‑1 (MQ=255)                                         tCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtc                                                                                                     <  1:166364/143‑1 (MQ=255)                                           cAAAAAGTTACGGACGTGGCTTGTCAGACGACGGCTGAAGCAAGAGACGCTCGCCGCAGGGGCGCTGCGGTCGCGGCTTAACTCTTGTTCCCTCGCCGGGCGTTAGCCTGCTCAGGTGCCGCGGAGCCCGCATCCACgggcgc                                                                                                   >  1:40065/1‑139 (MQ=255)                                           cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGTGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtct‑c                                                                                                  <  1:98264/143‑1 (MQ=255)                                           cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtct‑c                                                                                                  >  1:499061/1‑143 (MQ=255)                                           cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtct‑c                                                                                                  >  1:531171/1‑143 (MQ=255)                                           cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtct‑c                                                                                                  >  1:125294/1‑143 (MQ=255)                                            aaaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑Ca                                                                                                 <  2:1050949/143‑1 (MQ=255)                                            aaaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGAGTAACTCTTGTTCCCTGCACAGGTGTTAGCCTGAACAGGTCCCGCGGATCCCGAGTAAACGCTCT‑Ca                                                                                                 >  2:69800/1‑143 (MQ=255)                                              caaGTTAAGGGCGTTGCTTGTCAGACTACGGACACAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTATCTCTTTTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGc                                                                                               <  1:398346/142‑1 (MQ=255)                                                      ggACGTGGCTTGGCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATGAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAg                                                                                       <  2:168648/143‑1 (MQ=255)                                                         cGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCa                                                                                    <  1:306597/143‑1 (MQ=255)                                                            ggCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTc                                                                                 >  1:961988/1‑143 (MQ=255)                                                              cTTCTCACACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCg                                                                               <  1:114053/143‑1 (MQ=255)                                                               ttGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGa                                                                              <  1:639665/143‑1 (MQ=255)                                                                 gTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGCCT‑CAGCCCGTCAAGGCCCCCCGAcc                                                                            >  1:951942/1‑142 (MQ=255)                                                                  tCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACaa                                                                           <  1:178529/143‑1 (MQ=255)                                                                    aGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGACTAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACACGa                                                                         >  2:917726/1‑143 (MQ=255)                                                                       cgacgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAAt                                                                      <  2:368185/143‑1 (MQ=255)                                                                        gacgGATGAAGCTAGAGACGTTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCGGATCAGGTGCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAAtt                                                                     <  1:129734/143‑1 (MQ=255)                                                                         acgGATGCAGCAAGAGACGCTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTa                                                                    <  2:744246/143‑1 (MQ=255)                                                                             aTGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGccc                                                                >  2:686944/1‑143 (MQ=255)                                                                              tGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCACTTCACAGGTCTTAGCCTCATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCGCTCCGACAGGAACTGAGcccc                                                               >  1:210804/1‑143 (MQ=255)                                                                                  gCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑aa                                                          >  1:900622/1‑143 (MQ=255)                                                                                       agaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATggg                                                     <  2:26756/143‑1 (MQ=255)                                                                                       agaCGACCGCCGCAGGGGCGATGCGATAGCAGATTAACTCGTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTGAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATggg                                                     <  2:740275/143‑1 (MQ=255)                                                                                        tcgcaTCGCCCCTGCGGCGATGCGATAGCAGATTAACTATTGTTCCCTTCACCGGTATTAGCCTGATCACGTGCCGCGGATCGCGAATAAACGGTCTCCGGCCAG‑TAAGGCACTCCGACAAGAAATTAGCCCCAC‑CAAGgggg                                                    >  2:513003/5‑143 (MQ=255)                                                                                        gacgacCCCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATgggg                                                    <  1:701745/143‑1 (MQ=255)                                                                                            aTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTCTTAGCCTGATCAGGTTCCGCGGCTCCCGAATCCACGGTCT‑CAGCCCGTTAAGGCACTCAGACAAGAAAATAACCCCGC‑AAACGGGACacg                                                >  1:126716/1‑141 (MQ=255)                                                                                            aTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACGGGTATTAGCCTGAGCAGGGTGCGCGGATCCCGAATAAACGGTCT‑CAGCCGGTTAAGGAGCTCCGGAGAGAAGTTGGCCCCGC‑AGAGGGGGCAtt                                                >  1:95876/1‑143 (MQ=255)                                                                                            aTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCACTACCCAGCTTATTGACTGATGAGGTTTATAGTTTCCCTATGAAACGGTCT‑CAGCCCGTTACGGCACTCCGACAAGAAATTAGCCCCGC‑AGATGGGGCAtg                                                >  1:257478/1‑142 (MQ=255)                                                                                                ccgcTGGCGCGATGCGCTCGCATACTAGCTCTTGTTCCCGCCACAGGTATTCGCCTGATCAGCTTCCGCGGATCCCGAATACACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAAt                                            <  2:576786/143‑1 (MQ=255)                                                                                                ccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCAAGCCACACCACTCCGACAAGAACGGAGCCCCGC‑AAATGGGGCATTGagt                                            >  2:658028/1‑141 (MQ=255)                                                                                                ccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTAATAGCCTGATCAGGTTCCGCGGAGCCCGAATAAACGGTGT‑CAGCCAGTTAAGGCACTGCGACAAGAAATTAGGCGCGC‑AAATGGGGCAtgggat                                            >  1:227502/1‑138 (MQ=255)                                                                                                 cgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATg                                           >  2:201887/1‑143 (MQ=255)                                                                                                     ggggCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATGCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAAATGGGGCCTTGAATGtca                                        >  1:155765/1‑141 (MQ=255)                                                                                                        gcgatgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTa                                    <  1:966950/143‑1 (MQ=255)                                                                                                          gatgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGACAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTAcg                                  >  2:297798/1‑143 (MQ=255)                                                                                                                    cAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAgc                        >  2:357243/1‑143 (MQ=255)                                                                                                                     aGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAgcg                       <  2:125294/143‑1 (MQ=255)                                                                                                                     aGATTAACTCTGGTTCCCTTCACAGGCATTAGGCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGGCAAGAAATTAGCCCGGC‑ACGTCGGTCACTGAATGTAAATTACGCGTTAACAgcg                       <  2:907176/143‑1 (MQ=255)                                                                                                                          aaCTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGaa                  <  2:700155/143‑1 (MQ=255)                                                                                                                          aaCTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAGACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGaa                  <  1:641454/143‑1 (MQ=255)                                                                                                                             tctTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCCTTGAAAGTAAAGTACGCGTTAACAGCGCAGAACTc               >  1:623780/1‑143 (MQ=255)                                                                                                                              ctTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCa              <  2:85070/143‑1 (MQ=255)                                                                                                                              ctTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCa              <  1:651859/143‑1 (MQ=255)                                                                                                                               tcgTTCCCTTCGCAGGTATTCGCCTTATCAGGTTTCGCGGATCCCTAATCAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGATATTAGCCCCTC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCaa             <  1:321084/141‑1 (MQ=255)                                                                                                                                 ggtCCCTTCGCAGGCATGCGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGc           <  1:293057/141‑1 (MQ=255)                                                                                                                                     ccTTCACAGGTAGTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACGCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGa       <  2:40242/143‑1 (MQ=255)                                                                                                                                       ttCGCTGGTGTTTGCCTGGTCATGTTCCGCGGATCCCGATTAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTATATTACGCGTTAACAGTGCAGAACTCAAGCAGGATg     <  2:402893/143‑1 (MQ=255)                                                                                                                                          gCAGGTATTAGCCTGATGAGGTTCCGCGGATCCCGAATAAGCGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGttt  <  1:300816/143‑1 (MQ=255)                                                                                                                                          gCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATACACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGttt  <  2:14232/143‑1 (MQ=255)                                                                                                                                          gCAGGTACTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGttt  <  1:144/143‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                 TTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCT‑CAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGC‑AAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTT  >  NC_000913/4196067‑4196346

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