Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,413	0	A	C	40.4%	105.4 / 13.2	58	V370G (GTG→GGG)	thiP	thiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (30/4);  new base C (0/23);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 3.31e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTT  >  NC_000913/73265‑73563                                                                                                                                                      |                                                                                                                                                       aTAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCAcc                                                                                                                                                      >  1:532656/1‑151 (MQ=255)     gggATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGAc                                                                                                                                                  <  1:122318/151‑1 (MQ=255)           gCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGcc                                                                                                                                             <  2:661003/150‑1 (MQ=255)            ccATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGcc                                                                                                                                             >  1:898025/1‑149 (MQ=255)               ttAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCAt                                                                                                                                        >  2:919070/1‑151 (MQ=255)                 aaCGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCAttt                                                                                                                                      >  1:1091170/1‑151 (MQ=255)                 aaCGCATGGGTAAAAATCACAATGCCGCCAGCAGATGGTGGCAGGCCGATAGTGTTGTTGGGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCACCCCCGGCCCCGCCAGCAtt                                                                                                                                       <  2:135281/150‑1 (MQ=255)                   cGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCAtttt                                                                                                                                     >  2:428938/1‑150 (MQ=255)                    gCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCt                                                                                                                                   >  1:797897/1‑151 (MQ=255)                     cATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCt                                                                                                                                   >  1:429937/1‑150 (MQ=255)                         ggTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCcgc                                                                                                                              >  1:558753/1‑151 (MQ=255)                         ggTAAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCcgc                                                                                                                              <  1:396091/151‑1 (MQ=255)                             aaaaTCACAATGCGGTCAGCAGATGGTGGCAGGCCGATGGTGTTGTTGAGCA‑TAAAAAAAAGCCCGGTGGCCAGCACAACCCCCGGCATGGCGGGGATCAACCTGCCGCTCATCCCCAGCCCCCGACCCGCCAGCATTTTCTGCCGCGccc                                                                                                                          <  2:1131002/151‑1 (MQ=255)                             aaaaCCACAAGCCCGTCAGCAGATTGTGGCGGGCCGATGGTGTTGTGGAGCAGTAAAAAGAGG‑CCGGTAGCCAGCACAACCCCCGGCATGGCGAGTATCACCATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCcgcg                                                                                                                             <  2:154431/148‑1 (MQ=255)                                  cacaATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGt                                                                                                                     >  2:765663/1‑151 (MQ=255)                                     aaTGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcg                                                                                                                  >  1:793254/1‑151 (MQ=255)                                          cGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGAc                                                                                                              >  2:783765/1‑150 (MQ=255)                                                cagATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTcc                                                                                                        >  1:1058091/1‑150 (MQ=255)                                                 agATTGTGGCAGGCCGATAGTTTTGTTGAGCAGAAAAAAGAAG‑CCGGTAGCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCt                                                                                                          <  1:136218/147‑1 (MQ=255)                                                   aTTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCAt                                                                                                      >  1:812750/1‑149 (MQ=255)                                                    ttGTGGCGGGCCGATGGTGTTGTGGAGCAGAAAAAAGAAG‑CGGGTAGCCAGCACAATCCCCGGCTGGGCGAGGATCAACATGCCGCCCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATa                                                                                                     <  2:132659/149‑1 (MQ=255)                                                    ttGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCAt                                                                                                      <  1:866596/148‑1 (MQ=255)                                                         gCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcat                                                                                              >  2:593740/1‑151 (MQ=255)                                                          cAGGCCGATGGGTTTGTTGAGCAGAAAAAAGAAG‑CCGTTGGCCAGCACAACCCCCGGCAGGGCGAGGATAACCAGCCCGCTAATCTCCCCCCCCGGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcat                                                                                              <  2:970394/150‑1 (MQ=255)                                                             gcctatgttgttgtTGGGCAGAAAAAAGAAG‑CCGGTGGCCACCACATTCCCGGTCTTGGCGAGGATCAACATGCCGTACATCCCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATg                                                                                          <  1:122765/143‑1 (MQ=255)                                                              ccGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATgg                                                                                         >  1:541691/1‑151 (MQ=255)                                                                 aTAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTc                                                                                       >  1:512743/1‑150 (MQ=255)                                                                  tattgttggTGACCAGTAAAAAGAAG‑CCGGTAGCCCCCACAATCCCGGTCATGGCGAGGATCAACTTGCCGCTATTCACCAGCCCCTTACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCa                                                                                      <  1:231329/147‑1 (MQ=255)                                                                   aGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGc                                                                                    >  1:918654/1‑151 (MQ=255)                                                                   aGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGc                                                                                    >  1:1042302/1‑151 (MQ=255)                                                                    gtgtTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGc                                                                                    >  1:245612/1‑150 (MQ=255)                                                                       ttgttgAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGc                                                                                    <  1:566578/147‑1 (MQ=255)                                                                       ttgttgAGCAGAAAAAAGAAG‑CCGGTAGCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg                                                                                     <  1:417475/146‑1 (MQ=255)                                                                              gCAGTAAAAAGAAG‑CCGGTAGCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc                                                                             <  1:834928/147‑1 (MQ=255)                                                                                agaaaaaaGAGG‑CCGGTAGCCAGCACAACCCCCGGCAGGGCGAGGATCACCAGCCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACg                                                                              <  2:670929/141‑1 (MQ=255)                                                                                aGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAAt                                                                        >  2:461360/1‑150 (MQ=255)                                                                                   aaaaaGAGG‑CCGTTGGCCGCCACAACCCCCGGCATGGCGAGGACCACCTTGCCCCTCATCCCCACCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATa                                                                       <  1:301893/148‑1 (MQ=255)                                                                                    aaaaGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACAcc                                                                   >  1:856687/1‑151 (MQ=255)                                                                                         aag‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTgccg                                                              >  2:269650/1‑151 (MQ=255)                                                                                               ggTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAgc                                                         >  1:127592/1‑151 (MQ=255)                                                                                               ggTAGCCACCACAACCCCGGGCAGGGCGAGGACCAACATGCCGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAg                                                          <  2:1042267/150‑1 (MQ=255)                                                                                                gTACCCAGCACAACCCCCGGCAGGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAgcg                                                        <  2:566369/151‑1 (MQ=255)                                                                                                    ccAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAgcgc                                                       <  1:179804/148‑1 (MQ=255)                                                                                                          caATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCa                                               >  2:1025124/1‑150 (MQ=255)                                                                                                           aacccccGGCATGGCGAGGATCAACATGCCGCTCACCCCCAGCACCGGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAAc                                             <  1:681090/148‑1 (MQ=255)                                                                                                           aaTCCCCGGAATGGGGAGGATAAACAGCCCGCCCATCCCCACCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAAc                                             <  2:999939/151‑1 (MQ=255)                                                                                                              ccccGGCAGGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACg                                            <  1:654905/149‑1 (MQ=255)                                                                                                              ccccGGCAGGGCAAGGACCAACATCCTGCTCATCCCCACCCCCGGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAg                                          <  2:590129/151‑1 (MQ=255)                                                                                                                      tGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCcac                                    >  2:660059/1‑149 (MQ=255)                                                                                                                       ggCGAGGATCAACATGCCGCTCATCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCCGCCGCCAGCGCAATACGCAACGAGGTCcaca                                   <  1:653407/149‑1 (MQ=255)                                                                                                                       ggCGAGAATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCcaca                                   <  1:783969/149‑1 (MQ=255)                                                                                                                         cGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACACACAATACACCGGCCGCCAGCGCAATACGCAACGAGGTCCACAgcgc                               >  2:927758/1‑151 (MQ=255)                                                                                                                            ggATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCt                             <  1:144856/150‑1 (MQ=255)                                                                                                                                cAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCcac                        >  1:238551/1‑151 (MQ=255)                                                                                                                                    aTGCCGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCa                    <  1:314115/151‑1 (MQ=255)                                                                                                                                      gCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCAcc                  >  1:177460/1‑151 (MQ=255)                                                                                                                                           tCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTtg             >  1:430644/1‑151 (MQ=255)                                                                                                                                             aTCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTtgtg           >  1:71113/1‑151 (MQ=255)                                                                                                                                                      aCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACtt  >  1:33984/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       ATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAG‑CCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTT  >  NC_000913/73265‑73563

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