BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,458,272	(G)_5→4	coding (218/1311 nt)	b2342 ←	putative acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,458,268	0	G	.	99.4%	803.1 / ‑2.2	177	coding (222/1311 nt)	b2342	putative acyltransferase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base . (91/84); total (91/86)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

Only 100 of 210 total aligned reads displayed.

GCCTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTA  >  NC_000913/2458176‑2458356

gCCTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTg                                                                                          >  2:4750556/1‑93 (MQ=255)

  cttgAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggg                                                                                        >  2:6231161/4‑93 (MQ=255)

  cTGGAAACTGGTAGCGCAAGCGCGGCTGAAGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggg                                                                                        <  2:3586676/93‑1 (MQ=255)

   tGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTgggg                                                                                       >  1:119185/1‑93 (MQ=255)

     gAAACTGTTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggct                                                                                     <  1:1325859/93‑3 (MQ=255)

     gAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggct                                                                                     >  2:5424990/1‑91 (MQ=255)

     gAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggct                                                                                     <  1:2297576/93‑3 (MQ=255)

      aaaCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggctt                                                                                    <  2:453969/93‑4 (MQ=255)

      aaaCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggctt                                                                                    >  1:3071331/1‑90 (MQ=255)

      aaaCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggctt                                                                                    <  1:4248229/93‑4 (MQ=255)

       aaCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggcttc                                                                                   >  1:290546/1‑89 (MQ=255)

       aaCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggcttc                                                                                   >  2:3107617/1‑89 (MQ=255)

       aaCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTggggcttc                                                                                   <  1:4357107/93‑5 (MQ=255)

        aCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCa                                                                                 <  2:4963876/93‑1 (MQ=255)

        aCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCa                                                                                 <  2:40374/93‑1 (MQ=255)

        aCTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCa                                                                                 <  2:2013215/93‑1 (MQ=255)

         cTGGTAGCTCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAg                                                                                <  2:2462484/93‑1 (MQ=255)

         cTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAg                                                                                >  2:3172795/1‑93 (MQ=255)

         cTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAg                                                                                >  1:5205459/1‑93 (MQ=255)

           ggTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGc                                                                              >  2:828506/1‑93 (MQ=255)

            gTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTGCCTAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCa                                                                             <  1:6353157/93‑1 (MQ=255)

            gTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCa                                                                             >  2:338674/1‑93 (MQ=255)

              aGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCAtt                                                                           <  2:13032/93‑1 (MQ=255)

                cgcAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTg                                                                         >  1:4161851/1‑93 (MQ=255)

                  cAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTa                                                                       <  2:6690884/93‑1 (MQ=255)

                   acgcgcgcCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTAc                                                                      <  2:5264832/91‑1 (MQ=255)

                   aaGCGGGGCTGACTCTGTAAGCATCGGTATGTACTTTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTAc                                                                      <  2:1647569/93‑1 (MQ=255)

                   aaGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTAc                                                                      <  1:1820719/93‑1 (MQ=255)

                    aGCGCGGCTGACGCTGTAAGCATCGGGATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACg                                                                     >  2:1697583/1‑93 (MQ=255)

                     gcgcGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGa                                                                    >  1:2161479/1‑93 (MQ=255)

                      cgcgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGAc                                                                   <  2:5760450/93‑1 (MQ=255)

                      cgcgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGAc                                                                   <  2:5536299/93‑1 (MQ=255)

                       gcgGCTGACGCTGTAAGCATCGGTATTTACATTCATTCCCGTACCGAGAACAATTTCTCGCGCAATGTTGGGGG‑TTCAGGCATTTGTACGAcc                                                                  <  2:794520/93‑1 (MQ=255)

                       gcgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGAcc                                                                  <  1:1022092/93‑1 (MQ=255)

                        cgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCt                                                                 <  2:4147424/93‑1 (MQ=255)

                        cgGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCt                                                                 >  2:4167429/1‑93 (MQ=255)

                                gCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAg                                                         <  2:2034057/93‑1 (MQ=255)

                                gCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAg                                                         >  1:3448210/1‑93 (MQ=255)

                                gCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAg                                                         <  1:5119734/93‑1 (MQ=255)

                                   gTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGAcc                                                      >  2:5974528/1‑93 (MQ=255)

                                    tAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCa                                                     >  2:4933524/1‑93 (MQ=255)

                                      aGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgt                                                   <  1:4935360/93‑1 (MQ=255)

                                      aGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgt                                                   >  1:3928167/1‑93 (MQ=255)

                                         aTCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgttgt                                                >  1:4520766/1‑93 (MQ=255)

                                          tCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAgttgtt                                               >  2:6140412/1‑93 (MQ=255)

                                           cGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTc                                              <  2:2643938/93‑1 (MQ=255)

                                           cGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTc                                              <  1:1591663/93‑1 (MQ=255)

                                            ggTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCa                                             >  1:3512400/1‑93 (MQ=255)

                                             gTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATG‑TAGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCaa                                            <  1:5414633/93‑1 (MQ=255)

                                              tATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAAt                                           <  1:5800213/93‑1 (MQ=255)

                                               aTGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATc                                          <  1:629488/93‑1 (MQ=255)

                                               aTGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATc                                          >  2:2545738/1‑93 (MQ=255)

                                               aTGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATc                                          <  1:6344536/93‑1 (MQ=255)

                                                tGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCa                                         >  1:4566539/1‑93 (MQ=255)

                                                 gTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCAc                                        <  1:6593602/93‑1 (MQ=255)

                                                 gTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCAc                                        <  2:6457378/93‑1 (MQ=255)

                                                  tACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACt                                       <  1:2439318/93‑1 (MQ=255)

                                                   aCATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACtt                                      <  2:1267030/93‑1 (MQ=255)

                                                    cattcattCCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTc                                     <  2:5621077/93‑1 (MQ=255)

                                                     attcattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcg                                    <  2:1774279/93‑1 (MQ=255)

                                                     attcattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcg                                    >  1:800687/1‑93 (MQ=255)

                                                      ttcattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcgg                                   >  1:4893502/1‑93 (MQ=255)

                                                      ttcattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcgg                                   <  2:2191187/93‑1 (MQ=255)

                                                       tcattcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTcggc                                  >  2:4428706/1‑93 (MQ=255)

                                                          ttcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCggg                               >  1:4436685/1‑93 (MQ=255)

                                                           tcCCGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCgggg                              <  2:3849658/93‑1 (MQ=255)

                                                            cccGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGa                             <  1:2659124/93‑1 (MQ=255)

                                                            cccGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGa                             <  2:4245209/93‑1 (MQ=255)

                                                              cGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtc                           <  1:3590346/93‑1 (MQ=255)

                                                              cGTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtc                           >  1:4673596/1‑93 (MQ=255)

                                                               gTACCGAGAACAATTTCACGCGCAATGTT‑TGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGGTGTTCAATCACTTCGGCGGGGAtct                          >  1:4051723/1‑93 (MQ=255)

                                                               gTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtct                          >  1:3277473/1‑93 (MQ=255)

                                                               gTACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtct                          <  2:5458081/93‑1 (MQ=255)

                                                                tACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtctc                         >  1:4188045/1‑93 (MQ=255)

                                                                tACCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGAtctc                         >  1:6730408/1‑93 (MQ=255)

                                                                 aCCGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCg                        >  2:1965461/1‑93 (MQ=255)

                                                                  ccGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGc                       >  1:5854727/1‑93 (MQ=255)

                                                                  ccGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGc                       >  1:957833/1‑93 (MQ=255)

                                                                   cGAGAACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCt                      <  2:5213948/93‑1 (MQ=255)

                                                                    gagaACAATATCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTg                     >  2:6619017/1‑93 (MQ=255)

                                                                     agaACAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGc                    >  2:5880835/1‑93 (MQ=255)

                                                                      ggaCAATTTCACGCTC‑ATGGTTGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCgcggcg                   <  1:5611377/91‑1 (MQ=255)

                                                                        aCAATTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGAGCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTg                 >  2:4684709/1‑93 (MQ=255)

                                                                          aaTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGcc               >  1:5505390/1‑93 (MQ=255)

                                                                          aaTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGcc               <  1:6117898/93‑1 (MQ=255)

                                                                          aaTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGcc               <  1:1253061/93‑1 (MQ=255)

                                                                          aaTTTCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGcc               <  2:6051988/93‑1 (MQ=255)

                                                                              tCACGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCcagca           <  2:4250881/93‑1 (MQ=255)

                                                                                aCGCGCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGt         <  1:2725685/93‑1 (MQ=255)

                                                                                 cgcgCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGtt        <  2:1667292/93‑1 (MQ=255)

                                                                                 cgcgCAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGtt        <  1:6164242/93‑1 (MQ=255)

                                                                                  gcgcAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTc       <  2:2966249/93‑1 (MQ=255)

                                                                                  gcgcAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTc       <  2:5134823/93‑1 (MQ=255)

                                                                                    gcAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGc     >  2:4332404/1‑93 (MQ=255)

                                                                                    gcAATGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGc     >  2:3701953/1‑93 (MQ=255)

                                                                                     cAATGTT‑GGGGCTTCAGTCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGcc    <  1:6429249/93‑1 (MQ=255)

                                                                                      aaTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCt   >  1:698579/1‑93 (MQ=255)

                                                                                      aaTGTT‑GGGGCTTCAGGCAATTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCt   >  1:6517533/1‑93 (MQ=255)

                                                                                       aTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTa  >  1:82290/1‑93 (MQ=255)

                                                                                       aTGTT‑GGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTa  >  1:1917156/1‑93 (MQ=255)

GCCTGGAAACTGGTAGCGCAAGCGCGGCTGACGCTGTAAGCATCGGTATGTACATTCATTCCCGTACCGAGAACAATTTCACGCGCAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTA  >  NC_000913/2458176‑2458356

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 19 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`