Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,613,904	0	A	C	40.0%	102.5 / 14.0	60	V267G (GTA→GGA)	sad	succinate semialdehyde dehydrogenase, NAD(P)+‑dependent
Reads supporting (aligned to +/- strand):  ref base A (27/9);  new base C (0/24);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.09e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGC  >  NC_000913/1613755‑1614052                                                                                                                                                        |                                                                                                                                                     caCGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATa                                                                                                                                                      <  2:211034/150‑1 (MQ=255)  aCGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATAc                                                                                                                                                     <  1:607761/150‑1 (MQ=255)   cGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTTTCGCTGCACATACCt                                                                                                                                                   >  1:324509/1‑151 (MQ=255)   cGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCGCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAGCGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATa                                                                                                                                                      <  1:214380/148‑1 (MQ=255)      aaaTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTc                                                                                                                                                >  1:389847/1‑151 (MQ=255)       aaTCAAAACGAGCCATGGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTcc                                                                                                                                               <  1:400361/151‑1 (MQ=255)        aTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCg                                                                                                                                              <  2:981103/151‑1 (MQ=255)        aTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCg                                                                                                                                              <  2:302163/151‑1 (MQ=255)        aTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCg                                                                                                                                              >  2:432171/1‑151 (MQ=255)        aTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCg                                                                                                                                              <  1:274743/151‑1 (MQ=255)            aaaCGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTAt                                                                                                                                          >  2:648109/1‑151 (MQ=255)               cGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCt                                                                                                                                       >  1:640601/1‑151 (MQ=255)                gAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTg                                                                                                                                      >  1:420/1‑151 (MQ=255)                       tcccccGAGCGCGTTCTCTTCGTCACGAGGATCGCCCATCTCCAAGGCCGCCGCAGCTGCCACAAACCGCTCGGTAAA‑TG‑CCGCAGCAATTCCCTCTCCGATAATAAAGCGTTTTGCCGCTGCACATACCTGCCCGGTATTCTCATAAc                                                                                                                                 >  1:955908/5‑149 (MQ=255)                         gTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTcc                                                                                                                             <  2:1009664/151‑1 (MQ=255)                           ccGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCg                                                                                                                            <  2:823925/150‑1 (MQ=255)                                 gCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTAccg                                                                                                                     >  1:764704/1‑151 (MQ=255)                                   gTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTAccgcc                                                                                                                   >  1:457375/1‑151 (MQ=255)                                   gTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTAccgcc                                                                                                                   <  2:974748/151‑1 (MQ=255)                                      ctctTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCt                                                                                                                >  1:292599/1‑151 (MQ=255)                                       tctTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCtt                                                                                                               >  2:102199/1‑151 (MQ=255)                                         ttCGCCACGGGGATCGACCATTTTCAAGGCTGCCGCAGCTGCCGCAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGCATTCTGATAACGTCCGGCTACCGCCGCtct                                                                                                              >  1:409236/1‑148 (MQ=255)                                               aCGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGcc                                                                                                       >  2:123931/1‑151 (MQ=255)                                                cGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGcc                                                                                                       >  1:803661/1‑150 (MQ=255)                                                 ggggATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCa                                                                                                      <  2:289123/150‑1 (MQ=255)                                                 ggggATCGCCCATTTCCAAGGTTGCCGCACCGGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCCCATACCTGTCCGGTATTCTGATACCGTCCGGCTACCGCCGCTTTCACCg                                                                                                         <  2:226314/147‑1 (MQ=255)                                                 ggggACCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCa                                                                                                      <  2:923260/150‑1 (MQ=255)                                                      tCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTcc                                                                                                 >  1:340307/1‑150 (MQ=255)                                                                         ccccggCTCCCACAAACCGTTGGGTAAA‑TG‑CCAAAGAAATTCCCTCTTGGTAAAAAAAGCGTTTTGCCGCTCCACATCCCTGCCCGGTATTCTGATAACGTCGGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGcac                                                                             <  2:646193/146‑1 (MQ=255)                                                                         ccGCGGCTGCCACAAACCGTTGGGTAAA‑GG‑CCGAAGAAATCCCCTCTTCGAAAATAAAGCTTTTTGCCGCTGCACATCCCTGCCCGGTATTCTGAAAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAg                                                                                <  2:961272/148‑1 (MQ=255)                                                                          cGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGcac                                                                             >  2:59494/1‑150 (MQ=255)                                                                            tactcgCCAC‑AAACGTCTGGTAAATTC‑CCAAAGAAATCCCCTTTTCAAAAAAAAAGCTTTTTGCCTCCGCTCTTCCCTGCTGGGTATTCTTATACCTTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCTTTAAg                                                                                <  1:740726/140‑1 (MQ=255)                                                                             aGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAAt                                                                          >  1:763567/1‑150 (MQ=255)                                                                                  ccACAAAACGTTCGGTAAA‑TG‑CCGAAGAAATTCCCTCTTCGATAATAAAGCTTTTTGCCGCTGCACATCCCTTCCCGGTATTCTGATACCGTCCGGTTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATaaa                                                                       <  2:329644/148‑1 (MQ=255)                                                                                  ccAAAAAACGTTCGATAAT‑TC‑CCGAACCATTTCCCTCTTCAAAAAAAACGCTTTTTGCCGCTGCCCACCCCTCTCCGATATTCTAAAACCCACCGACACCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATa                                                                         <  1:201988/146‑1 (MQ=255)                                                                                      aaactgtgcGGTAAA‑TG‑CCGAGACAATTCCCTCTTGAATAAAAAAGTGTTTGGCCTCTCCACATCCCTGCCCGTTATTCTAATGAGTGCGGGCAACCGCCGGTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGAt                                                                  <  1:610540/141‑1 (MQ=255)                                                                                       aaaCCTTCGGTAAA‑TGCCCGAAGAAATCCCCTCTCCAAAAAAAAAGATTTTTGCCCCTCCACAACCCTTCCCGTTTTTATAAGAACGTCGGGCACCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGa                                                                   <  1:714394/148‑1 (MQ=255)                                                                                          cGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATACCGTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGaa                                                             <  2:700541/150‑1 (MQ=255)                                                                                          cGTTCGGAAAA‑TG‑CCGAAGAATTTCCCTCTTATAAAAAAAAGCTTTTTGCCGCTGCACATCCCTGCCCGTTTTTGTGAAAACGTCCGCATACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGa                                                              <  1:458063/149‑1 (MQ=255)                                                                                                  aaa‑TG‑CCAAAGCAATCCCCTCTCCAATAATAAAGCGTTTTGCCGCTCCACATCCCTGCCCGGTATTCTGAAAACCTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCa                                                      <  2:760511/149‑1 (MQ=255)                                                                                                  aaa‑TG‑CCAAAGAAATCCCCTCTTCAATAAAAAAGCTTTTTGCCGCTCCACATCCCTGCCCGGTATTCTGATAACTCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCa                                                      <  2:194238/149‑1 (MQ=255)                                                                                                      tG‑CCGAAGCAAGTCCCTCTTCGATAATAAAGCGTTTGGCCGCTGCACAGACCTGTACGGTATTCTGATACAGTCTGGCTACCGCCGCTTTCACAGCCAGATCCAGACCGTCATCGTTAAGCACAATAAACGTGGCCGAAGCGCCCAGTTcc                                                 >  1:399748/1‑151 (MQ=255)                                                                                                      atgccaaggaatttCCCTCTTATTTAAAAAGGCTTTTTGCCGCTCCCCTTCCCTGCCCGGAATTCTGAAAACGTCCGGCTCCCGCCTCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTc                                                  <  2:562349/139‑1 (MQ=255)                                                                                                               cAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGc                                          >  2:764621/1‑150 (MQ=255)                                                                                                                aaTCCCCTCTTCGAAAAAAAAGCTTTTTGCCGCTGCACATCCCGTTCCGGTATTCTAAAACCGCCCGTCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCAt                                        <  2:69542/151‑1 (MQ=255)                                                                                                                    cccTCTTGGAAAATAAGGCGTTTGGCCGCTCCCCATCCCTGCCCGTTATTCTGATAACGCCCGCCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCAt                                        <  2:86899/147‑1 (MQ=255)                                                                                                                     ccTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCAtttttt                                   >  1:1136423/1‑151 (MQ=255)                                                                                                                       tctTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCAtttttt                                   >  2:154924/1‑149 (MQ=255)                                                                                                                         tccGATAAAAAAGCTTTTTGCCGCTGCACATCCCTGTCCGGTTTTCTGAAAACGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTc                                  <  2:930088/146‑1 (MQ=255)                                                                                                                           caaaaaaaaaGCTTTTTCCCGCTCCACATCCCTGCCGGGTATTCTAAAAACGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTg                              <  2:1062071/143‑1 (MQ=255)                                                                                                                             ataataAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCg                           >  1:28550/1‑151 (MQ=255)                                                                                                                               aataaGGCTTTTTGCCGCTGCACATCCCTGTCGGTTATTCTAATACCGCCCGCCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCt                         <  2:435400/151‑1 (MQ=255)                                                                                                                               aataaAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCt                         >  2:62622/1‑151 (MQ=255)                                                                                                                               aaaaaaGCGTTTTGCCGCTGCACATCCCTGCCGGTTTTTCTAAAACCGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCg                           <  2:1039528/146‑1 (MQ=255)                                                                                                                                        ttttGCCGCTCCACATCCTTTTCGATTTTTCCAATAACGTCGGCCTCCCGCCGTTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAg                     <  1:585347/146‑1 (MQ=255)                                                                                                                                        ttttCCCGCTCCACATCCCTGCCGGGTATTCTAATAACGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGcc                   <  2:324699/148‑1 (MQ=255)                                                                                                                                         tttGCCGCTGCACATCCCTGCCGGGTTTTCTAATACCGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCtgt                <  2:406583/150‑1 (MQ=255)                                                                                                                                          ttGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGCCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCtgtg               <  2:328209/150‑1 (MQ=255)                                                                                                                                                cTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAAt        >  1:304951/1‑151 (MQ=255)                                                                                                                                                 tGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAAt        >  1:473657/1‑150 (MQ=255)                                                                                                                                                      aTACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGc  >  1:808642/1‑151 (MQ=255)                                                                                                                                                      aTACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGc  >  2:262392/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                     CACGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TG‑CCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGC  >  NC_000913/1613755‑1614052

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