Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,664,896	0	T	G	39.4%	102.5 / 17.1	70	G169G (GGT→GGG)	yfhR	putative S9 family prolyl oligopeptidase
Reads supporting (aligned to +/- strand):  ref base T (8/32);  new base G (26/0);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.15e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGCCAGCGTCAGCCCGATCCC  >  NC_000913/2664748‑2665042                                                                                                                                                     |                                                                                                                                                    aaGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑Ga                                                                                                                                                  <  2:1165927/151‑1 (MQ=255) aaGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑Ga                                                                                                                                                  <  1:869890/151‑1 (MQ=255)      aCGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTCTGGATGTTATTGGTCGGGGG‑GATCGt                                                                                                                                              >  1:869846/1‑150 (MQ=255)       cGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAAAATTCTGGGTGTTATTTGTCGGGGG‑GATCg                                                                                                                                               >  2:967970/1‑148 (MQ=255)           gTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGGGGATATTTTGGGTGTTTTTGGTCGGGGG‑GGTCGTg                                                                                                                                             >  2:374509/1‑146 (MQ=255)           gTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGc                                                                                                                                        <  2:561772/151‑1 (MQ=255)              ccAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGGGGATATTCTGGGTGGTATTTGGCGGGGG‑GGTCGTGAAgg                                                                                                                                         >  2:394008/1‑147 (MQ=255)                   ccGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGGTATTCTGGGAATTATTGGGCGGGGG‑GGCCGTGTAgg                                                                                                                                         >  1:1075349/1‑142 (MQ=255)                   ccGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGCATATTCTTGATGTTTTTGGTCGGGGG‑GATCGTGAACGCCGACGTg                                                                                                                                 >  1:260731/1‑150 (MQ=255)                   ccGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGc                                                                                                                                <  1:46577/151‑1 (MQ=255)                          gCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTCTTGGTGTTATTGGGTGGGGG‑GATCGGGAAGGCATaagcgggg                                                                                                                              >  2:316918/1‑139 (MQ=255)                          gCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTCTGGGTGTTTTTTGTCGGGGG‑GGGCGTGAAGGCATACGTGCGGGGagc                                                                                                                         >  2:1194825/1‑149 (MQ=255)                          gCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATc                                                                                                                         <  1:335775/151‑1 (MQ=255)                                cGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTCTGGGTGTTATTGGTCGGGGG‑GATCGTGAAGGCATACGTGCGGGGATCCTCGa                                                                                                                    >  1:78109/1‑150 (MQ=255)                                cGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGGTGTTTTTGGTCGGGGG‑GATCGTGGAGGCATACGTGCGGGGATTCTCGa                                                                                                                    >  1:1196642/1‑150 (MQ=255)                                      gCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAAAAATCTGGGTGTTATTTGGTGGGGG‑GGTCGGGGAGGCATACGGGCGGTGATCCTCGACTcc                                                                                                                >  1:1245203/1‑148 (MQ=255)                                       cAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGGGAATATTTTGGGTGTTATTGGTCGGGGG‑GGTCGTGAAGGCATACGTGCGGGGATCCTCGACTcccc                                                                                                              >  1:819678/1‑147 (MQ=255)                                        aaaGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGGTGTTATTGGTCGGGGG‑GATCGTGAAGGCATACGTGCGGGGATCCTCGACTCCACAtt                                                                                                           >  1:65126/1‑151 (MQ=255)                                        aaaGAGCCATCAATGTGGTGCGCCATCGCAGTGATGTCAAACCACAACGTCTGGTGCTGTTCGGGCAGAGCATTGGCGGGGCGAAAAATTTGGACGGTCTTGGTGGCGGA‑GGACGGGGAGGCCTAGCTGCGGGGAGCCACGACTcc                                                                                                                >  2:958799/1‑146 (MQ=255)                                         aaGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCCAGTATTCTGGATGTTATTGGTCTGGGT‑GGTCGTGAAGGGCTACGTCCGGGGATCctccactccaacttt                                                                                                          >  1:1276539/1‑146 (MQ=255)                                              ccATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCctc                                                                                                     <  1:684892/151‑1 (MQ=255)                                               cATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGG‑GGTCGTGAAGGCATACGTGCGGTGGTCCTCGACTCCCCATTTTCctct                                                                                                    >  1:46729/1‑151 (MQ=255)                                                    aTGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGTGCGGAAATTTTAGATGGTTTTTGTTGGGGG‑GATCGAGAAGGCATACCTGGGGTTATCCTTTACTTCACACTTGTCTCtttt                                                                                                 >  2:922695/1‑147 (MQ=255)                                                     tgtgGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTCTGGATGTTATTGGTCCGGGG‑GATCGTGAAGGGATACGTGCGGGGATCCTCGACTCCACATTTGCCTCtttt                                                                                                 >  2:115234/1‑146 (MQ=255)                                                        ggtgCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAAcc                                                                                           <  2:876965/151‑1 (MQ=255)                                                        ggtgCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAAc                                                                                            <  2:112965/150‑1 (MQ=255)                                                        gctacgcCATCGCAGTGATCTCATCCCACTACGTTTTGTGCTGTTCGGGCAGAGTATTTGGGGGGCGCATATTCTGCGTGTTATTGGTCCGGGT‑GATCCTGAAGGGATCCGTGCTTTGATCCTCGACTCCGCAGTTGCCCCCTATGCatg                                                                                            <  2:661903/146‑3 (MQ=255)                                                          tgCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATACTTTGGGTGTTATTGGTCGGGGG‑GGTTGTGAAGGCGTACGTGCGGTGGTCCCCGACTCCACATTTTGCTCTTAttctaac                                                                                           >  2:273892/1‑143 (MQ=255)                                                           gcgcCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGGGAATATTTTGGGGGTTTTTGGGCGGGGG‑GGGCGTGAAGGGATAAGGGGGGTGGTCCTCGGACCCCCATTTGTCTCTTTTGCAAccc                                                                                          >  2:765744/1‑148 (MQ=255)                                                           gcgcCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCAt                                                                                         <  1:579334/150‑1 (MQ=255)                                                           gcgcCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATc                                                                                        <  1:644208/151‑1 (MQ=255)                                                             gcCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGGGGATAATCTGGGTGTTATTTGTGGGGGG‑GGTCGTGGAGGCATACGTGCGGGGGTCCTCCACTCCACATTTGCCTCTTATTCAACCATCg                                                                                       >  2:1201364/1‑150 (MQ=255)                                                               cATCGCAGTGATGTGAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGcca                                                                                    <  2:870273/151‑1 (MQ=255)                                                                aTCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTCTGGATGGTTTTGGGCGGGGG‑GATTGTGAAAGCATACGTGCGGGGATCCTCGACTCCACCTTTTCCTCTTTTGTAACCCTCGcc                                                                                     >  1:1215846/1‑149 (MQ=255)                                                                aTCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaa                                                                                   <  2:1147410/151‑1 (MQ=255)                                                                  cGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAACATTTTGGTCGTTATTGGTCCGGGTGGCTCAT‑TAGGCGAACGGGCGGTGGTCCCCGCCTCCTCATTTTTCTCTTATTTAACCAACAccaa                                                                                   >  1:1031789/1‑149 (MQ=255)                                                                    cAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGGGCATATTCTGGATGTTTTTGGTGGGGGG‑GATCGTGAAAGCATACCGGCGGGGTTCCTCGACTCCACCTTTTCCTTTTATTCCACCATCGcccacc                                                                                 >  2:592945/1‑149 (MQ=255)                                                                       tGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATAATCTGGATGTTATTGGTAGGGGT‑GGTCGTGGAGGCATACGTGCGGTGTTCCTCGACTCCACATTTGCCTCTTATTCAAACAACCCCAACCaaa                                                                              >  2:232195/1‑149 (MQ=255)                                                                            tAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGG‑GATCGTGAAGGCATACGTGCGGGGATCCTCGACTCCACATTTTCCTCTTATGCAACCATCGCCAACCAAATGATccc                                                                       >  1:1286565/1‑151 (MQ=255)                                                                             aaaCCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATcccc                                                                      <  1:141672/151‑1 (MQ=255)                                                                                cccACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGGGGATATTTTGGATGGTATTTGGCGGGGG‑GGTCGGGAAAGCATACGTGCGGGGATCCTCCACTCCACATTTTCCTCTTATTCAACCATCGCCAACCAAAAGATaccc                                                                      >  1:1080515/1‑148 (MQ=255)                                                                                    caACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTg                                                               <  1:625124/151‑1 (MQ=255)                                                                                         tCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTCTGGGTGTTATTGGTCGGGGG‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAAACAAATGATCCCCGGCAGGGGCtac                                                          >  1:970458/1‑151 (MQ=255)                                                                                          cTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGG‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCtact                                                         >  1:1052497/1‑151 (MQ=255)                                                                                            ggTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCtactt                                                        <  1:521277/150‑1 (MQ=255)                                                                                            ggTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCtact                                                         <  2:289437/149‑1 (MQ=255)                                                                                               gCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCtacttact                                                     <  1:269534/150‑1 (MQ=255)                                                                                               gCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTacttactt                                                    <  2:426071/151‑1 (MQ=255)                                                                                                 tGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGGGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGGGGCTacttactt                                                    >  1:620539/1‑149 (MQ=255)                                                                                                  gTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTtg                                                   <  2:826213/149‑1 (MQ=255)                                                                                                     cGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGGCGGGGG‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTTCCTCTTATTCAACCATCGCCAACCAAATGATACCCGGCAGTGGCTacttactt                                                    >  1:316348/1‑145 (MQ=255)                                                                                                       ggCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGagag                                            <  1:766152/151‑1 (MQ=255)                                                                                                            agTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATGCGTGCGGTGATCCTCGACTCCTCATTTGCCTCTTATGCAGCCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTa                                         <  2:85323/149‑1 (MQ=255)                                                                                                                   gCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGa                                <  2:1134659/151‑1 (MQ=255)                                                                                                                    cGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGa                                >  2:596771/1‑150 (MQ=255)                                                                                                                      gggCGAATATTCTGGATGTTATTGGTCGGGGC‑GATCGTGAAGGCATACGTGCGGTGATCCTCGCCTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCCGCAGTGGCTACTTACTTGATGAGAGTTACAGCCGCGtaa                              >  2:148571/1‑150 (MQ=255)                                                                                                                         cGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATta                          <  1:374523/151‑1 (MQ=255)                                                                                                                          gAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATtat                         <  2:990940/151‑1 (MQ=255)                                                                                                                          gAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTGGCAGCGGCGAAAATtat                         <  2:709173/151‑1 (MQ=255)                                                                                                                            atatTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATtata                        <  2:363986/150‑1 (MQ=255)                                                                                                                             tatTCTGGATGTTATTGGTCGGGGC‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATc                      >  2:458411/1‑151 (MQ=255)                                                                                                                                   ggATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTTCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTCCTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGcc                   <  1:1115108/148‑1 (MQ=255)                                                                                                                                   ggATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGCCAGc                <  1:614866/151‑1 (MQ=255)                                                                                                                                       gTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGCCAGCGTCa            <  2:42990/151‑1 (MQ=255)                                                                                                                                         tATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGCCAGCGTCAGc          >  2:355428/1‑151 (MQ=255)                                                                                                                                            tGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGCCAGCGTCAGccc        <  2:781105/150‑1 (MQ=255)                                                                                                                                            tGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGCCAGCGTCAGCCCg       >  2:590297/1‑151 (MQ=255)                                                                                                                                              gTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCCCATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGCCAGCGTCAGCCCGAt     <  1:949820/151‑1 (MQ=255)                                                                                                                                               tCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGCCAGCGTCAGCCCGATc    >  1:745776/1‑151 (MQ=255)                                                                                                                                                  gggA‑GATCGTGAAGGCATACGCGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAACGATTCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGCCAGCGTCAGCCCGTTccc  >  2:1249141/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGT‑GATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTATATCGCCAGCGTCAGCCCGATCCC  >  NC_000913/2664748‑2665042

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: ‑