Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,177	0	A	C	43.4%	125.2 / 13.8	77	G308G (GGT→GGG)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base A (36/7);  new base C (0/33);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 5.84e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGT  >  NC_000913/2921028‑2921316                                                                                                                                                        |                                                                                                                                            aGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑Tcccac                                                                                                                                            <  2:315909/151‑1 (MQ=255)    gcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑Tcccaccca                                                                                                                                         <  1:49722/151‑1 (MQ=255)    gcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑Tcccaccc                                                                                                                                          >  1:678145/1‑150 (MQ=255)    gcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑Tcccaccc                                                                                                                                          >  1:633321/1‑150 (MQ=255)      cATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGc                                                                                                                                       >  2:1344700/1‑151 (MQ=255)        tAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGcac                                                                                                                                     >  2:965448/1‑151 (MQ=255)           aGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACAcc                                                                                                                                  >  1:827603/1‑151 (MQ=255)           aGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACAcc                                                                                                                                  >  1:461995/1‑151 (MQ=255)              cGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACAccgcc                                                                                                                               <  1:814307/151‑1 (MQ=255)              cGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACAccgcc                                                                                                                               >  2:189102/1‑151 (MQ=255)               gatgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACAccgc                                                                                                                                >  2:1052317/1‑149 (MQ=255)                   atcatcCACTCAACGTTGACGTAGTTGCAGACCCCCATCCCCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATaa                                                                                                                          <  1:522087/151‑1 (MQ=255)                      atcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAaaac                                                                                                                       >  2:1319192/1‑151 (MQ=255)                      atcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAaaac                                                                                                                       >  2:1441223/1‑151 (MQ=255)                              aaCTTTGACGTAGTTGCAAACCCCCACCCCCACCGACACCAACATGCCCATTACGATCGGTTTTTTACGCGCAATGTC‑AAGCGATCCCTGGCGGCGCACCCGCCAGTCGGAAATAA‑CCCCCCCCACCCCACCGCCGATAaaaccgcaa                                                                                                                  <  1:45862/148‑1 (MQ=255)                              aaCGTGGACGTAGTTGCAAAACCCCATCCCCATCGACAGCAACATGCCCATTCCGATCGGTGTTTTACGCGCAATTTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGAAAAAAA‑CCCCCCCCAGCACACCGCCGATAaaaccgcaaacc                                                                                                               <  2:1198694/151‑1 (MQ=255)                                cGTTGACGTAGTTGCAGAACACCATCCCCATCGACAGCAACATGCCCATTACGATCGTTTTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGAAAATAA‑CCCCACCCAGCACACCGCCGATAaaaccgcaaaccg                                                                                                              <  2:731436/150‑1 (MQ=255)                                cGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAaaaccgcaaaccgc                                                                                                             >  1:184070/1‑151 (MQ=255)                                   tGACGTAGTTGCAAAACACCATCCCCATCGACAGAAACATGCCCATTACGATCGGTGTTTTCCGCGAAATTTT‑CAGCGATCCCGTGCGGCGCACCAGCCAGTCGGAAATAA‑CCCCCCCCAGCACACCGCCGATAaaaccgcaaaccg                                                                                                              <  1:295996/147‑1 (MQ=255)                                    gACGTGGTGGCAGAACCCCACCCCCACCGACAGCAACATGCCCATTACGATCGGTTTTTTACGCGCAATGTT‑CAGCGACCCCGGGCGGCGCATCACCCAGTCGAAAAAAA‑CCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcgg                                                                                                          <  1:1227458/150‑1 (MQ=255)                                      cGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaac                                                                                                       >  1:1365680/1‑151 (MQ=255)                                       gtttttGCAGAACACCATCACCAGCGACAGCAAAATGCCCATTATTGTCGGTTTTTTGCGCGAAATGTT‑CAGCTATCCCGTGCGGGGAACCACCCGGTCAGAAGTAT‑TCCCCCCCGCCACCCCGCCGATAAAACCGCAAACCGCcgg                                                                                                          <  2:905051/143‑1 (MQ=255)                                       gTAGTTGCAGAACCCCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGAAATGTT‑CAGCGACCCCGGGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCtggaac                                                                                                       <  2:642873/150‑1 (MQ=255)                                          gTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgga                                                                                                    >  1:637595/1‑150 (MQ=255)                                           ttaaaGACCACAATCCCCACCGACAGCAACACGCCCATTGCGAGCGGTGTTTTACGCGCAATGTC‑CAGGGCCCCCGGGGGGCGCTGCGCCGAGGGGGAAAAAA‑CCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacg                                                                                                      <  2:551509/143‑1 (MQ=255)                                           ttGCAGAACCCCATCCCCATCGACAGCAACATCCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGACCCCGTGCGGCGCATCAGCCAGTCGAAAATAA‑CCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacggaa                                                                                                   <  2:263738/150‑1 (MQ=255)                                            tGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGc                                                                                                 >  2:621257/1‑151 (MQ=255)                                            tGCAAAACCCCACCCCCACCGACAGCAACATGCCCTTTACGATCGGTTTTTTACGCGAAATGTC‑CAGCGACCCCGGGCGGCGCATCAGCCAGTCGAAAATAA‑CCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacggaa                                                                                                   <  2:59071/149‑1 (MQ=255)                                               ggAAGACCATCACCATCGCCAGCAACATGCCCATTACGATCGGCGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAC‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcac                                                                                              >  1:440418/2‑151 (MQ=255)                                               aGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcac                                                                                              >  1:868337/1‑151 (MQ=255)                                                         ccccatcGACAGCAACAGGCCCTTTACGATGGGTGTTTTCCGCGAAATGTT‑CAGCGACCCCGTGCGGCGCATCAGCAAGTCGAAAAAAC‑CCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGc                                                                                     <  2:366252/148‑1 (MQ=255)                                                         ccccagcGACAGCAACATCCCCTTTCCGATCGGTTTTTTACGCGAAAGTTC‑CAGCGACCCGGGGCGGCGCACCAGCCAGTGGAAAAAAA‑CCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCt                                                                                    <  2:1336363/149‑1 (MQ=255)                                                           ccaccGACAGAAACATCCCCATTCCGAGGGGTTTTTCCCGCGAAGTGTT‑CAGCGCCCCCGGGCGGCCCACCAGCCGGGCGAAAAAAA‑CCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGcc                                                                                        <  1:1121471/145‑1 (MQ=255)                                                           caacggACACCAACAGGCCCATTCCGATGGTTGTTTTCCCCGAAATGTT‑CAGCGACCCGGGGCGGCGCACCGGCCAGTCGAAAATAA‑CCCCCCCCAGCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCtt                                                                                   <  1:1408208/145‑1 (MQ=255)                                                             atcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                 >  2:761411/1‑150 (MQ=255)                                                               cGACAGAAACATGCCCATTCCGACCGGTGTTTTACGCGCATTGTT‑CAGCACCCCCGGGCGGCGCACCAGCCAGTCGAAAATAA‑CCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCa                                                                                <  2:490071/149‑1 (MQ=255)                                                                gACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGa                                                                              >  1:1426385/1‑150 (MQ=255)                                                                   acaaaCATCCCCATACCGAGCGGTTTTTTCCG‑GCAATTTTCCGGCGGCCCCGGGGGGCCCACCGGCCAGAGGAAAATAA‑CCCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCttt                                                                                  <  1:808200/140‑1 (MQ=255)                                                                     cAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc                                                                        >  2:808856/1‑151 (MQ=255)                                                                     cAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc                                                                        >  2:299221/1‑151 (MQ=255)                                                                      aaCATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCGCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc                                                                        >  2:478427/1‑150 (MQ=255)                                                                           cccaaTTCCGATCGGTGTTTTACGCGAAAGGTG‑CAACGTCCCCGTGCGGCACCCCAGCAAGTGGAAAAAAC‑CCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcg                                                                  <  1:194277/147‑1 (MQ=255)                                                                            cccttcaCGAGCGGTGTTTTACCCGCAATGTT‑CAGCGACCCGGTGCGGCGCATCCCCCAGTCGGAAAAAA‑CCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAGCGGAAGCCACAAAGCCCGCTTTCAGAATCg                                                                          <  1:321077/136‑1 (MQ=255)                                                                             catttACGATCGGTTTTTTCCCCGAAATTTC‑CACCGCCCCCGGGCGGCGCACCAGCCGGCCGAAAAAAA‑CCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgcg                                                                <  1:1169030/148‑1 (MQ=255)                                                                               aTTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTg                                                              >  2:1141742/1‑151 (MQ=255)                                                                                ttACGAGGGGTTTTTTCCGCGCAATTTC‑CAGCGACCCCGGGCGGCGCATCCGCCAGTGGAAAAAAC‑CCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTg                                                              <  1:825061/150‑1 (MQ=255)                                                                                    gATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGc                                                         >  1:10390/1‑151 (MQ=255)                                                                                     aTCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGc                                                         >  2:774627/1‑150 (MQ=255)                                                                                      tCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAc                                                       >  1:214845/1‑151 (MQ=255)                                                                                         gtgtTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                      >  2:615327/1‑149 (MQ=255)                                                                                         gtgaTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAACCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                      >  1:970981/1‑149 (MQ=255)                                                                                          tgtTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                      >  1:1384287/1‑148 (MQ=255)                                                                                              ttACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaa                                                >  2:279881/1‑150 (MQ=255)                                                                                              ttACGCGAAATTTC‑CACCGACCCCGGGCGGCGCACAACCCAGTCGAAAATAA‑CCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaa                                                <  2:1222745/150‑1 (MQ=255)                                                                                              ttACGCGAAATGTT‑CAGCGACCCCGTGCGCCGCACCAGCCAGTCGAAAATAA‑CCCCACCCACCACACCGCCAATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaa                                                <  1:852874/150‑1 (MQ=255)                                                                                               tACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAc                                              >  2:357228/1‑151 (MQ=255)                                                                                                 cgcgCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCg                                            >  2:1276240/1‑151 (MQ=255)                                                                                                   cgcAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCggg                                          >  1:1378075/1‑151 (MQ=255)                                                                                                    gcAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCggg                                          >  1:784528/1‑150 (MQ=255)                                                                                                      attGTA‑CAGCAATCCCGGGGGGCGCATCAGCCAGCCGAAAATAA‑TCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaaa                                               <  2:744566/141‑1 (MQ=255)                                                                                                          tt‑CAGCGAGCCGGTGCGGCTCATCAGCCAGTCGGAAAAAA‑GCCCACCCAGCACACCCCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAg                                    <  1:589639/150‑1 (MQ=255)                                                                                                                cGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaat                               >  2:1000166/1‑150 (MQ=255)                                                                                                                 gtccccGGGCGGCGCATCAGCCAGTCGAAAATAA‑CCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGta                                 <  2:1261916/144‑1 (MQ=255)                                                                                                                 gATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaata                              >  1:1258115/1‑150 (MQ=255)                                                                                                                  accccGTGCGGCGCATCACCCAGTCGAAAATAA‑CCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaa                                <  2:762399/145‑1 (MQ=255)                                                                                                                     cggggcggcgCATCACCCAGCCGAAAATAA‑CCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaa                         <  2:441317/147‑1 (MQ=255)                                                                                                                     ccGGGCGGCGCACCAGCCAGTCGAAAAAAA‑CCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaa                         <  2:139220/151‑1 (MQ=255)                                                                                                                      cGTGCGGCGAATCAGCCAGCGGAAAATAA‑CCCCCCCAAACACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaag                        <  2:172920/151‑1 (MQ=255)                                                                                                                       gTGCGGCGCATCATCAAG‑CGGAAAAAACCCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagt                       <  2:857047/151‑1 (MQ=255)                                                                                                                       gTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaag                        >  2:121126/1‑150 (MQ=255)                                                                                                                        tGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagt                       >  2:1139786/1‑150 (MQ=255)                                                                                                                         gcgtcgCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCACCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagta                      <  1:1126378/150‑1 (MQ=255)                                                                                                                           ggcgCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaagt                   >  1:693344/1‑151 (MQ=255)                                                                                                                           ggcgCATCAGCCAGTCGAAAATAA‑CCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaag                    <  2:1259647/150‑1 (MQ=255)                                                                                                                               cagcAGCGAGCCGGAAAAAT‑CCCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaagt                   <  2:1011634/144‑1 (MQ=255)                                                                                                                               cATCAGCCAGTCGAAAATAA‑CCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaagt                   <  2:911519/147‑1 (MQ=255)                                                                                                                                  cAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCg            >  1:1232058/1‑151 (MQ=255)                                                                                                                                  cAGCCAGTCGGAAAAAA‑CCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCg            <  1:1441183/151‑1 (MQ=255)                                                                                                                                     ccAGTCGAAAAAAA‑CCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGtt          <  2:387631/150‑1 (MQ=255)                                                                                                                                           ggAAATAA‑ACCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATa      <  2:538312/148‑1 (MQ=255)                                                                                                                                           gaaaaaaa‑CCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAg   <  2:1021487/149‑1 (MQ=255)                                                                                                                                             aaaaaacccccccccaGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt  <  1:151002/143‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                            AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTT‑CAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAA‑TCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGT  >  NC_000913/2921028‑2921316

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