Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,118	0	A	C	42.1%	99.0 / 11.9	59	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (23/10);  new base C (0/24);  total (23/36)
Fisher's exact test for biased strand distribution p-value = 2.86e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTT  >  NC_000913/1620971‑1621268                                                                                                                                                    |                                                                                                                                                       aCTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAccc                                                                                                                                                     >  1:8512/1‑151 (MQ=255)  cTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAcccc                                                                                                                                                    >  2:634132/1‑151 (MQ=255)       ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAG‑CCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCttt                                                                                                                                              >  1:15981/1‑151 (MQ=255)           gTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                            <  2:1033024/150‑1 (MQ=255)            tGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgc                                                                                                                                          >  1:528562/1‑151 (MQ=255)            tGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgc                                                                                                                                          >  1:522512/1‑151 (MQ=255)               gACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAc                                                                                                                                       >  2:28814/1‑151 (MQ=255)                    aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                                  <  1:232932/151‑1 (MQ=255)                     aaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTa                                                                                                                                 <  1:587088/151‑1 (MQ=255)                               ggAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTg                                                                                                                        >  1:346956/1‑150 (MQ=255)                                   aaGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAAttt                                                                                                                   <  1:112143/151‑1 (MQ=255)                                   aaGCGACCAGCATAAAGCGAAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAAttt                                                                                                                   <  2:533232/151‑1 (MQ=255)                                      cGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCAt                                                                                                                >  1:737793/1‑151 (MQ=255)                                        aCCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAACCCCCACTTTGATGACCCAAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTCCCGCCTCTGAATTTCAt                                                                                                                <  2:212833/149‑1 (MQ=255)                                         ccAGCATAAACCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCCCTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCt                                                                                                              <  2:775942/150‑1 (MQ=255)                                          cAGCATAAAGCGCAACCCGGCCAGCACCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATCCGCCCACTCCCAGTAGCGCCACCCCCCCACCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtt                                                                                                             <  2:528955/150‑1 (MQ=255)                                             cATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                         >  1:910410/1‑151 (MQ=255)                                              aTAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCCCTTTGATGACCACAAAATTTAGCCCCCATACGACCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTg                                                                                                          <  1:598921/149‑1 (MQ=255)                                                aaggccaaacccggccggCAACAGCGGTGCCATTTTATAAAGCCCCCCTTGTATGCCCAAAAAATTAAGCCCCCATCCGACCACTGCCATGAGCGCCAACCCCCCATCTTTCCCCGACATTCTACCGCCTCTGAATTTCATCTTTTg                                                                                                          <  1:641245/137‑1 (MQ=255)                                                aaaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCCCTTTGATGACCCCAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                         <  1:685080/148‑1 (MQ=255)                                                            ggCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAg                                                                                          <  1:103409/151‑1 (MQ=255)                                                            ggCCAGCATCAGCGGTGCCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGACCACTCCCAGTAGCTCCACCCCCCCATCTTTCCGCACCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                             <  2:574786/148‑1 (MQ=255)                                                             gCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                         >  2:814978/1‑151 (MQ=255)                                                              ccAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                                        >  1:1033708/1‑151 (MQ=255)                                                               cgcatccggCGGTGGCATATAATGAAGCCCGATTTGGATGGCCGCAAAATTGAGCCCCCATACGCCCACGACAAGGAGCTCCATCTCCCCACCTTTTCCCAACATTGGACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAg                                                                                          <  2:458315/140‑1 (MQ=255)                                                               cAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                       >  1:523608/1‑151 (MQ=255)                                                               cAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGACCCCTCCCAGTGGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAg                                                                                          <  1:955963/148‑1 (MQ=255)                                                                   accAGCGGTGGCAGGTTATGAGGCCCCTCTTTGATGCCCACAAAATTGAGCCCCCTTACGGCCATTACCAGTAGCGCCAACACCCCATCTTTTCGCGATATTCAACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAtt                                                                                   <  1:108885/149‑1 (MQ=255)                                                                   aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                                    <  1:170345/150‑1 (MQ=255)                                                                     cAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAttt                                                                                  >  2:594553/1‑150 (MQ=255)                                                                      aGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                                >  1:422800/1‑151 (MQ=255)                                                                       gCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACt                                                                               >  2:887565/1‑151 (MQ=255)                                                                         ggTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                             >  2:652867/1‑151 (MQ=255)                                                                         ggTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATACGCCCACTACCAGTAGCGCCACCACCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                              <  1:244725/150‑1 (MQ=255)                                                                            ggCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCt                                                                          <  1:246422/151‑1 (MQ=255)                                                                               aTGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                        >  2:383932/1‑150 (MQ=255)                                                                                 gTTATGAAGCCCCACTTTGATGCCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTaa                                                                      <  1:979209/150‑1 (MQ=255)                                                                                    aTGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGt                                                                  >  1:821902/1‑151 (MQ=255)                                                                                       aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAg                                                                   <  2:698522/147‑1 (MQ=255)                                                                                       aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCTCTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGt                                                                  <  1:1075829/148‑1 (MQ=255)                                                                                       aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                                <  1:490638/150‑1 (MQ=255)                                                                                          ccccATTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                             <  1:766538/150‑1 (MQ=255)                                                                                          ccccACTTTGATGACCACAAAATTTACCCCCCATACGCCCCCTCCCAGTAGCCCCAACCCCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                            <  2:820915/151‑1 (MQ=255)                                                                                           cccTCTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                           <  2:158241/151‑1 (MQ=255)                                                                                           cccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                           <  1:896537/151‑1 (MQ=255)                                                                                           cccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                           >  1:463128/1‑151 (MQ=255)                                                                                                tttGATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGtt                                                         <  1:694241/148‑1 (MQ=255)                                                                                                          aaaaaaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCCCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTa                                             <  1:285189/148‑1 (MQ=255)                                                                                                                    gCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                   >  2:326724/1‑150 (MQ=255)                                                                                                                     cccccATCCGACCGCTACCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAt                                    <  1:787386/148‑1 (MQ=255)                                                                                                                     cccccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTa                                  >  2:992104/1‑150 (MQ=255)                                                                                                                         cATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcg                             >  1:499402/1‑151 (MQ=255)                                                                                                                                  aCTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAg                     <  2:1010884/150‑1 (MQ=255)                                                                                                                                      ccAGTAGCGCCAACCCCCCATCTTTTCGCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc                   <  2:756345/148‑1 (MQ=255)                                                                                                                                      ccAGAAGCCCAAACGCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGc                 <  1:1026088/150‑1 (MQ=255)                                                                                                                                            gcgcCAACCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtata             <  2:924665/148‑1 (MQ=255)                                                                                                                                            gccccAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtata             <  2:66588/145‑1 (MQ=255)                                                                                                                                             cgcCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACtt          <  2:15980/150‑1 (MQ=255)                                                                                                                                               ccAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg       >  2:70684/1‑151 (MQ=255)                                                                                                                                                   cacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtt    <  1:687995/150‑1 (MQ=255)                                                                                                                                                    acCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtttt  >  2:869865/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                       ACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTT  >  NC_000913/1620971‑1621268

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