Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,174,610	0	A	C	41.3%	102.5 / 14.8	64	W186G (TGG→GGG)	ycfT	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (27/10);  new base C (0/26);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.07e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTG  >  NC_000913/1174461‑1174752                                                                                                                                                        |                                                                                                                                               ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccac                                                                                                                                               <  2:617391/151‑1 (MQ=255)      cACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACgg                                                                                                                                          <  2:443060/151‑1 (MQ=255)      cACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACgg                                                                                                                                          <  2:101553/151‑1 (MQ=255)        ccgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCg                                                                                                                                        <  1:246180/151‑1 (MQ=255)          gaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTg                                                                                                                                      <  1:586972/151‑1 (MQ=255)           acagCCAGACCGGTAAGCAAAGAAGCTTTCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATAATAGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGAAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCg                                                                                                                                        >  2:130753/4‑148 (MQ=255)            ccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTggg                                                                                                                                    <  1:937461/151‑1 (MQ=255)                   aaCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaa                                                                                                                              >  2:693025/1‑150 (MQ=255)                    aCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAt                                                                                                                            >  2:1162/1‑151 (MQ=255)                             aaaGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGcc                                                                                                                    >  1:837300/1‑150 (MQ=255)                                  agCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTc                                                                                                              >  2:77340/1‑151 (MQ=255)                                  agCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTc                                                                                                              >  1:75208/1‑151 (MQ=255)                                   gCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTc                                                                                                              >  2:982280/1‑150 (MQ=255)                                         agcagATGGCGGCGCAGCGGCCCCTCTTTAACACAGGTCATTATTTTTGCCCAAAA‑CCAGGCGCAAA‑GGCTGTAATAAGGAAATTTGCGGATCACACTTTTCATCCCCCCCCCCGGCGGGGGAACAAAATTAACCGCCACACTcagca                                                                                                          <  2:1008357/148‑1 (MQ=255)                                         agcagAGGGCGGCGCACCGCCCCCTCTTAAACACAGGTTATTTTTTTTGCGCAAAA‑CCAGGCGCAAA‑GGCTGTAATAAGGAAAATGGCGGCTCACCCTTTTCATCCCCCCCCACGGCGGGGGAACAAAATAAACCGCCACACTCAGCAGTa                                                                                                       <  2:620190/151‑1 (MQ=255)                                               tGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaac                                                                                                  >  2:423252/1‑150 (MQ=255)                                                 gcggcgCAGCGGCCCCTTTTTAACACAGTTAATTATTGTTGCGCAAAA‑CCATGCGCCAA‑GGCTGTAAAAAGGCAAATGGCGGACCACACTTTTCTTTCCCCCCCGCGGCGGGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAg                                                                                               <  1:888937/151‑1 (MQ=255)                                                     cgcAACGGCCCCTTTTAAACACAGGTCATTATTG‑TGCCCCAAACCCAGGGGCAAG‑GGTTGAAAAAAGGAAAATTGCGGACCACCCTTTTCTTCCCCCCCCACGGCGGGGGAAAGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGa                                                                                           <  2:847050/151‑1 (MQ=255)                                                     cgcAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCg                                                                                            >  1:552904/1‑150 (MQ=255)                                                       cAACGCCCCCTCTTAAACACAGGTATTTATTGTTGCCCAAAC‑CCACCCGCCAA‑GGCTGTAAAAAGGAAAATTGCGGATCCCACTTTTCATCCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCg                                                                                            <  2:814946/148‑1 (MQ=255)                                                        accGGCCCCTTTTAAACACAGGTCATTATTGTTGCGCCAAC‑CCATGCGCAAA‑GGCTGTAATAAGGAAAATTGCGGATCACACTGTTCATTCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGaa                                                                                          <  1:866157/147‑1 (MQ=255)                                                        accGGCCCCTCTTTAACACAGGTCTTTTTTGTTGCGCAAAA‑CCATGCGCCAA‑GGCTGTAAAAAGGAAAATTGCGGATCACACTGTTCATCCCCCCCCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGa                                                                                           <  1:16095/146‑1 (MQ=255)                                                        aaCGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAAt                                                                                         >  1:557352/1‑150 (MQ=255)                                                             cACCTTTTTAACACGGGTCATTATTGTTGCGCAAAA‑CCAGGCGCCAA‑GGCTGTAATAAGACAAATGGCGGATCACTCTGTTCATTCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTg                                                                                     <  1:533216/149‑1 (MQ=255)                                                               ccTCTTTAACACAGGTCTTTATTTTTGCGCAAAA‑CCAAGCGCAAG‑GGCGGAAAAAGGGCAAATTGCGAATCACCGTTTTCATCCCCCTCCAGGGCGGGGGACCAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAg                                                                                 <  1:386589/151‑1 (MQ=255)                                                                ctctTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgc                                                                                >  2:938108/1‑151 (MQ=255)                                                                   tttAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                             >  1:845854/1‑151 (MQ=255)                                                                      aaCACAGGTCATTATTGTTGCGCAAAA‑CCTGGCGCAAG‑GGCTGTAAAAAGGAAAATTGCGGATCACCCTGTTCATTCCCCCCCGCGGGGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCa                                                                            <  2:1061184/149‑1 (MQ=255)                                                                        cacgggtattttttgtggCGCCAAC‑CCAGGCGCCAG‑GGTTGTAAAAAGGCAAATTGCGGATCACCCTGTTCATCCCCCACCCCGGCGTGGGAACAAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                        <  1:5409/139‑1 (MQ=255)                                                                         acggttgtttatgGTTGCGCAAAC‑CCTGGCACCAG‑GGGTGTAAAAAGCAAATTTGCGAATCCCACTGTCCTTCCCCCCCCACGGCGTGAGAGCGAATTTACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                        <  2:880255/142‑1 (MQ=255)                                                                          caGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGc                                                                      >  2:384261/1‑151 (MQ=255)                                                                                catttttttgccgcAAA‑CCCGGCCCCAAGGGTTGAAATAAGGCAAATGGCGGATCCCCCTTTTCTTCCCCCCCCAGGGGGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                                   <  2:596164/136‑1 (MQ=255)                                                                                atttttGTCGCGAAAAA‑CCCTGCACAAG‑GGTTGTAAAAAGGAAATTTGCGGATCACCCTGTTCATTCCCCCCCCGGGCGTGAGAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaa                                                                  <  1:472698/149‑1 (MQ=255)                                                                                 ttattGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAt                                                               >  2:275637/1‑151 (MQ=255)                                                                                 ttattGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAt                                                               >  2:283489/1‑151 (MQ=255)                                                                                     tgttgCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtttt                                                            >  2:784000/1‑150 (MQ=255)                                                                                      gttgCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTa                                                           >  2:783825/1‑150 (MQ=255)                                                                                          cgcAAAC‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACCCTTTTCATCCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATa                                                       <  1:378683/150‑1 (MQ=255)                                                                                             cAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCAc                                                   >  1:752779/1‑151 (MQ=255)                                                                                                   cATGCGCCAA‑GGCTGTAATAAGGAAAATTGCGGATCACACTGTTCATCCCCCCCCCCGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                              <  1:24228/151‑1 (MQ=255)                                                                                                    aTGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                              >  1:415394/1‑150 (MQ=255)                                                                                                       cgcCAG‑GGCTTAAAAAAGAAAAATTGCGTATCACACTGTTCATTCCCCCCCACGGCGGGGGAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGAc                                          <  2:819018/151‑1 (MQ=255)                                                                                                         caaa‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATCCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACa                                         <  1:723687/148‑1 (MQ=255)                                                                                                          cAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGATCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAt                                       >  1:952417/1‑151 (MQ=255)                                                                                                           aa‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAtt                                      >  2:426338/1‑151 (MQ=255)                                                                                                           aa‑GGCTGAAAAAAGGAAAATTGCGGACCACACTGTTCATCCCCCCCCACGGGGGGAGAAGGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACa                                         <  2:276464/148‑1 (MQ=255)                                                                                                              ggttGTAAAAAGGAAAATTGCGGATCCCACTGTTCATCCCCCACCACGGGGGGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTaa                                    <  2:84685/148‑1 (MQ=255)                                                                                                              ggctaaagaaaGGAAAATTGCGAATTACATTTTTCTTCCCCCCCCCCGGTTGGAGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                               <  2:126433/131‑1 (MQ=255)                                                                                                               gCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAAc                                   >  1:954112/1‑151 (MQ=255)                                                                                                                 tGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTAAAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGc                                 >  2:491789/1‑151 (MQ=255)                                                                                                                  gTAAAAAGAAAAATGGGGGATCACACTGTTCATTCCCCCCCAGGGGGGGGAAACAAAATAAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCa                                <  2:1019937/151‑1 (MQ=255)                                                                                                                       aaGGAAAATTGCGGACCCCACTTTTCATTCCCCCCCACGGCGTGGGAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAg                           <  2:811476/151‑1 (MQ=255)                                                                                                                        aGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                             >  2:9103/1‑148 (MQ=255)                                                                                                                            aaaTTGCGGATCCCACTTTTCATCCCCCCCCCCGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATa                        <  2:576667/149‑1 (MQ=255)                                                                                                                            aaaTTGCGGACCCCCTTGTTCATCCCCCACCCCGGCGTGGGAAAGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                      <  2:76720/151‑1 (MQ=255)                                                                                                                             aaTTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAc                       >  2:205531/1‑149 (MQ=255)                                                                                                                              aTTGCGGATCCCCCTTTACATCCCCCCCCACGGGGGGGGAAAGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcac                    <  2:146042/151‑1 (MQ=255)                                                                                                                              aTTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                      >  1:621406/1‑149 (MQ=255)                                                                                                                                     aTCACCCTGTTCATCCCCCCCCACGGCGTGGAAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTgg             <  2:352590/151‑1 (MQ=255)                                                                                                                                       cacaCTGTTCATCCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcaca                   <  1:255084/143‑1 (MQ=255)                                                                                                                                        acacTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGc          >  2:970366/1‑151 (MQ=255)                                                                                                                                         cccttttcttcccccccCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTg              <  1:591412/135‑1 (MQ=255)                                                                                                                                         cccTGTTCTTCCCCCCCCCCGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTg           <  1:259234/147‑1 (MQ=255)                                                                                                                                             gTTCATCCCCCACCCCGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTg        <  1:276427/148‑1 (MQ=255)                                                                                                                                              ttcaccccccACCAGGGCGTGGAAAAGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTg        <  2:152384/141‑1 (MQ=255)                                                                                                                                                 aTTCCCCCCCACGGCGTGGGAAAGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTg  <  2:383804/150‑1 (MQ=255)                                                                                                                                                    ccacaccagGGTGGGGAAACGAAGTAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTGTACATACCACGAAATAGGCAATTAACGCATACAGTTACCACAAGCTGGTGCTGGCGGt   <  2:484556/143‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                               CCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAA‑CCATGCGCCAA‑GGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTG  >  NC_000913/1174461‑1174752

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