Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,460,777	0	A	C	41.7%	96.2 / 13.7	61	T321P (ACC→CCC)	gspF	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base A (24/11);  new base C (1/24);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 2.18e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGTATTTTGCAATCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTAACGCTCTCCATCTTCGAGCCAGC  >  NC_000913/3460628‑3460925                                                                                                                                                      |                                                                                                                                                     tAAGTATTTTGCAATCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAAc                                                                                                                                                     >  1:45285/1‑151 (MQ=255)  aaGTATTTTGCAATCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAAc                                                                                                                                                     <  2:137924/150‑1 (MQ=255)     tATTTTGCAATCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGc                                                                                                                                                  >  1:818468/1‑150 (MQ=255)        tttGCAATCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAAt                                                                                                                                               <  1:1168240/150‑1 (MQ=255)          tGCAATCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAAttt                                                                                                                                             <  1:93655/150‑1 (MQ=255)           gCAATCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTc                                                                                                                                           >  1:437074/1‑151 (MQ=255)             aaTCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTccc                                                                                                                                         <  1:1086465/151‑1 (MQ=255)              aTCCAGCGTCGTCCCTCGGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTGGATATTCGCCAGCGTCTGGCAAATGAGGCAGAGAAAGATCGCCAGCGTAACAGCATTCATCTTTCGCTGGAACAACCCGCAATTTTCccg                                                                                                                                        >  2:263394/1‑151 (MQ=255)              aTCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTccc                                                                                                                                         >  2:74504/1‑150 (MQ=255)               tCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCccgc                                                                                                                                       >  2:460117/1‑151 (MQ=255)                 cAGCGGCGTCCCTCTGCTGGAGGGGATGAATTTGTCCACCGAAAGCCTCAACACCCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCccgccg                                                                                                                                     <  2:1193375/151‑1 (MQ=255)                   gcggcgTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCga                                                                                                                                    <  2:987370/150‑1 (MQ=255)                    cggcgTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAACCCGCAATTTTCCCGCCga                                                                                                                                    <  2:74065/149‑1 (MQ=255)                        gTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGc                                                                                                                              >  1:1119668/1‑151 (MQ=255)                        gTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGc                                                                                                                              >  2:375831/1‑151 (MQ=255)                                 ctgGATGGGATGAATTTGCCCCCCAAAACCCTCAACACCCCCAAAATCCGCCAGCGTCTGGAAAAGGCGGAAGAGAACTTTCGCCAGGGTAAAAGCATTCATTTTTCGCTGAAACAACCCGAATTTTTCCCGCCGATGATGCTCTAc                                                                                                                         <  2:417504/147‑1 (MQ=255)                                       gggatgAATTTGCCCACCAAAAGCCCAAACAACCTCGAATTCCGCCAGCGTCGGGCAAATGCGGCAGAGACCGTTCGCCAGGGTACCAGCATTCATCTTTCGCTGGAACAACCCGCAATCTTCCCGCCGATGATGCTCTACAtggtgg                                                                                                                  <  1:33426/148‑1 (MQ=255)                                             attttGCCCCCCGAAACCCTAAACACCCCCAAAATCCCCCAGCGTCTGGAAAAGGCGGCAGAGAACGTTCGCCGGGGAAACAGCTTTCTTTTTTCGCGGAAACAACCCGCATTTTTCCCGCCGATGATGCTCTACATGGTGGCctc                                                                                                              <  1:652674/144‑1 (MQ=255)                                              ttttGTCCCCCAAAACCCTCAACACCCTCGAAATGCGCCGGCGTCTGAAAATGGGGGCAGAAACCTTTCGCCGGGTAAACACCATCCCTCTTTCTGTGGAACAACCCTAATTTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTgg                                                                                                           <  2:344400/147‑1 (MQ=255)                                               ttttccccccaaaggCCTAAACACCCTGAAAATTCCCCGGCGTCGGGAAAATGCGGCAGAGACCGTGCGCCGGGTAAACAGCATTCTTCTTTCGTGGGAAAAACCCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTg                                                                                                            <  1:582597/132‑1 (MQ=255)                                               tttGTCCCCCGAAAGCCTCAACAACCTCAAAATTCGCCAGCGTCTGGAAAAGGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAACCCGCATTTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGa                                                                                                        <  1:486860/150‑1 (MQ=255)                                                  gTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGaaaa                                                                                                     >  2:1036474/1‑150 (MQ=255)                                                      aCCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGc                                                                                                 >  2:1181482/1‑150 (MQ=255)                                                      aCCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCg                                                                                                >  2:806000/1‑151 (MQ=255)                                                      aCCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCg                                                                                                >  1:534321/1‑151 (MQ=255)                                                       ccaaaaGCCACAACACCCTCAAATTCCGCCAGCGTCGGGAAAAGGGGGCAGAGAACTTTCGCCGGGATAACAGCATTCAGTTTCCGGTGGAAAAACCCGAAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCg                                                                                                <  2:464754/147‑1 (MQ=255)                                                         gAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGc                                                                                             >  2:719082/1‑151 (MQ=255)                                                          aaaGCCTCAACACCCTGAAAATTCGCCGGCGTCTGGCAAATGCGGCAGAGACCGTTCGCAAGGGTAACAGCATTCTTCTTTCGCTGGAACAAACCGCATTTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCgg                                                                                               <  1:601995/148‑1 (MQ=255)                                                          aaaGCCACAACACCCTCGAAATTCGCCGGCTTCTGAAAAATGCGGCAGAGACCGTCCGCCAGGTTAACACCATTCATTTTCCGCTGAAACAACCCAAATTTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCggg                                                                                              <  2:831116/149‑1 (MQ=255)                                                           aaGCCTCAACACCCTCAAATTTCGCCAGCTTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATTTTTCGCTGAAACAACCCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAg                                                                                           <  1:1174391/151‑1 (MQ=255)                                                              ccTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTc                                                                                        >  2:1190870/1‑151 (MQ=255)                                                               cTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCg                                                                                       >  1:892784/1‑151 (MQ=255)                                                                     acccTCGAAATTCGCCAGCGTCTGGCAAATGCGGGAGAGAACTTGCGCCAGGTTAACAGCATTCATTTTTCGCTGGAAAAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACAt                                                                                 <  1:455953/149‑1 (MQ=255)                                                                            aaaTTCGCCGGCGTCTGGAAAAGGCGGCAGAGAACGTTCGCCGGGGTAACAGCTTTCATTTTTCGCTGAAACAACCCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATg                                                                            <  1:1136943/149‑1 (MQ=255)                                                                            aaaTTCGCCGGCGTCTGGAAAAGGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATTTTTCGCTGAAACAACCCGAATTTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATgg                                                                           <  1:447997/150‑1 (MQ=255)                                                                            aaaTTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGt                                                                          >  1:1012699/1‑151 (MQ=255)                                                                                   ccggCGTCGGGCAAATGGGGCAGAGAACGTTCGCCAGGGAAACAGCATTCATCTTTCGCTGAAACAACCCGCATTTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCaga                                                                      <  1:970944/145‑1 (MQ=255)                                                                                   ccAGCGTCTGGAAAATGCGGAAGAGACCGTCCGCCGGGGAAACAGCATTCTTCTTTCGCGGGAAAAACCCGCATTTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGcc                                                                   <  2:23790/151‑1 (MQ=255)                                                                                    cAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCg                                                                  >  2:885233/1‑151 (MQ=255)                                                                                        gcccggCAAATGAGGGAGGAAACGTGCCTCGGGGTAACTGCATTCTTTTTTCGTTGGAACAACCCAAATTTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGcc                                                                   <  1:854636/142‑1 (MQ=255)                                                                                           tGGAAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCTTTCATCTTTCGCGGGAAAAACCCGAAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGTGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAg                                                               <  1:884091/147‑1 (MQ=255)                                                                                               aaaTGCGGCAGAGAACGTTCGCCGGGGTAACAGCATTCATCTTTCGCGGGAACAACCCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAg                                                       <  1:378003/151‑1 (MQ=255)                                                                                               aaaGGCGGCAGAGAACGTCCGCCGGGAAAACAGCATTCTTTTTTCGCTGAAACAACCCAAAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAg                                                       <  2:1066314/151‑1 (MQ=255)                                                                                                 aTGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAgg                                                      >  2:175760/1‑150 (MQ=255)                                                                                                     ggCAGAGACCTTTCGACAGGTTAACATCTTCCTTTTTTCGCTGAAACAAACCTCAATTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGgag                                                    <  1:96248/148‑1 (MQ=255)                                                                                                     ggCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAACCCGCAATTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGaca                                                 <  1:892802/151‑1 (MQ=255)                                                                                                                  cGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATc                                    >  1:137212/1‑151 (MQ=255)                                                                                                                       gggTAACAGCATTCTTTTTCCGCTGAAACAAACCGCATTTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATc                               <  1:1233288/151‑1 (MQ=255)                                                                                                                       gggTAACACCTTTCATCTTTCGCGGGAACAACCCGCATTTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGa                                 <  1:915010/149‑1 (MQ=255)                                                                                                                       gggAAACAGAATTCATTTTTCGCTGAAACAACCCGAATTTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGAt                                <  1:751537/150‑1 (MQ=255)                                                                                                                         gTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGc                             >  1:854189/1‑151 (MQ=255)                                                                                                                          tAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGcc                            >  2:277563/1‑151 (MQ=255)                                                                                                                           aaCAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCt                           >  1:457931/1‑151 (MQ=255)                                                                                                                              agaaTTCATTTTTCGCTGAAAAAAACCGCATTTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTa                         <  2:75361/147‑1 (MQ=255)                                                                                                                              aGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTa                         >  2:218842/1‑150 (MQ=255)                                                                                                                                  ttcttcttcCGCTGGAACAACCCGCAATTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTa                         <  1:1062227/146‑1 (MQ=255)                                                                                                                                        ttcggcgggAAAAACCCGAAATTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTAACGCTCTCCa               <  2:294874/143‑1 (MQ=255)                                                                                                                                        ttccGCTGAAACAACCCGCAATTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTAACGCtctc                 <  2:820052/145‑1 (MQ=255)                                                                                                                                         ttCGCTGAAAAAATCTGCAATTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTAACGCTCTCCAt              <  1:467366/150‑1 (MQ=255)                                                                                                                                           cgcggaaacaacCCGAAATTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTAACGCTCTcc                <  2:842321/140‑1 (MQ=255)                                                                                                                                            gcggaaacaacCCGTATTTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTAACGCTCTCCATCt            <  1:968990/144‑1 (MQ=255)                                                                                                                                              tGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTAACGCTCTCCATCTTCg         >  2:1048569/1‑150 (MQ=255)                                                                                                                                                  aCAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGACAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTAACGCTCTCCATCTTCGAGcc     >  1:901762/1‑150 (MQ=255)                                                                                                                                                    aaaCCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTAACGCTCTCCATCTTCGAGCCAGc  >  1:1014545/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TAAGTATTTTGCAATCCAGCGGCGTCCCTCTGCTGGATGGGATGAATTTGTCCACCGAAAGCCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGCAATTTTCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGATAACCAGGAGACACTCCAACAAAATCGGATCGCCTTAACGCTCTCCATCTTCGAGCCAGC  >  NC_000913/3460628‑3460925

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