Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,581	0	A	C	43.1%	92.6 / 17.7	58	V314G (GTA→GGA)	thiP	thiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (27/6);  new base C (0/25);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 4.74e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTC  >  NC_000913/73432‑73730                                                                                                                                                       |                                                                                                                                                      tCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTa                                                                                                                                                       <  1:956509/150‑1 (MQ=255) tCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAc                                                                                                                                                      >  1:61295/1‑151 (MQ=255)   tGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAccc                                                                                                                                                    >  1:842131/1‑151 (MQ=255)       gcgcCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATc                                                                                                                                                >  1:23047/1‑151 (MQ=255)          cccGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTGCGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATcgac                                                                                                                                             <  2:595660/151‑1 (MQ=255)                gTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCAc                                                                                                                                       >  1:325565/1‑151 (MQ=255)                gTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCAc                                                                                                                                       >  2:732186/1‑151 (MQ=255)                      gACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAg                                                                                                                                 <  2:493729/151‑1 (MQ=255)                                 agcagcATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGgc                                                                                                                       >  1:594311/1‑150 (MQ=255)                                       atcatGTCCAGCACCCCCCACAATCCCCCCGCCGCC‑GCCCAATTCCGAAAGGAGGTCCACCGCCCCTCCCACACCCCCGTTGGCGCCGCCACTCCCGCCAACTGGCGATTTACCCCACCGACGACCCCCGCCAGTAACGGTGGCAGCaa                                                                                                                  <  2:919931/149‑1 (MQ=255)                                       atcatGGCCAGCCCTCCGCACAACCCCCCGGCCGCCAGCGCAA‑TACGAAACGGGCCCCACAGCCCCTCCCACAGCCCCGGTTGTGCCAGCCCTCCCGGAACCGGGCGATTTCCCCCATCGCCAATCACCGCCAGTAACGGTGgcagca                                                                                                                   <  2:986151/148‑1 (MQ=255)                                        tcatGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAc                                                                                                                 >  1:17233/1‑149 (MQ=255)                                                     tACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGc                                                                                                  >  2:28803/1‑151 (MQ=255)                                                      aCGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGc                                                                                                  >  2:281768/1‑150 (MQ=255)                                                         cacaATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGcaca                                                                                               >  2:553681/1‑150 (MQ=255)                                                            aacacccccccccccAGCACAA‑TACACAACGAGCCCCACAGCGCCGCCCACAGCCCCGGTTGTGCCGGCTCTTCCGGCAACTGGCGATTTCCCCCACCGACAATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAtt                                                                                            <  2:628786/137‑1 (MQ=255)                                                            aaTCCACCTGCCGCAAGCGCAA‑TCCCCACCGAGCCCCACGCCCCCGGCCACACCACCGGTTGGGCCAGCATTTCCGGAAACGGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAtt                                                                                            <  1:28268/150‑1 (MQ=255)                                                              tACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTa                                                                                           >  2:720444/1‑149 (MQ=255)                                                              tACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAac                                                                                         >  2:626074/1‑151 (MQ=255)                                                               acacCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAac                                                                                         >  1:897761/1‑150 (MQ=255)                                                                cccccccccccAGCGCAA‑TACGCAACGGGCCCCACCGCGCCTGCCACAGCCCGGGTTGGGCCAGCACTTCCGGCAACGGGCGTTTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAac                                                                                         <  2:459466/140‑1 (MQ=255)                                                                       cgccAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCg                                                                                >  1:663953/1‑151 (MQ=255)                                                                        gccAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGCGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGc                                                                               >  2:493051/1‑151 (MQ=255)                                                                        cgccaggcaaaaacccAACGGGGCCCACACCCCCTCCCACAGCACCGTTTGTGCCGCCACTCCCGGAACCGGGCTATTTCCCCCACCGACAATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCg                                                                                <  2:1026845/136‑1 (MQ=255)                                                                         ccggcgcAA‑ACCACAACGGGCCCCACCGCCCCCCCCACACCCCGGGTTGTGCCGCCACTCCCGGCAACGGGCGATTTCCCCCACCGCCGATCCCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACAcc                                                                                      <  2:686648/140‑1 (MQ=255)                                                                            gcgcAA‑TCCCCAACGAGCCCCACCGCCCCCCCCACGCCCCCGGTTGGCCCAGCACTTCCGGCACCGGGGGTTTTCCCCCATCGACATTCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCaa                                                                             <  2:244284/149‑1 (MQ=255)                                                                               gcaaaccGCAACGGGGTCCAAAGCCCCTGCCACACCACCGGTGGTGCCAGCCCTCCCGGCAACTGGGGATTTCCCCCATCGACGAACACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCaaa                                                                            <  1:997449/142‑1 (MQ=255)                                                                                aa‑TACGCAACGAGGCCCACAGCGCCTGCCACAGCCCGGGTGGGGCCAGCACTCCCGGCAACGGGCGATTTCCCCCACCGCCGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGc                                                                         <  1:169098/149‑1 (MQ=255)                                                                                    ccGCACCGAGGCCCACGCCGCCTGCCACAGCCCCGGTTGGCCCAGCCCTCCCGGCATCTGGCGTTTTCCCCCATCGACGATCTCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCt                                                                     <  1:659665/149‑1 (MQ=255)                                                                                       cAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGc                                                                 >  2:465176/1‑151 (MQ=255)                                                                                           ggggcccacgcccccTCCCACAGCCCCGGTTGGGCCAGCACTCCCGGCAACGGGCGATTCCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGa                                                              <  2:247557/137‑1 (MQ=255)                                                                                             ggcccAGAGCGCCTGCCACAGCCCGGTTTGCGCCAGCACTTCCGGCAACGGGCGATTTCCCCCACCGACGTTCTCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGa                                                           <  2:194807/148‑1 (MQ=255)                                                                                                ccccggccccTGCCACAGCCCGGTTTGTGCCAGCACTTCCGACAGTTGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGa                                                           <  2:963839/140‑1 (MQ=255)                                                                                                ccACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCg                                                        >  1:574888/1‑151 (MQ=255)                                                                                                   cggcgcCTCCCACAGCACCGGTTGGGCCAGCCCTCCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTcc                                                     <  1:567293/149‑1 (MQ=255)                                                                                                     gccccTGCCACAGCACCGGTTGGGCCAGCTCTCCCGGCAACTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCg                                                    <  1:1061638/147‑1 (MQ=255)                                                                                                      cgcCTGCCACAGCCCCGGTTGCGCCAGCCCTCCCGGCAACTGGAGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCg                                                    <  2:891717/149‑1 (MQ=255)                                                                                                       gcCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGt                                                 >  1:608890/1‑151 (MQ=255)                                                                                                         cTGCCCCAGCCCCGTTTGTGCCAGCACTTCCGGCAACTGGCGATTCCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCggg                                                  <  1:908124/148‑1 (MQ=255)                                                                                                            ccACAGCCCGGTTTGCGCCAGCATTCCCGGCAACGGGCGATTTCCCCCAGCGACGAACACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTcc                                                     <  2:33313/142‑1 (MQ=255)                                                                                                                  cccgggtggggCCACCACTCCCGGCACGGGGGGATTTCCCCCATCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGGTATGCAGACGATCGTCCGGGTCGCGCCa                                          <  1:161922/137‑1 (MQ=255)                                                                                                                   aCCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCtt                                     >  2:232904/1‑151 (MQ=255)                                                                                                                         tgcgCCAGCCCTCCCGCCAACTGGCGTTTTCCCCCATCGCCGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGc                                        <  1:664789/139‑1 (MQ=255)                                                                                                                           tgCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCagcagc                             >  1:494403/1‑151 (MQ=255)                                                                                                                            gCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGCGGCAGCAACAGCAGCAGCGCCAGCACACTTAACACCGGGTCGCAAATGCGACTAAGCAGACGATCGTCCGGGTCGCGCCAGCCTTgcagcaac                             >  2:445103/1‑148 (MQ=255)                                                                                                                             ccaccccTTCCGGCAACTGGGATTTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgcagca                               <  1:858667/141‑1 (MQ=255)                                                                                                                              cAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtg                          >  2:136526/1‑151 (MQ=255)                                                                                                                               aGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtgg                         >  1:373419/1‑151 (MQ=255)                                                                                                                                gCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtggt                        >  1:58265/1‑151 (MQ=255)                                                                                                                                   cTTCCGGAAACTGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtg                          <  1:661213/146‑1 (MQ=255)                                                                                                                                     tCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCgg                   >  2:1002135/1‑151 (MQ=255)                                                                                                                                      ccGGCACCTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCggg                  <  1:980578/151‑1 (MQ=255)                                                                                                                                        ggCAACTGGCGATTTCCCCCATCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCggg                  <  2:663935/149‑1 (MQ=255)                                                                                                                                        ggCAACGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgc                 <  2:494389/150‑1 (MQ=255)                                                                                                                                         ggaccgggCGATTTCCCCCTTCGACAATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGcc                     <  1:748701/139‑1 (MQ=255)                                                                                                                                               ggCGTTTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGc         <  1:606893/151‑1 (MQ=255)                                                                                                                                                gCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGcc        >  1:800779/1‑151 (MQ=255)                                                                                                                                                   aTTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCtt      >  1:183964/1‑150 (MQ=255)                                                                                                                                                    tttACCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATgg          <  1:1054363/145‑1 (MQ=255)                                                                                                                                                       aCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTc  >  1:101419/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      TCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAA‑TACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTC  >  NC_000913/73432‑73730

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