Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,115,396	0	A	C	68.8%	75.7 / 56.4	61	L32R (CTG→CGG)	mdtG	putative drug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (9/10);  new base C (22/20);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

CGGTCGGCGAGTCCACCCCAAAACGGTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATA  >  NC_000913/1115252‑1115541                                                                                                                                                 |                                                                                                                                                  cGGTCGGCGAGTCCACCCCAAAACGGTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGa                                                                                                                                             <  2:55413/151‑1 (MQ=255)         gAGTCCACCCCAAAACGGTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAggcggc                                                                                                                                      >  2:61396/1‑150 (MQ=255)                ccccAAAACGGTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTa                                                                                                                               <  1:742849/150‑1 (MQ=255)                ccccAAAACGGTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTaa                                                                                                                              <  1:904785/151‑1 (MQ=255)                       aCGGTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAAc                                                                                                                        <  1:473224/150‑1 (MQ=255)                         ggTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGc                                                                                                                     >  1:739007/1‑151 (MQ=255)                            gAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTa                                                                                                                  <  1:405955/151‑1 (MQ=255)                               gCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGcc                                                                                                               >  1:212982/1‑151 (MQ=255)                               gCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGcc                                                                                                               >  2:25865/1‑151 (MQ=255)                                cGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATTTTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGcc                                                                                                               >  1:777786/1‑150 (MQ=255)                                    ggCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGc                                                                                                           <  1:360664/150‑1 (MQ=255)                                     gCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGcg                                                                                                          >  2:638899/1‑150 (MQ=255)                                      ccGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGcgac                                                                                                        >  2:358271/1‑151 (MQ=255)                                      ccGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGcgac                                                                                                        >  2:366715/1‑151 (MQ=255)                                      ccGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGcgac                                                                                                        <  2:732387/151‑1 (MQ=255)                                         aaaaTAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGAt                                                                                                     >  2:176723/1‑151 (MQ=255)                                         aaaaTAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGAt                                                                                                     >  1:768730/1‑151 (MQ=255)                                          aaaTAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATc                                                                                                    <  1:417347/151‑1 (MQ=255)                                            aTAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATc                                                                                                    >  2:245955/1‑149 (MQ=255)                                            aTAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAg                                                                                                  <  1:451567/151‑1 (MQ=255)                                              aaaaaTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAgg                                                                                                 <  2:108930/150‑1 (MQ=255)                                              aaaaaTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAgg                                                                                                 <  1:182444/150‑1 (MQ=255)                                               aaaaTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGgtt                                                                                               <  2:760622/151‑1 (MQ=255)                                                aaaTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGgttt                                                                                              <  1:183644/151‑1 (MQ=255)                                                 aaTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGgtttc                                                                                             >  2:394489/1‑151 (MQ=255)                                                  aTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGgtttcg                                                                                            <  1:757080/151‑1 (MQ=255)                                                  aTGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGgtttcg                                                                                            >  2:331052/1‑151 (MQ=255)                                                   tGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGgtttcgt                                                                                           >  1:837546/1‑151 (MQ=255)                                                    gTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGgtttcgt                                                                                           <  2:137499/150‑1 (MQ=255)                                                      aaTGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGgtttcgttt                                                                                         <  1:132712/150‑1 (MQ=255)                                                         gCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGgtttcgtttc                                                                                        >  2:411017/1‑148 (MQ=255)                                                         gCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAg                                                                                     <  1:795172/151‑1 (MQ=255)                                                           tGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGtt                                                                                   >  2:2003/1‑151 (MQ=255)                                                                  aaTACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAggg                                                                            >  1:745545/1‑151 (MQ=255)                                                                              caTATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTca                                                                <  1:912438/151‑1 (MQ=255)                                                                              caTATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTc                                                                 >  1:158117/1‑150 (MQ=255)                                                                               atatTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTcac                                                               >  2:137600/1‑151 (MQ=255)                                                                               atatTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTcac                                                               >  2:770999/1‑151 (MQ=255)                                                                                    cAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGt                                                          <  1:328374/151‑1 (MQ=255)                                                                                         cGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACAt                                                     <  1:97419/151‑1 (MQ=255)                                                                                         cGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACAt                                                     >  1:751122/1‑151 (MQ=255)                                                                                            aGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGc                                                  <  2:913509/151‑1 (MQ=255)                                                                                            aGTGACCGGGAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACGGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAg                                                   <  2:788092/150‑1 (MQ=255)                                                                                                aCCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATc                                              >  2:776794/1‑151 (MQ=255)                                                                                                 ccGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATcc                                             >  2:760394/1‑151 (MQ=255)                                                                                                  cGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCg                                            <  1:337510/151‑1 (MQ=255)                                                                                                    gTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCt                                          >  2:789478/1‑151 (MQ=255)                                                                                                     tAACGCCAAGCTGCTCAACGTAAAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTg                                         >  2:402936/1‑151 (MQ=255)                                                                                                         gCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTgg                                     >  2:632501/1‑151 (MQ=255)                                                                                                          ccAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGt                                    <  1:327200/151‑1 (MQ=255)                                                                                                           cAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTg                                   >  1:42861/1‑151 (MQ=255)                                                                                                                 gctCAACGTAGAGGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCa                              <  1:531426/150‑1 (MQ=255)                                                                                                                          agagGGGTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTg                     >  2:646660/1‑150 (MQ=255)                                                                                                                           gagGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGaa                   <  1:661516/151‑1 (MQ=255)                                                                                                                              gggTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACAt                <  1:118231/151‑1 (MQ=255)                                                                                                                               ggTAAGAAGGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATa               <  2:59557/151‑1 (MQ=255)                                                                                                                                  aagaagGGCATTACCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATaaga            <  2:756787/151‑1 (MQ=255)                                                                                                                                  aagaagGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATa               >  2:723966/1‑148 (MQ=255)                                                                                                                                   agaagGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATaaga            >  1:7229/1‑150 (MQ=255)                                                                                                                                       gggCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGaaa       <  2:669295/151‑1 (MQ=255)                                                                                                                                              aCCCGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATa  >  1:478508/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CGGTCGGCGAGTCCACCCCAAAACGGTGAGGCGATGGCCGAAAATAAAAATGTAATGCTGAAGACAATACCGGACCACATATTCAGGGCGGAGTGACCGGTAACGCCAAGCTGCTCAACGTAGAGGGGTAAGAAGGGCATTACCAGACTGAAGGCGGCACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATA  >  NC_000913/1115252‑1115541

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