Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,662	0	T	A	41.7%	98.3 / 21.8	61	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (4/31);  new base A (25/0);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 4.57e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508519‑3508809                                                                                                                                                 |                                                                                                                                                     gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTcacac                                                                                                                                                <  2:936222/150‑1 (MQ=255)             ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACA‑CGGGAGccccc                                                                                                                                   >  2:232328/1‑150 (MQ=255)               agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACA‑CGGGGGccccccc                                                                                                                                 >  1:946915/1‑148 (MQ=255)                gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGGAAAAATTACACA‑CGGTAGGCCCcac                                                                                                                                 >  2:875110/1‑149 (MQ=255)                 agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATc                                                                                                                              <  2:309310/151‑1 (MQ=255)                      aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGaaa                                                                                                                         <  2:857894/151‑1 (MQ=255)                       atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGGAAAAATCACACA‑CCGGAGCCCCCCCATCAGaaa                                                                                                                         >  2:589573/1‑150 (MQ=255)                       atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATAACAAA‑CGGGAGCCCCCCCATCCGaaa                                                                                                                         >  1:581040/1‑150 (MQ=255)                           gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCgtg                                                                                                                    <  2:290008/151‑1 (MQ=255)                           gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGGGGGGGAAAAATTACAACA‑CGG‑AGCCCCCCCCTCAGCAAGCggg                                                                                                                    >  1:188247/1‑149 (MQ=255)                            cAGATGCTCATGGGCGAGTGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGAGCATTCTGTTGCGGTGTC‑GCGCGTAAAAAGTCACACC‑CGGTAGCCCCCACATCAGAAAGCgt                                                                                                                     <  1:1021480/149‑1 (MQ=255)                             aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGGAAAAATCACAAA‑AGGGAGCCCCCCCCTCAGAAAACggggg                                                                                                                  >  1:124775/1‑150 (MQ=255)                               aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTgc                                                                                                                <  1:800465/151‑1 (MQ=255)                                    cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGC‑GCGCGGGAAAAATTACACC‑CGGGAGCCCCCCCATCAGACAGCgggggggg                                                                                                               >  1:326732/1‑143 (MQ=255)                                           gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                    <  2:933875/151‑1 (MQ=255)                                           gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGC‑GCGCGGGAAAAATCACACA‑ACGTAGCCCCCCCATCAGAAAGCgcgggggg                                                                                                               >  2:574396/1‑136 (MQ=255)                                            aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg                                                                                                   <  1:923061/151‑1 (MQ=255)                                                   aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGGGAAAAACACACA‑AAGTAGCCCCCCCCACACAAAGCGGGGGGCGcgcaaccgactgcccggc                                                                                             >  1:50382/1‑136 (MQ=255)                                                   aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                                                            <  2:926313/151‑1 (MQ=255)                                                      ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                         <  2:877625/151‑1 (MQ=255)                                                       cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACA‑CGGGAGCCCCCACCTCAGAAAGCGGGGGGCGCACACCGAACTGCg                                                                                                 >  1:939182/1‑142 (MQ=255)                                                          cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGGAAAAATCACACA‑CGGTAGCCCCCACATCAGAAAGAGTGGGGCGCACACCGAGCTGCGTAGCGATCATcc                                                                                     >  2:214585/1‑151 (MQ=255)                                                            gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGGAAAAATCAAACA‑CCGGAGCCCCCACATCAGAAAACGGGGTGGGGACACCGGGCTGGGCAGGGGTCATccc                                                                                    >  2:833089/1‑150 (MQ=255)                                                              tCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACA‑CGGTAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCGAGCTGCGTAGCGATCATcccct                                                                                  >  2:877432/1‑148 (MQ=255)                                                                aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGG‑GCGCGGAAAAAATCACAAA‑CCGGAGCCCCCCCATCAGAAAAGGGGGTGCCCACACCCAGCTGCGTAGCGAACATccc                                                                                    >  2:817014/1‑146 (MQ=255)                                                                 aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                                                              <  2:264564/151‑1 (MQ=255)                                                                   aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAt                                                                            <  1:790583/151‑1 (MQ=255)                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                       <  2:1025992/150‑1 (MQ=255)                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                      <  1:831189/151‑1 (MQ=255)                                                                            gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                                                                   <  1:432064/151‑1 (MQ=255)                                                                             tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GGGCGGGAAAAATTACAAA‑AGGGAGCCCCCCCATCAGAAAAGGGGGGGGGGACACCCAGCGGCGTAGGGAACCTCCCTTTCTGATaaaa                                                                        >  1:736524/1‑144 (MQ=255)                                                                              gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGAAc                                                                  <  2:507064/150‑1 (MQ=255)                                                                                   cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATTAAAAA‑CGGGAGCCCCCACATCAGGAAACGGGGGGCGCCACCCCCGCTGCGGAGCGATCCTACCCTTCTGATAAAACACCACCgggcc                                                            >  1:827243/1‑148 (MQ=255)                                                                                    tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGGGGGAAAAATAACCCACCGGAGGCCCCCAC‑TCATAAACCGGGGGGCGGACACCGCCCTGCGGCGGGATCATCCCTTGCTGTTGAACcaccaccagcgcc                                                           >  1:654908/1‑144 (MQ=255)                                                                                       aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGGAAAAATAAAAAA‑CAGGAGCCCCCACATCCGAAAGCGGGGTGGGCACACCGAGCTGCGTAGCGATCATCCCTTTGTGATAAACCAACAACGGcac                                                            >  1:926751/1‑147 (MQ=255)                                                                                       aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcaccac                                                         <  1:858229/150‑1 (MQ=255)                                                                                         aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                      <  2:124858/151‑1 (MQ=255)                                                                                            ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                                   <  1:648538/151‑1 (MQ=255)                                                                                            ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                                   <  2:663444/151‑1 (MQ=255)                                                                                             gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                                  <  2:275793/151‑1 (MQ=255)                                                                                              agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATTACACA‑CGGTAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCTGCGGAGCGGTCATCCCCTGCTGATAAACCAGCACCGGCACCACCTCGCggg                                                 >  2:454164/1‑151 (MQ=255)                                                                                                 gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACA‑AGGTAGCCCCCCCATCAGAAAGCGGGGTGCGCACACCCAGCTGCGTAGCGATCAACCCTTGCTGATAAACCAGCACTTGcaccccctttcgggg                                                >  2:729575/1‑138 (MQ=255)                                                                                                   ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg                                            <  2:119381/151‑1 (MQ=255)                                                                                                     aaagaaaagaGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtg                                          <  2:1015320/151‑1 (MQ=255)                                                                                                      aagaaaagaGCGACCATTTTGTGGCGGTGGG‑GCGCGGGAAAAATCACAAA‑CGGGAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCCCGCTGCGGAGCGGTCATCCCTTGCTGATAAAACAACACTGGCCCCCCTTCGCGGGGAtggtgg                                         >  1:1013970/1‑151 (MQ=255)                                                                                                          aaagaGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACA‑CAGGAGCCCCCCCCTCCGAAAGGGGGGGGGGCACACCCAGCCGCGTAGCGATCCTACCCTGCTGATAAACCAACACTGGCACCACCTCGCGGGGATGGGGGGtg                                      >  1:50920/1‑150 (MQ=255)                                                                                                              agCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATTACACA‑CGGGAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACCTCGCGGGGATGGGGGGTGTGAcc                                 >  2:351158/1‑151 (MQ=255)                                                                                                                    cATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGt                           <  1:501849/151‑1 (MQ=255)                                                                                                                       tttGTGGCGGTGGT‑GCGCGGAAAAAATCCCCCA‑CGGGAGGCCCCACCTCAGGAAGCGGGGGGGGCACACCGAGCTGCGTAGCGATCATTCCTTTCTGATAAACCAACCCTGGCCACCCTTCTCGGGCAGGGGGGCGGTGagccaaggtggg                        >  1:721621/1‑140 (MQ=255)                                                                                                                        ttGTGGCGGTGGG‑GCGGGGGAAAAATTACAAA‑AGGTAGGCCCCACATTAGAAAGCGGGGGGGGCAAACCGAGCTGCGTAGGGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGATGGTGGGTGTGGCCAATGGGCgg                        >  1:619316/1‑150 (MQ=255)                                                                                                                         tgtgGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGa                       <  1:761450/150‑1 (MQ=255)                                                                                                                          gtggcggtggc‑gCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGAt                      <  2:55061/150‑1 (MQ=255)                                                                                                                            ggcggtggc‑gCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                   <  1:129001/151‑1 (MQ=255)                                                                                                                              cggtggc‑gCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg                 <  1:417560/151‑1 (MQ=255)                                                                                                                               ggtggc‑gCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGc                >  2:1031779/1‑151 (MQ=255)                                                                                                                                gtggc‑gCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc               <  2:322708/151‑1 (MQ=255)                                                                                                                                       cgcgGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc         <  2:318064/151‑1 (MQ=255)                                                                                                                                         cgGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGTTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCCGTGACCAATGGTCGGATCGTTGCCGTGAt          <  1:417208/148‑1 (MQ=255)                                                                                                                                          ggAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc      >  2:667996/1‑151 (MQ=255)                                                                                                                                             aaaaTTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt   >  1:253075/1‑151 (MQ=255)                                                                                                                                              aaaTTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc  <  2:947281/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508519‑3508809

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