Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,060,264	A→C	8.7%	intergenic (+34/‑183)	bipA → / → yihL	ribosome‑dependent GTPase, ribosome assembly factor/putative transcriptional regulator YihL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,060,264	0	A	C	8.7%	132.0 / 7.4	63	intergenic (+34/‑183)	bipA/yihL	ribosome‑dependent GTPase, ribosome assembly factor/putative transcriptional regulator YihL
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base C (3/3);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-02

GATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGC  >  NC_000913/4060123‑4060406                                                                                                                                              |                                                                                                                                                gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt                                                                                                                                                >  2:565620/1‑143 (MQ=255)      tgccGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                                                           <  1:531481/140‑1 (MQ=255)      tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAACCGATCGTCGGCGCGCCAACCGCGCACCGAAAGAGGATTAATTTCGTTTTTTTGTCATAATCGACCTGCGCGCGATgtcgg                                                                                                                                           >  1:309794/1‑139 (MQ=255)       gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGCCGGAAAACGCACGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTGTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                                                                          >  1:765802/1‑143 (MQ=255)       gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAgg                                                                                                                                          <  1:242091/143‑1 (MQ=255)        aCGAACTGGTAGAAGTGGCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCGGCGCCAACCGCGCGCCGAAAGACGATTAATTTCGTTGTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                                                                         >  1:74736/1‑143 (MQ=255)          gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGAAGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAGTAATTTCGTGTGTTAGTCATCAAAAACCGCCCAGCGATGGCAggctc                                                                                                                                       >  2:816600/1‑140 (MQ=255)           aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                                                                      >  1:62852/1‑143 (MQ=255)              tGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGCCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAAGAATGGCGTGTTTTAGTCATACAAAACCTGCCAGCGATGGGAGGttgtctt                                                                                                                                   >  2:756324/1‑138 (MQ=255)               gggaGAAGGGACTCCGACCTCTATCCGTATTCGTAAACGTCAACTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttat                                                                                                                                  <  2:332134/140‑3 (MQ=255)               ggTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCCCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttat                                                                                                                                  <  1:400596/143‑3 (MQ=255)                gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATACAAAACCTGCCAGCGATGGCAGGgttttattt                                                                                                                                 >  2:829910/1‑143 (MQ=255)                gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttt                                                                                                                                 <  1:62697/143‑1 (MQ=255)                 tAGAAGTGACTCCGACCTCTAGCCGTATTCGTAAACGGCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                                                                <  2:854754/143‑1 (MQ=255)                  aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttttc                                                                                                                               >  1:230236/1‑141 (MQ=255)                      gTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTTTTCTTTTTTTACCTGCCAGCGATGGCAGGTGTTTTCTTGCCTGt                                                                                                                           >  2:989595/1‑143 (MQ=255)                         aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGt                                                                                                                        <  2:447758/143‑1 (MQ=255)                           tCCGACCGCTATCCGTATTCGTAAGCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGa                                                                                                                      <  1:366115/143‑1 (MQ=255)                             cttccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgc                                                                                                                    <  1:181815/140‑1 (MQ=255)                              ggccTCTATCCGCATTCGTAAGCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcg                                                                                                                   <  1:593015/141‑1 (MQ=255)                                ccgctcgccGGATGCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGca                                                                                                                 <  1:855052/136‑1 (MQ=255)                                ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGCTGGCAGGTTTTTTTTTGCCTCTCGTGAGCGca                                                                                                                 >  1:524776/1‑143 (MQ=255)                                 ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcac                                                                                                                <  2:671167/143‑1 (MQ=255)                                  tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCGAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcaca                                                                                                               >  1:795495/1‑143 (MQ=255)                                    tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                                                             >  1:968928/1‑143 (MQ=255)                                    tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGAGTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGATTATTTTTGCATGTCGTGAGCCCacaca                                                                                                             >  2:611623/1‑142 (MQ=255)                                          tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTATTTAGTCATAAAAAACCTGCCAGCGATGGCAGGGTTCTTTTTGCCTGTCGTGCGCGGCCGCCAATGcc                                                                                                       >  2:1027364/1‑142 (MQ=255)                                           aTTCGTAATCCTCACCTGACGGAAAACGATCGTCCCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAGAAAAGACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                      <  2:411956/143‑1 (MQ=255)                                           aTTCGTAAACGTCACCTGACGGAAAACGCTAGTCGCCCCACCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                      <  1:724897/143‑1 (MQ=255)                                           aTTCGTAAACGTCACCCGACCCAAACCGCTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                      <  2:386231/143‑1 (MQ=255)                                                   aCGTCACCTCACGGGAAACCGTCGTCGCCGCGACAACCACGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAGCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTcc                                                                                             <  2:832523/143‑1 (MQ=255)                                                    cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAATCCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTccc                                                                                            <  2:154017/143‑1 (MQ=255)                                                     gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCt                                                                                           <  2:62852/143‑1 (MQ=255)                                                     gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCt                                                                                           >  2:841087/1‑143 (MQ=255)                                                      tCACCTGACGGAACACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCtt                                                                                          <  2:1076421/143‑1 (MQ=255)                                                      tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACTAATGCGAGTAT‑CTACCtt                                                                                          >  1:828305/1‑143 (MQ=255)                                                      tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGGTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCACTGCGAGTGC‑CTcccct                                                                                          >  2:137687/1‑141 (MQ=255)                                                           tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCCCACCAATGCGAGTAT‑CTTCCTTCTAAg                                                                                     >  2:18193/1‑143 (MQ=255)                                                               ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGTTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCtg                                                                                 >  1:246504/1‑143 (MQ=255)                                                                  aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAACCGACCAACTGCGCTATTTAGTCATAAAAAACCTGCCAGCGATGGCAGGGTGTATGTTGCCAGACGTGAGCGCACGCAAATGCGAGGAT‑CTCCCTTATCAGTcgcggg                                                                              >  2:81459/1‑138 (MQ=255)                                                                       aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTTAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGt                                                                        <  1:839394/143‑1 (MQ=255)                                                                       aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGGTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAGGTCTGTGGGTACg                                                                         >  2:1065666/1‑143 (MQ=255)                                                                             gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTTCGTCCAtg                                                                   >  1:410067/1‑142 (MQ=255)                                                                              ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTTTCTTGAGAGCACACCAATGCGAGTAC‑GTCCGTTATCAGGAGGTGGTTACGCCACTac                                                                  >  2:577360/1‑142 (MQ=255)                                                                                       ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                         >  2:776239/1‑143 (MQ=255)                                                                                        cgcgCAGCGAAAGACGATTAACTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt                                                        <  1:352602/143‑1 (MQ=255)                                                                                        cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt                                                        <  2:274006/143‑1 (MQ=255)                                                                                        cgcgCACCGAAAGACGATTAATTTCGTTTCTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACCCCAATGCGTGTAT‑CTGCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt                                                        <  1:295560/143‑1 (MQ=255)                                                                                                                    ttttAGTCATAAAGAACCTGCCAGCGATGGCAGGTTTTGTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCCATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGt                            <  1:905275/143‑1 (MQ=255)                                                                                                                      ttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGCGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGttc                          >  1:200041/1‑142 (MQ=255)                                                                                                                        cgTCATAGAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAt                        <  2:810783/142‑1 (MQ=255)                                                                                                                            aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGct                    >  1:156445/1‑143 (MQ=255)                                                                                                                            aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAA‑CTCCCTTATAGGTCTGTGGTTACGTCAATAGAGAGCTGTGTCACAATTATCTGCAAAGTCATACGCCGttgctctct                    >  2:402922/1‑136 (MQ=255)                                                                                                                              aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGcttt                  <  1:841087/143‑1 (MQ=255)                                                                                                                               aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttc                 <  1:998424/143‑1 (MQ=255)                                                                                                                                 aaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCATGAGCGCACACCAATGCGAGTAT‑CTACCTTATAAGTATGTGGTTACGTCCATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGATAATTGctttctt               >  2:248787/1‑143 (MQ=255)                                                                                                                                   aCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttt             >  1:870480/1‑143 (MQ=255)                                                                                                                                    ccTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCGATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCACAGTCTTACGACGTTAACTGCTTTCttttt            >  2:444484/1‑143 (MQ=255)                                                                                                                                     ctcccAGCGACTGCAGGTCTCTTTTTGCCTGTCGTGAGCTCACACCACTGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt           <  2:637551/140‑1 (MQ=255)                                                                                                                                     cggCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt           <  1:460815/141‑1 (MQ=255)                                                                                                                                        ccAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc        <  2:230236/143‑1 (MQ=255)                                                                                                                                        ccAGCGATGGCAGGTTTTTCTCCGCCTGTCGTGCGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc        <  1:539465/143‑1 (MQ=255)                                                                                                                                        ccAGCGATGGCAGGTTTGTTTTTGCCTGTCGTGAGCGCACACCAATG‑TACAATAATATGTCATTATTTTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc        <  1:492249/143‑1 (MQ=255)                                                                                                                                              aTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagc  <  1:613113/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                GATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAT‑CTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGC  >  NC_000913/4060123‑4060406

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