Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	40.0%	79.0 / 11.0	56	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (8/25);  new base G (22/0);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 4.79e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCG  >  NC_000913/3508514‑3508804                                                                                                                                                     |                                                                                                                                                    aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑gaaaa                                                                                                                                                 >  1:198353/1‑148 (MQ=255)   gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑G‑C‑C‑GGG‑GGAAATTca                                                                                                                                             >  2:359074/1‑150 (MQ=255)    ttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑G‑CGG‑GGG‑AAAAAtta                                                                                                                                              >  2:221283/1‑148 (MQ=255)     tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGATGC‑G‑G‑G‑GGG‑aaaaaataa                                                                                                                                             >  2:439023/1‑144 (MQ=255)        ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑G‑C‑G‑CGG‑GAAAAAacacacac                                                                                                                                        >  2:665817/1‑150 (MQ=255)         tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACg                                                                                                                                       <  1:257148/150‑1 (MQ=255)          ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAAATCacacaccg                                                                                                                                      >  1:304917/1‑148 (MQ=255)            aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAg                                                                                                                                   <  2:418794/151‑1 (MQ=255)             caGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑G‑C‑G‑CGG‑GAAAAATTACAAAAGGt                                                                                                                                     >  2:344030/2‑148 (MQ=255)                ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGGGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGcccc                                                                                                                               <  2:606764/151‑1 (MQ=255)                   gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑G‑C‑G‑CGG‑GAAAAATTACACACGGGAGccccc                                                                                                                              >  1:302264/1‑149 (MQ=255)                          cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGa                                                                                                                      <  2:58567/150‑1 (MQ=255)                           aaTATGCAGATGCTCATGGGCGAGGGTTTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGTTGGC‑G‑G‑G‑CGG‑AAAAATTCTCACACGGTAGCCCCCGCATCAGaaa                                                                                                                    <  2:324472/151‑1 (MQ=255)                           aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACGCACGGTAGCCCCCACATCAGaaa                                                                                                                    <  2:563509/151‑1 (MQ=255)                             tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTCGCGG‑G‑G‑GGG‑AAAAATTACCAC‑CGGTAAGCCCCCCATTAGgaa                                                                                                                    >  2:207140/1‑149 (MQ=255)                                gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑G‑G‑CGG‑GAAAATACACAAACAGGTGACCCCCCACCAGAAAGCggg                                                                                                               >  1:212891/1‑149 (MQ=255)                                  aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑GAAAAATCAAACACCGGAGACCCCCCCTCAGAAAACgggg                                                                                                              >  1:648287/1‑150 (MQ=255)                                      gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAGACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcgc                                                                                                         <  2:417384/151‑1 (MQ=255)                                       cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑GAAAATTTACACACCGTAGCCCCCCCATCAGAAAACgtggtg                                                                                                            >  1:301155/1‑147 (MQ=255)                                            gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑GAAAAATCACACACGGTAGCCCCCCCCACAGAAAACGGGGGGCGCACAcc                                                                                                    >  1:683999/1‑150 (MQ=255)                                             ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑GAAAAATTACAAACGGTAGCCCCCCCATCAGAAAGAGTGGGGCGCACAcc                                                                                                    >  1:823366/1‑149 (MQ=255)                                                         aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑GCCGG‑GAAAATT‑AAACACCGGAGACCCCACATCTAAAAGGGGGGGGGGGAAACCCGGCTGCGGAGCg                                                                                       >  2:189974/1‑150 (MQ=255)                                                            cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑GAAAAATTACACACGGGAGCCCCCCCATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGgtc                                                                                    >  1:788049/1‑150 (MQ=255)                                                             gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑G‑G‑CGG‑GAAAAATTACACACGGGAGCCCCCCAATCAGAAAGAGTGGGGGGCGCAccccgccgccttccg                                                                                       >  1:185143/1‑133 (MQ=255)                                                              tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGT‑G‑C‑G‑CGG‑TGAAATTCAAACACGGGAGTCCCCCCATCAGGAAGCGGGGCGGGCTCACCGAGGTGccgaatgat                                                                                     >  1:240740/1‑139 (MQ=255)                                                                 gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCt                                                                              <  2:853884/151‑1 (MQ=255)                                                                    cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑G‑C‑G‑CGGAAAAAATTAACA‑ACGGGAACCCCCCCAACCAAAAACGGGGGGCGGCCACCCAGATGCggtggggttatccc                                                                               >  1:649261/1‑134 (MQ=255)                                                                       aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGa                                                                        <  2:372725/151‑1 (MQ=255)                                                                           cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑G‑C‑G‑GGG‑GAAAAATCACAAACCGGAAACCCCACATCAGAAAAGGTGGGGGGGACACCCAGCTGCGTCGCGAGCATCCCTTTCttgt                                                                       >  2:165541/1‑145 (MQ=255)                                                                           cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑GAAAAATCAAAAACCGTTGACCCCACAACAAAAAGCGGGGGGCGGACACCGAGGTGCGTAAGGATAATCCCCTGCTGATaa                                                                     >  1:448303/1‑150 (MQ=255)                                                                            aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑GAAAAATTACACACGGTAGCCCCCACCTCAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTTCTGATaaa                                                                    >  1:475826/1‑150 (MQ=255)                                                                                ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGC‑GAC‑G‑CGG‑AAAAATT‑ACACAACGTACCCCCCCCCCCAGAAAACGAGGGGCGCCCAACCACCTGCGTCGCGAACCTCCCCTGCTTATAAAACa                                                                 >  2:249839/1‑149 (MQ=255)                                                                                ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                                                               <  1:821086/151‑1 (MQ=255)                                                                                  tgtCTCCTGCAAACCGGAGGGAGGAAAAAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCTCATCAGAAAGCGTTGTGCGCACACCTAGCTGCGTAACGATCATCCCTTGCTGATAAACCAGCa                                                              <  2:614755/150‑1 (MQ=255)                                                                                      tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑GGG‑GAAAATTCACAAACGGGAGCCCCCCCATCAGCAAGGGTGGTGGGGACACCCAGCTGCGGAGCGATCACCCCTTGCTGATTAACAAACACTgg                                                          >  2:267946/1‑150 (MQ=255)                                                                                      tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGG‑G‑C‑G‑CGG‑GAAAAATAACACAAGGGGGCCCCCCCCTACGAAAGCGGGGGGGGGGCCCCGAGCTGGGTGGCGAGCATCCCCTTCTgctcaaacaccccc                                                            >  2:104086/1‑135 (MQ=255)                                                                                      tcCTGCAAACCGGAGGGAAGAAAGAAAAGAGCGAACATTTTGTGGCGGTGGC‑G‑C‑G‑GGG‑AGAAATTCACGCACGGAAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCGTGCTGATAAACCAGCACTg                                                           <  2:171187/149‑1 (MQ=255)                                                                                       cgtGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                                                        <  2:142331/149‑1 (MQ=255)                                                                                       ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                          <  1:282010/149‑1 (MQ=255)                                                                                       ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                                                        <  1:814993/151‑1 (MQ=255)                                                                                       ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGAGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                                                           <  2:22748/148‑1 (MQ=255)                                                                                          gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca                                                     <  1:800835/151‑1 (MQ=255)                                                                                             aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑G‑G‑CGG‑GAAAAATCACAAACGGGACGCCCCCCACTAGAAAGCGGGGGGGGCGCACCCAGCCGGGTTGCGAACATCCCTTTTTGAAAAAACAGCACCGGc                                                         >  2:605708/1‑144 (MQ=255)                                                                                                cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑G‑C‑G‑CGG‑GGAAAATCACACACGGGAACCCCCCCCTCAGAAAGCGGGGGGCGGCCACCCAGCCGCGGGGCGGAAATCCCCTTCTGCTAAAACAACACCTGcaacac                                                    >  2:192723/1‑146 (MQ=255)                                                                                                cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑GAAAAATAACAAACGGTCGCCCCCCCCACCGAAAGCGTGGGGCGCACAACCAGCTGCGGAGCGGACATCCCTTTCTGATAAACCAGCACTGGcaccaccttcc                                               >  2:806705/1‑146 (MQ=255)                                                                                                 ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑GAAAAATTACACACGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGAACATCCCTTGCTGATAAAACAGCACCGGCACCACCTc                                                 >  1:127942/1‑148 (MQ=255)                                                                                                 ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                              <  2:799315/151‑1 (MQ=255)                                                                                                     ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGT‑G‑C‑G‑CGG‑GGAAAAACAAAAACGGCAGCCCCCACACCAAAAAGCGGGGTGCGCACACCGAGCCGCGGGGCGATCAACCCTTTCTGATAAAACAGCACCGGGACCACCTCCCggggg                                          >  2:683275/1‑149 (MQ=255)                                                                                                                   agCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAACCCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                            <  2:411431/151‑1 (MQ=255)                                                                                                                    gCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                           <  2:437008/151‑1 (MQ=255)                                                                                                                      gACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAAt                         <  2:802924/151‑1 (MQ=255)                                                                                                                          aTTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                     <  1:207199/151‑1 (MQ=255)                                                                                                                             ttGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                   <  2:294240/150‑1 (MQ=255)                                                                                                                               gtggcggtggc‑g‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGATCAATGGTCGGAt                 <  1:147111/150‑1 (MQ=255)                                                                                                                                tggcggtggc‑g‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg               >  1:408943/1‑151 (MQ=255)                                                                                                                                tggcggtggc‑g‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg               >  1:849840/1‑151 (MQ=255)                                                                                                                                         c‑g‑c‑g‑cgG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg       <  2:758601/150‑1 (MQ=255)                                                                                                                                                 cgG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg  >  1:253539/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑G‑C‑G‑CGG‑AAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCG  >  NC_000913/3508514‑3508804

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