Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,941	0	A	C	44.3%	104.4 / 10.6	61	E392D (GAA→GAC)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (28/6);  new base C (0/27);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 7.09e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAAT  >  NC_000913/3026791‑3027089                                                                                                                                                        |                                                                                                                                                      aGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGaa                                                                                                                                                       <  2:723878/151‑1 (MQ=255) aGCCTACAAAGTATGGCAATTACAATGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGaa                                                                                                                                                       <  1:698303/151‑1 (MQ=255)     tACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAAcc                                                                                                                                                     >  2:1016307/1‑149 (MQ=255)         aaGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGc                                                                                                                                              >  2:721529/1‑151 (MQ=255)            tATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCa                                                                                                                                            <  1:500474/150‑1 (MQ=255)            tATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCAc                                                                                                                                           >  2:730066/1‑151 (MQ=255)             aTTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACt                                                                                                                                          >  1:1224395/1‑151 (MQ=255)                 cAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGc                                                                                                                                      >  1:951086/1‑151 (MQ=255)                            cTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTg                                                                                                                           >  2:253092/1‑151 (MQ=255)                              aTCCCCTCGCGGCATATGAAGCTTTTTACCTGGCCACGCCCGGCGGAGCAAAATCTCTGGGCCTTGACGATTTTATTGGCAAC‑TTTTTACCTGGCAAAGGGGCTGTTTTGGGGGTGATGAACCCC‑GCCGCCACTCCGCTACAGCAGCTg                                                                                                                           <  2:681349/149‑1 (MQ=255)                                cGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCt                                                                                                                       >  1:1066070/1‑151 (MQ=255)                                  cccccgggCATATGAACGGTTTTACCTGGCCCCCCCGGGGGGGGGGAAATCCCGGGGCCTTGACGTTTTGATTGGCCAC‑TTTTTCCCGGGAAAAGAGGCTGTTTTCGGGGTGATGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTa                                                                                                                      <  1:362472/144‑1 (MQ=255)                                  ccTCTCGGCATATGAGGCTTTTTCCCTGCCCACCCTCGCCGGAGCGAATTCTCTGGCCCTTGACGTTTTGTTTGCCAAT‑TTTTTCCCTGGAAAAGAGGCTTTTTTCTTGGTGTTGGACCCC‑ACCGCCACTCCGCTACAGCAGCTgcgc                                                                                                                        <  2:1092367/148‑1 (MQ=255)                                  ccTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCGGGGCCTTGACGATTTGATTGGCAAC‑TTTTTCCCTGGAAAAGGGGCTGTTTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTAt                                                                                                                     <  1:206361/151‑1 (MQ=255)                                   ctctCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATg                                                                                                                    >  1:396951/1‑151 (MQ=255)                                    tctcGGCATATGATGCGTTTTACGTGCCCCCGCTCGGCGGGCCGAAGTCTCTGGGCTTTAGCGATTTTATTGGCACC‑TTTTTACCTGGCAAAGAGGGTGTTTGGGTGTTGATGGACCCC‑CCCGCCACTCCGCTACAGCAGCTGCGCTa                                                                                                                      <  2:671457/148‑1 (MQ=255)                                        ggCATATGAAGCGTTTTCCCTGGCCGCGCGGGGCGGGGCAAAATCTCGGGGCTTTGACGTTTTGATGGGAAAT‑TTTTTACCTGGAAAAGAGGCTGATTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAc                                                                                                               <  1:612365/151‑1 (MQ=255)                                            tatGAAGCGTTTTCCCGGGCCACGCTGGGGGGAGCGAAATCTCTGGGCCTGGACGATTTGATGGGAAAC‑TTTTTACCTGGAAAAGAGGCGGATTTCGGGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActct                                                                                                            <  1:353495/150‑1 (MQ=255)                                              tGAAGCTTTTTCCCTGGCCCCGCTCGGGGGAGGGAAATCTCGGGGCTTGGACGTTTTATTTGGCAAT‑TTTTTCCCTGGCAAGGGGGTTGATTTCGTGGTGTTGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGt                                                                                                          <  1:717018/150‑1 (MQ=255)                                               gAAGCGTTTTACCTGGCCACGCTGGGGGGAGCGAAATCTCGGGGCCTTGACGATTTGATTGGCAAT‑TTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGtt                                                                                                         <  1:505588/150‑1 (MQ=255)                                                aaGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTc                                                                                                       >  2:530037/1‑151 (MQ=255)                                                aaGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCG‑CACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCt                                                                                                      >  2:908486/1‑151 (MQ=255)                                                    gTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                                    >  2:38133/1‑150 (MQ=255)                                                     ttttCCCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCTTTGACGATTTGATTGGCAAC‑TTTTTCCCTGGAAAAGAGGCTGATTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTa                                                                                                   <  2:753663/150‑1 (MQ=255)                                                     ttttACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                                    >  1:937831/1‑149 (MQ=255)                                                     tttacccTGGCCACCCTGGGCGAAGCAAACTCTCTGGGCCGTGACGATTTAATTGGCACC‑TTTTTCCCTGGAAAAGAGGTTTTTTTCGTGGTGTTGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                                    <  1:1016203/144‑1 (MQ=255)                                                      tttCCCGGGCCACGCTCGGCGGGGCAAAATCTCTGGGCTTTGACGTTTTGATTGGCACC‑TTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAg                                                                                                  <  1:932547/150‑1 (MQ=255)                                                      tttACCGGCCCCCGCTCGGCGGAGCAAAATCTCTGGGCCTTGACGTTTTGATTGGCAAC‑TTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGAGGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAg                                                                                                  <  1:566555/150‑1 (MQ=255)                                                       ttcccTGGCCCCGCTCGGGGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTCCCTGGCAAAGGGGCTGTTTTCGTGGTGATGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt                                                                                                 <  2:529854/147‑1 (MQ=255)                                                            tGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACa                                                                                            >  1:89206/1‑150 (MQ=255)                                                               ccccGCTCGGCGGAGCGAATTCTCTGGGCCTTGACGATTTGATTGGCACC‑TTTTTCCCGGGCAAAGAGGCTGTTTTGGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACa                                                                                            <  2:753782/144‑1 (MQ=255)                                                                   gCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc                                                                                    >  1:1215952/1‑151 (MQ=255)                                                                   gCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc                                                                                    >  1:130816/1‑151 (MQ=255)                                                                          ggAGCGAAATCTCTGGGCTTTGACTATTTGATTGGCATC‑TTTTCACCTGCCAAAGAGTCTGATTTGGTGTTGATGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatga                                                                             <  2:255487/151‑1 (MQ=255)                                                                             gCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatg                                                                           >  2:106346/1‑150 (MQ=255)                                                                               aaattcccggggCCTTGACGATTTGATTGGCAAT‑TTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACg                                                                        <  1:788722/142‑1 (MQ=255)                                                                                aaaTCTCGGGGCTTTGACGATTTGATTGGCAAT‑TTTTTCCCTGGCAAGGAGGCTGATTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatg                                                                           <  2:85248/147‑1 (MQ=255)                                                                                 aaTCTCTGGGCTTTGACGATTTGTTTGGCAAT‑TTTTTCCCTGGCAAAGAGGCTGTTTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGtt                                                                      <  2:444160/151‑1 (MQ=255)                                                                                    ctctGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTggg                                                                   >  1:998474/1‑151 (MQ=255)                                                                                    ctctGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTggg                                                                   >  1:431777/1‑151 (MQ=255)                                                                                            ctggagattttgtttgGCAAC‑TTTTTCCCGGGCAAAGGGGCTGATTTGGTGGT‑ATTGAACCCAACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGAc                                                            <  1:572289/142‑1 (MQ=255)                                                                                                 cGATTTGATTGGCAAC‑TTTTTCCCTGGCAAAGAGGCTGATTTCGGGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTc                                                       <  2:1146989/150‑1 (MQ=255)                                                                                                   aTTTGATTGGCACC‑TTTTAACCGGGAAAAGAGGCTTATTTCGTGGTGATGAAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCg                                                      <  2:772108/149‑1 (MQ=255)                                                                                                   aTTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGAt                                                    >  2:497298/1‑151 (MQ=255)                                                                                                     ttttttGGCAAATTTTTTCCCTGGCAAAGAGGCTGATTTTGGGTTATTGCACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGa                                                     <  1:266258/145‑1 (MQ=255)                                                                                                       gATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTAc                                                >  2:941883/1‑151 (MQ=255)                                                                                                             cAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCAcc                                          >  2:837437/1‑151 (MQ=255)                                                                                                              aaC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCt                                         >  1:962689/1‑151 (MQ=255)                                                                                                                c‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTAc                                       >  1:50458/1‑151 (MQ=255)                                                                                                                  tttttACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACg                                      >  1:513094/1‑151 (MQ=255)                                                                                                                   ttttCCCGGGAAAAGAGGCTGTTTTCGTGGTGATGAACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACg                                      <  1:511365/150‑1 (MQ=255)                                                                                                                        ccTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATg                                >  2:3145/1‑151 (MQ=255)                                                                                                                        ccTGGAAAAGGGGCTGTTTTCGGGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGAt                                 <  1:730001/150‑1 (MQ=255)                                                                                                                           ggCAAAGAGCCTGATTTCGTGTTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgt                              <  1:90830/150‑1 (MQ=255)                                                                                                                             cAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgt                           >  1:739656/1‑151 (MQ=255)                                                                                                                               aaGAGGCTGATTTCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCt                         <  2:739720/151‑1 (MQ=255)                                                                                                                                      tGATTTCGTGGTGAGGCACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgt                      <  1:816214/147‑1 (MQ=255)                                                                                                                                         tttCGTGGTGATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACg                 <  1:1231410/149‑1 (MQ=255)                                                                                                                                         tttCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGa                >  1:372535/1‑150 (MQ=255)                                                                                                                                               ggtgATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCaa          <  1:346421/150‑1 (MQ=255)                                                                                                                                               ggtgATGGACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAAc         <  1:192920/151‑1 (MQ=255)                                                                                                                                               ggtgATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAAc         >  1:970496/1‑151 (MQ=255)                                                                                                                                                     ggACCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACta       <  1:735229/147‑1 (MQ=255)                                                                                                                                                      aaaCCC‑GCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaata    <  1:721464/148‑1 (MQ=255)                                                                                                                                                       aaCCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaataat  >  1:412839/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                      AGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCC‑ACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAAT  >  NC_000913/3026791‑3027089

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