Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,601,176	0	T	C	44.3%	82.5 / 42.1	71	E304G (GAA→GGA)	ftsX	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (30/9);  new base C (0/31);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 3.89e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGG  >  NC_000913/3601028‑3601326                                                                                                                                                        |                                                                                                                                                        ttattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTcc                                                                                                                                                       >  2:438291/1‑151 (MQ=255)  tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTcc                                                                                                                                                       <  2:1004012/150‑1 (MQ=255)  tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTcc                                                                                                                                                       <  1:1356863/150‑1 (MQ=255)   attCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCg                                                                                                                                                      <  2:153058/150‑1 (MQ=255)      cAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑aa                                                                                                                                                  <  2:320365/150‑1 (MQ=255)      cAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AAc                                                                                                                                                 <  2:534078/151‑1 (MQ=255)       aGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AAcc                                                                                                                                                >  1:1420145/1‑151 (MQ=255)        ggCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AAccg                                                                                                                                               <  2:1153529/151‑1 (MQ=255)            tAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AAccgccg                                                                                                                                            >  1:386383/1‑150 (MQ=255)                      tAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCgc                                                                                                                                 >  2:1408132/1‑151 (MQ=255)                       aaaTGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCgc                                                                                                                                 >  1:858829/1‑150 (MQ=255)                       aaaTGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAACGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCgc                                                                                                                                 >  1:1291369/1‑150 (MQ=255)                                    tgacAACCCCCGCCGCCCCCCAGCTAATCATCGAACACACCAGAAGCAATACCAGGCATTCATCAAATAATAACCAATTAATATCAACCTCCGTCCCAAAAC‑CC‑TGCC‑CCCCTCCCCA‑AACCGCCGATGACAATCGCAGCACCAGAAttt                                                                                                                    <  2:1072858/147‑1 (MQ=255)                                        aatccccccccccacccacCCAACAATCGATCACACACACAAAAATACCCGGCATACAGCAAATAACAACCTATTAATAAAAATCCCCGTCCCCAAAC‑CC‑TGTC‑CCACCCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGa                                                                                                                <  2:291379/139‑1 (MQ=255)                                           cccccctcccccccAGACAAACACCGAACACACCAGCAGAAAGAACAGGCATTCATAAAATAAAACCCCATTAATAACAACCTCCCTCCCAAAAA‑CC‑TGTG‑CCACTCCCCC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGa                                                                                                                <  2:45738/136‑1 (MQ=255)                                           ccACGCCCCCACCCAGCAAAACATAGAGCACACCAGCAGAAAAAGCAGACATACAAGAAAAAAAAAGCCATTGATAAAAACCTTC‑TTACAAAAA‑CC‑TGGGCCCCCCCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaaa                                                                                                            <  2:536417/151‑1 (MQ=255)                                            cccGCTCCCCCCCAGCAAATCATCGACAATCCCAGCAGAAATAGCAGGCATTCATCGAATAAAAAGCCATTGATATC‑AACTTTCTTCCAAAAA‑CC‑CGTG‑CCCCCTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTc                                                                                                                   <  1:556929/141‑1 (MQ=255)                                               gCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAtt                                                                                                         >  1:441755/1‑150 (MQ=255)                                                   ccccccAGCAAATCACCGAACACACCAGCAGAAATAGCAGCCATTCATAGAATAAAAACCAATTAATAACAACCTTCTTCCCAAAAC‑CC‑TGTC‑CCCCTCCCAC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaaaa                                                                                                           <  2:898892/141‑1 (MQ=255)                                                    cccccAGCCAATCACGGCGAATACCACAAGAAATGGCGGCCATTCATCGAATGATACGCCATTAATTAAAAACTTCTCCCAAAAAA‑CC‑TGAC‑CCGCTACCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATg                                                                                                     <  2:177560/147‑1 (MQ=255)                                                     aCCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGac                                                                                                   >  1:1182726/1‑150 (MQ=255)                                                      cccAACAAATCATCGAGCATCCCACCACAAATAGCAGGATTTCATAAAATAATAAGCAATTAATATAAACCTTC‑TTCCAAAAACCC‑TGTC‑CCCCTCCCAC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGacaa                                                                                                 <  1:1000385/151‑1 (MQ=255)                                                           caaaTCACCGACCATACCAGCAGCAAAAGCAGCCATTAATAAAATAATAAGCCTTTAATACAAACCTTCGTCCCAAAAC‑CC‑TGGC‑CCCCTCCCGA‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATg                                                                                            <  1:1031662/149‑1 (MQ=255)                                                              atcatcGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgcgc                                                                                         >  2:91192/1‑151 (MQ=255)                                                                    gAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGaa                                                                                    >  2:30019/1‑150 (MQ=255)                                                                           ccAGAAGCAATAGCAGACTTTCATCAAATGATAAGCCTTTAATATAAACCTTCGTCCCAAAAC‑CC‑TGTC‑CCCCTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCca                                                                             <  1:244259/150‑1 (MQ=255)                                                                               cagcaATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCcagcag                                                                         >  1:1173287/1‑150 (MQ=255)                                                                               cagcaATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGt                                                                        >  1:842074/1‑151 (MQ=255)                                                                                agcaATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTg                                                                       >  1:397214/1‑151 (MQ=255)                                                                                agaaaTAGCGGGCATTCATAAAATGAAAACCAATTAATATCAACCTTCTTCCCAAAAC‑CC‑TGTC‑CCACTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATccc                                                                              <  1:1143294/141‑1 (MQ=255)                                                                                agaaaTAGCAGGCATAAAGCAAATAAAAACCAATTGAAAAAAATCTTCGTCCCAAAAC‑CC‑TGTC‑CCACTCCCGA‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCcagca                                                                          <  2:744648/145‑1 (MQ=255)                                                                                  caATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgcg                                                                     >  1:631459/1‑151 (MQ=255)                                                                                  aaaTAGCGGGCTTCAATCGATTGATAAGCCATTTATAACAACCTCCTTTCAGAAAC‑CC‑GGTC‑CCCCTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgcg                                                                     <  1:1181589/150‑1 (MQ=255)                                                                                   aaaaGCAGGCATTCATCAAATAATAACCAATTAATATCAAACTTCCTCCCAAAAC‑CC‑TGTC‑CCCCTCCCGA‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCca                                                                             <  2:23964/139‑1 (MQ=255)                                                                                   aaaaGCAGGCATTCATCAAATAAAAAGCATTTAATAACAACCTTCTTCCCAAAAC‑CC‑TGTC‑CCCCTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgcgc                                                                    <  1:1026090/148‑1 (MQ=255)                                                                                         ggtctttCATATAATGAGAACCTATTGATAAAAAACTTC‑TCCCAAAAACCC‑TGCC‑CCTTGTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCAt                                                              <  1:1403270/146‑1 (MQ=255)                                                                                           gCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATAc                                                            >  1:995264/1‑151 (MQ=255)                                                                                               tCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAgg                                                         >  1:1149683/1‑150 (MQ=255)                                                                                               tCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAgg                                                         >  1:856662/1‑150 (MQ=255)                                                                                                 aacaaataatAACCCATTAATATAACACTTCTTCCCAAAAC‑CC‑TGTC‑CCCCTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGaa                                                       <  2:941416/146‑1 (MQ=255)                                                                                                 aTCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAAc                                                      >  2:203241/1‑151 (MQ=255)                                                                                                   cGAATAATAAGCCATTGATA‑AAAACCTCCTTCCAAAAACCC‑TGTG‑CCACTCCCCC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAAc                                                      <  1:636170/149‑1 (MQ=255)                                                                                                    gAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACggg                                                   >  2:956481/1‑151 (MQ=255)                                                                                                     aaTGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgc                                                  >  1:719096/1‑151 (MQ=255)                                                                                                     aaTGATAAGCCATTGATATAAACCTTCTTCCCAAAAC‑CC‑TGTC‑CCCCTCCCAC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACgg                                                    <  1:110812/149‑1 (MQ=255)                                                                                                     aaTGAAAACCCATTGATATCAACCTTCTTCCCAAAAC‑CC‑TGTC‑CCCCTCCCCC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAAct                                                     <  1:172273/148‑2 (MQ=255)                                                                                                      ataaaaccccATTAATATCAATCTCCGTCCCAAACC‑CC‑TGAC‑CCACCCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACgg                                                    <  2:902777/140‑1 (MQ=255)                                                                                                       tGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgcgc                                                >  1:899439/1‑151 (MQ=255)                                                                                                        aaTAAGCCATTAATACAAACCTTCTTCCCAAAAC‑CC‑TCTG‑CCCCCCCCAC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGaa                                                       <  2:982734/142‑1 (MQ=255)                                                                                                         aaaaGCAATTAATAAAAAACGTCGTACCAAAAC‑CC‑TGCC‑CCACCCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCa                                               <  2:386373/148‑1 (MQ=255)                                                                                                          tAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAgg                                             >  2:1155638/1‑151 (MQ=255)                                                                                                          tAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAgg                                             >  1:186065/1‑151 (MQ=255)                                                                                                            agctattaaaaaaaaaCTTCG‑TCCAAAAA‑CCTTGTC‑CCTCTCCCGA‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGAt                                           <  2:873797/138‑1 (MQ=255)                                                                                                            aGCCATTTATACAAAACTCCGTCCCAAAAC‑CC‑GGTC‑CCACTCCCAC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGAt                                           <  2:1224669/151‑1 (MQ=255)                                                                                                             gCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATg                                          >  1:77039/1‑151 (MQ=255)                                                                                                              ccATTAATAACAACCTTCGTCCAAAAAC‑CC‑TGTC‑CCCCTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGAt                                           <  1:412927/149‑1 (MQ=255)                                                                                                                    atatCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCAt                                   >  1:45187/1‑151 (MQ=255)                                                                                                                     atacAAACCTCCTTCCCAAAACCC‑TGTC‑CCCCTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCa                                    <  1:1193725/147‑1 (MQ=255)                                                                                                                      atCAAACTTCTTCCCAAAAC‑CC‑TGTC‑CCACTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATc                                  <  1:662704/150‑1 (MQ=255)                                                                                                                         aaccttcttCCCAAAAC‑CC‑TGCG‑CCCCTCCCAA‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATg                                          <  1:1006569/136‑1 (MQ=255)                                                                                                                          accttcttCCCAAAAC‑CC‑TGTC‑CCCCTCCCAC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTc                              <  1:340255/148‑1 (MQ=255)                                                                                                                          aaCCTCCTTCCCAAAACCC‑TGTC‑CCCCTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATc                                      <  1:299245/143‑1 (MQ=255)                                                                                                                                gTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCa                        >  2:1150051/1‑150 (MQ=255)                                                                                                                                  tCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATc                     >  2:607310/1‑151 (MQ=255)                                                                                                                                    caaaac‑CC‑TCTC‑CCCCTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCt                                 <  2:312332/135‑1 (MQ=255)                                                                                                                                    cGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATc                     >  1:502455/1‑149 (MQ=255)                                                                                                                                      aaac‑CC‑TGTC‑CCCCTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCa                           <  1:453347/141‑1 (MQ=255)                                                                                                                                      aaaa‑CC‑GGTC‑CCACTCCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCa                           <  2:184478/141‑1 (MQ=255)                                                                                                                                        aa‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTc               >  2:863234/1‑151 (MQ=255)                                                                                                                                        aa‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTc               >  1:172678/1‑151 (MQ=255)                                                                                                                                                   ccACTCCCCC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGtt       <  1:438303/151‑1 (MQ=255)                                                                                                                                                       ttCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATg   >  1:1144305/1‑151 (MQ=255)                                                                                                                                                        tCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATgg  >  1:905642/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                        TTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGC‑AACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGG  >  NC_000913/3601028‑3601326

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