Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	44.1%	86.9 / 23.6	59	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/5);  new base C (0/26);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 6.86e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACG  >  NC_000913/2984293‑2984590                                                                                                                                                      |                                                                                                                                                      aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACAcc                                                                                                                                                     <  1:163184/151‑1 (MQ=255) aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACAcc                                                                                                                                                     >  2:567722/1‑151 (MQ=255)    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAGATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTg                                                                                                                                                  <  2:820608/151‑1 (MQ=255)     tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTgg                                                                                                                                                 >  1:147755/1‑151 (MQ=255)      gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTgg                                                                                                                                                 <  1:589551/150‑1 (MQ=255)           gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGccccc                                                                                                                                            <  2:835011/150‑1 (MQ=255)             gaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCccgcc                                                                                                                                         <  2:702161/151‑1 (MQ=255)                   aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATAc                                                                                                                                    >  1:176570/1‑150 (MQ=255)                   aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATAc                                                                                                                                    >  2:552108/1‑150 (MQ=255)                     aaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACaa                                                                                                                                  >  2:330205/1‑150 (MQ=255)                      aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACaa                                                                                                                                  >  2:1111858/1‑149 (MQ=255)                            ggacccaaaaaaaTGTATGCCGCCTTGCCAGCCATCATTATTCAACAAATAAATGCCGGCAGATTTTTAGCGGTAGGGAGAAAATGCAGATCTTTTCGTCACTTTAAGTGGTA‑AATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGc                                                                                                                          <  1:495908/142‑1 (MQ=255)                                ccaaaaaaaTGTATGCCGCCTTGCCACCCATCATTATTCAACAGATTAAGGACGGAAGATATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTAAAGGGTC‑AATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAg                                                                                                                       <  2:374374/145‑1 (MQ=255)                                 ccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                     >  2:945844/1‑151 (MQ=255)                                   aaaaaaTGTATGTCGTCTCCCCAGCCATCATTATTCAAAAGATAAAGGCCGGAAAATATTTAGCGGAAAGGAGAAAATGAAGAGCTATTCGCCCTGTTAAATGGTC‑AATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                   <  1:613443/149‑1 (MQ=255)                                       aaTGTATGCCGACTGGCCAGCCATCATTATTCAAAAGATAAATGACGGCAGATATTTATGGGAAAGGAGGAAATGCAGAGTTTTTCGTCACGTTCAATGGTA‑ACTCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                               <  1:1040238/151‑1 (MQ=255)                                        atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                              >  1:967055/1‑151 (MQ=255)                                         tgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                             >  2:95751/1‑151 (MQ=255)                                              gtCGACTTCGCAGCCATTATTATTAAAAAGATTAATGCGGGAAGATTTGTCTGGGTAAGGAAAAAATCAAGAGTTTTTCCCCAATTAAAGGGGTC‑AGTGCCCCCCCCTGACCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                        <  1:1189338/151‑1 (MQ=255)                                                     tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                                 >  1:788242/1‑151 (MQ=255)                                                     tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                                 >  1:1100587/1‑150 (MQ=255)                                                       ccagccaTCATTATTCAAAAAATTAAGGCCGGCAGATTTTTATGGGAAAGAAGAAAATGAAGAGTTTTTCGTCACTTTAAATGGTC‑AATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                       <  1:640629/143‑1 (MQ=255)                                                            catcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                          >  2:358757/1‑151 (MQ=255)                                                              tttttattaAACAATTAAAGGACGACATATATTTTTGGGAAGGGAGAAAATGCAGAGTTTTTCGCCAGTTTAAGTGGTC‑AAGCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                         <  1:1058523/147‑1 (MQ=255)                                                                attattCAACAGATTAATGACAGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                      >  2:240282/1‑151 (MQ=255)                                                                 ttattCAACAGATTAAGGACGGAAGATATTTATGGGTAAGGAGAAAATGAAGAGTTTTTCGTCACGTTAAATGGTA‑AATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                      <  2:147791/150‑1 (MQ=255)                                                                 ttattCAACAGATAAATGGCGGCAGATATTTTCCGGTAGGGAGGAAAGGCGGACCTATTCGTCCGGTTCAGTGGAC‑AACCCCACCCCCCGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                      <  2:55017/150‑1 (MQ=255)                                                                 ttattCAACAAATTAATGACGGCAATTTTTTATGGGAAAGGAGAAAATGCAGAGTTTTTCGTCACGTTCAATGGTC‑AATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                        <  1:53706/148‑1 (MQ=255)                                                                 ttataCAACAGATTAATAGCGGCATATATTTGGCGGAAAGGAGGAAATGAAGAGCTTTTCGCCACGTTCAATGTTA‑ACTCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                      <  1:1070273/150‑1 (MQ=255)                                                                 ttataCAACAGATAAATGGCGGCAGATATTTTGCGGTAAGGAGAAAATGCAGGGTTTTTCGTCACTTTAAAGGGTC‑AATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                      <  1:444342/150‑1 (MQ=255)                                                                   attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                   >  1:806617/1‑151 (MQ=255)                                                                      cAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                  >  2:478394/1‑149 (MQ=255)                                                                       aaCAAATAAATAACGGTAGATATTTACCGGAAGGGGAAAAATGCAGAGTTTTTCGTCATTTAAAATGGTC‑ACCCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                                <  1:1027329/150‑1 (MQ=255)                                                                        aCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                                >  2:222019/1‑149 (MQ=255)                                                                        aCAGATTAATGACGGCAGATATTTAGCGGAAGGGAGAAAATGAAGAGCTATTCGCCGTGTAAAGGGGTA‑ATTCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                 <  1:833462/148‑1 (MQ=255)                                                                          aGATTAAGGACGGCAAATATTTACCGGTAAGGAGAAAATGCAGAGCTATTCTTCACTTTCAATGGTA‑AATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                                                              <  2:411890/149‑1 (MQ=255)                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                          >  2:865514/1‑151 (MQ=255)                                                                               aaggACGGAAGATATTTATCGGAAAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAAGGTTC‑ATTCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                        <  2:1076100/147‑1 (MQ=255)                                                                               aaTGACGGCAGATTTTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACTTTCAATGGTC‑AATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                         <  1:404964/149‑1 (MQ=255)                                                                                aTGACGACAAATTTTTATGGGAAAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAAGGTTA‑AATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                        <  2:86867/149‑1 (MQ=255)                                                                                  gACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                                    >  1:276606/1‑151 (MQ=255)                                                                                  gACGGAAGATATTTATCGGAAAGGAGGAAATGCAGAGTTTTTCGTCCCGTTAAATGGCA‑AATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                                    <  1:1190804/151‑1 (MQ=255)                                                                                             tttATCGGTAGGGAGGAAATGCAGAGTTTTTCGCCACGTTAAATGGTC‑AATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                           <  1:567582/149‑1 (MQ=255)                                                                                                aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                      >  1:19326/1‑119 (MQ=255)                                                                                                aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                        >  1:904316/1‑149 (MQ=255)                                                                                                   ggTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                                   >  2:217458/1‑151 (MQ=255)                                                                                                        ggagaaAATGCAGAGTTTTTCTTCACTTTCAAGGGTC‑AATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                               <  2:752230/150‑1 (MQ=255)                                                                                                             aaaTGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                         >  1:635787/1‑151 (MQ=255)                                                                                                             aaaGGCGGAGCTTTTCGTCCCGTTAAATGGTA‑AACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                            <  1:999226/148‑1 (MQ=255)                                                                                                              aaTGCAGACCTATTCTACACTTTC‑ATGGCCAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa                                        <  2:788619/151‑1 (MQ=255)                                                                                                                     agCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                                 >  2:756591/1‑151 (MQ=255)                                                                                                                     actttttCGTCACGTAAAATGTTC‑AATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                      <  1:707165/141‑1 (MQ=255)                                                                                                                        ttttCGTCCCGTTAAAGGGAC‑ATTCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                                 <  1:749705/146‑1 (MQ=255)                                                                                                                         aTTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                             >  1:639626/1‑151 (MQ=255)                                                                                                                           tCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                           >  2:627908/1‑151 (MQ=255)                                                                                                                           tCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                           >  2:136554/1‑151 (MQ=255)                                                                                                                               cACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                       >  1:953035/1‑151 (MQ=255)                                                                                                                                  gTTAAATGGTC‑AACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccac                     <  2:788413/150‑1 (MQ=255)                                                                                                                                               aTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCattgacattgac        >  2:823641/1‑151 (MQ=255)                                                                                                                                                   cacaCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg     >  2:984090/1‑150 (MQ=255)                                                                                                                                                    acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc   >  1:890644/1‑151 (MQ=255)                                                                                                                                                     cacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACg  >  1:463863/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACG  >  NC_000913/2984293‑2984590

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