Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,941	0	A	C	42.1%	89.4 / 10.5	57	E392D (GAA→GAC)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (25/8);  new base C (0/24);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.55e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATA  >  NC_000913/3026791‑3027090                                                                                                                                                        |                                                                                                                                                      aGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGaa                                                                                                                                                       <  2:676648/151‑1 (MQ=255)   ccTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAAcc                                                                                                                                                     >  2:313466/1‑151 (MQ=255)              ttGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTc                                                                                                                                          >  1:684461/1‑151 (MQ=255)                gCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTcc                                                                                                                                         >  1:1026899/1‑150 (MQ=255)                          ggCTACCGCCTCTCGGCATATGAGGCTTTTTCCCGGCCCACGCTCGGCGGAGCAAAATCTCGGGGCCTTGACTATTTGATTGGCAC‑CTTTTTCCCTGGCAAAGAGGTTGATTTCGTGGTGATGAAACCCACCGCCACTCCGCTAcagcag                                                                                                                               <  1:251841/150‑1 (MQ=255)                            ctaccccctccggccATAGGAGGCTTTTTCCCTGCCCCCGCCCGGCGGAGCAAAATCTCGGGGCCTTGCCGTTTTGTTTGCCAC‑TTTTTCCCCTGCAAAAGAGGCTTTTTTCGTGGGGATGGACCCCACCGCCACTCCGCTAcagc                                                                                                                                 <  1:576656/132‑1 (MQ=255)                            cTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTg                                                                                                                            >  1:569565/1‑151 (MQ=255)                             cctCGCCTCTAGGCATACGAAGCGTTATACCTGGCCACGCTCGGCGTAGCGAAATCTCTGGGCCTTGAAGGTCTAATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTCCGTGGTGATGGAACCCACCGCCACACCGCTAcagca                                                                                                                                >  1:908237/3‑146 (MQ=255)                              aTCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTgcg                                                                                                                          >  1:910656/1‑151 (MQ=255)                              aTCGCCTCTCGGCAAAGGAAGCTTTTTCCCTGGCCCGGCCCGGGGGAGCGAAATCTCGGGGCTGGAACGTTTTGAGTGGAAA‑TTTTTTACGTGAAAAAGGGGCTGATTCGGTGGTGATGGACCCAACCGCCACTCCGCTACAGCAGCTgc                                                                                                                           <  2:766819/150‑1 (MQ=255)                                cGCCTCGCGGATTATAAAGCTTTTTACCTGGCCACGCTCGGGGGAGCAAAATCTCGGGCCCGTGACGATTTGATTGACAA‑CTTTTAACCTGGCAAGGAGGTTGATTTCGTGGGGATGGAACCCACCGCCACTCCGCTACAGCAGCTgc                                                                                                                           <  2:290997/148‑1 (MQ=255)                                 gCCTCCCGGCTTAGGAAGTGTTTTCCCGGCCCCCGCTCGGGGGGGCAAATTCTCGGGGCTTTGACGTTTTGTTGGGCAT‑TTTTTTCCCTGGAAAAGGGCCTATTTTCGTGGGGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCt                                                                                                                        <  1:912287/150‑1 (MQ=255)                                  ccTCTCGGCATATGAAGCTTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCTTTGACGATTTGATTGGCAA‑TTTTTTCCCGGGCAAAGAGGTTGATTTCGGGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTAt                                                                                                                      <  2:649474/151‑1 (MQ=255)                                  ccTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTAt                                                                                                                      >  1:475337/1‑151 (MQ=255)                                  ccTCGCGGCATATGAGGCGTTTTCCCGGGCCACGCTCGGCGGGGCAAAATCTCGGGGCTTTGACGTTTTGATGGGCAA‑TTTTTTCCCGGGAAAAGGGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTgc                                                                                                                           <  1:411730/146‑1 (MQ=255)                                        ggaatatGAAGCGTTTTACCTGGCCACGCGGGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTCCCTGGAAAAGAGGCTGATTTCGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAc                                                                                                                <  2:321357/148‑1 (MQ=255)                                        ggCATATGAAGCTTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAT‑CTTTTTACCGGGAAAAGAGGCTGATTTGGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATg                                                                                                                     <  2:715393/146‑1 (MQ=255)                                        gccATATGAGGCGTTTACCCTGGCCCCGTGCGGGGGAGAAAAATCCCTGGGCCTAGAGGATTTGTTGGGCAA‑CTTTTTACCTGGAAAAGAGGTTGATTTCGGGGGGGTGAACCCCGCCGCCACTCCGCTACAGCAGCTGCGCTATGa                                                                                                                    <  2:434247/145‑1 (MQ=255)                                         gCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAct                                                                                                               >  2:788456/1‑151 (MQ=255)                                         gCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAct                                                                                                               >  2:841971/1‑151 (MQ=255)                                            tatGAGGCTTTTTACCTGGCCACGCTCGGCGGAGCAAAATCTCGGGGCTTTGACGTTTTGATTGGCAT‑CTTTTTCCCTGGCAAAGGGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActct                                                                                                             <  1:863502/150‑1 (MQ=255)                                              tGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGtt                                                                                                          >  2:198642/1‑151 (MQ=255)                                               gAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGttt                                                                                                         >  1:185721/1‑151 (MQ=255)                                                 aGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAAAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCt                                                                                                       >  1:88582/1‑151 (MQ=255)                                                            tGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACa                                                                                             >  1:412270/1‑150 (MQ=255)                                                             ggCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACaa                                                                                            >  2:299900/1‑150 (MQ=255)                                                               ccccgccggggggAGCAAAACCTCGGGCCCTTGACGTTTTTA‑TGGCAACCTTTTTCCCTGGAAAAGGGTTTGTTTCGGGGGTTAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAtt                                                                                         <  1:494968/140‑1 (MQ=255)                                                               ccACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAtt                                                                                         >  2:107933/1‑151 (MQ=255)                                                                  cGCTCGGCGGAGCGAAATCTGGGGGCCTTGACGATTTGATGGGCAA‑CTTTTTACCTGGCAAAGAGGTTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGtt                                                                                      <  1:645195/151‑1 (MQ=255)                                                                   gCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc                                                                                     >  1:949408/1‑151 (MQ=255)                                                                    cTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGAAAAGAGGTTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGtt                                                                                      <  2:188449/149‑1 (MQ=255)                                                                       ggcggAGGAAAATGTCTGGGCGTTACCAATTCGAGTGACAC‑CTTTTTACCTGGCAAAGTGGCTGATTCCGTGGAGATGGCCCCCACCGCCACTCCGCTACAGCAGCTGCGTTATGACAACTCTGTTTCTTTAGTCGACaaa                                                                                           <  2:856753/141‑1 (MQ=255)                                                                       ggcggAGCAAAATCTCGGGGCTTTGACGATTTGATTGGCAA‑CTTTTTCCCGGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGtt                                                                                      <  1:225450/146‑1 (MQ=255)                                                                                    ctctGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTggg                                                                    >  2:382346/1‑151 (MQ=255)                                                                                        gggcgtagccgttttGATTCGAAC‑CTTTTCACCTGGCAAGGGGGGTGATTCCGTGGTGAAGCACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGc                                                                   <  2:317152/136‑1 (MQ=255)                                                                                        gggCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGAAAAGAGGCTGATTTCGTGTTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTgg                                                                     <  2:874951/146‑1 (MQ=255)                                                                                        gggCCTTGACGATTTGATTGGCAA‑CTTTTTACCGGGAAAAGGGGCTGTTTTCGTGGTGATGAAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTgg                                                                     <  2:799886/146‑1 (MQ=255)                                                                                            cTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGAcc                                                            >  2:1170595/1‑151 (MQ=255)                                                                                            cTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGAcc                                                            >  1:39144/1‑151 (MQ=255)                                                                                              tGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGt                                                          >  1:108333/1‑151 (MQ=255)                                                                                                aCGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTc                                                        >  2:86908/1‑151 (MQ=255)                                                                                                 cgttttGATTGGCAA‑TTTTTCACCTGGAAAAGAGGCTGATTGCGTGGTGATGGAACCAACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTc                                                        <  2:896147/147‑1 (MQ=255)                                                                                                    tttGATTGGCAA‑TTTTTTCCCGGGAAAAGGGGTTTATTTCGTGGTGATGGAACCCGCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCg                                                       <  2:539077/148‑1 (MQ=255)                                                                                                      tGATTGGCAA‑CTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATc                                                    <  2:670415/149‑1 (MQ=255)                                                                                                              aa‑CTTTTTCCCGGGAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCAc                                            <  1:20147/149‑1 (MQ=255)                                                                                                                   ttttACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGt                                      >  1:133001/1‑151 (MQ=255)                                                                                                                    tttccctgcaatagagGTTGTTTTGGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGt                                      <  1:769215/138‑1 (MQ=255)                                                                                                                     ttACCTGGAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGtt                                     <  2:412271/150‑1 (MQ=255)                                                                                                                       aCCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGAt                                  >  2:923513/1‑151 (MQ=255)                                                                                                                       aCCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGAt                                  >  2:961940/1‑151 (MQ=255)                                                                                                                           ggaaaaGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATgg                                <  2:1017551/146‑1 (MQ=255)                                                                                                                           gccAAAGAGGTTGTTTCTGTGGTATTGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtc                              <  1:1182264/149‑1 (MQ=255)                                                                                                                           gaaaaaGAGGTTTATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtc                              <  1:94191/148‑1 (MQ=255)                                                                                                                              aaaGAGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgtc                           <  1:391514/151‑1 (MQ=255)                                                                                                                              aaaGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGCCAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgt                            >  1:1176176/1‑150 (MQ=255)                                                                                                                                   ggcttttttCGGGGTGGTGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtg                      <  1:536050/145‑1 (MQ=255)                                                                                                                                    gCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtgt                     <  2:172200/151‑1 (MQ=255)                                                                                                                                                  gagggAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACta        <  1:381941/147‑1 (MQ=255)                                                                                                                                                        aCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaataata  >  2:53170/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                      AGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAA‑CTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATA  >  NC_000913/3026791‑3027090

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