BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,722,918	T→A	N247Y (AAC→TAC)	b1647 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,722,918	0	T	A	99.5%	762.7 / NA	210	N247Y (AAC→TAC)	b1647	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); major base A (96/112); minor base C (1/0); total (98/112)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

TCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACA  >  NC_000913/1722830‑1723007

tCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAacc                                                                                       >  2:4014348/1‑93 (MQ=255)

 cTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAacca                                                                                      >  2:1147180/1‑93 (MQ=255)

  tGCTTCGACAGCTGCCCGTTCGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccac                                                                                     >  1:6540757/1‑93 (MQ=255)

   gCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTCAaccacc                                                                                    >  1:1640326/1‑93 (MQ=255)

   gCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                    >  1:1570523/1‑93 (MQ=255)

    cTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCt                                                                                   <  1:19596/93‑1 (MQ=255)

    cTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCt                                                                                   <  2:3842309/93‑1 (MQ=255)

     ttCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAATCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTg                                                                                  >  1:4700887/1‑93 (MQ=255)

     ttCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTg                                                                                  <  1:694483/93‑1 (MQ=255)

     ttCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTg                                                                                  <  2:3634686/93‑1 (MQ=255)

       cGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGtt                                                                                <  2:1533364/93‑1 (MQ=255)

         aCAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCa                                                                              <  1:595176/93‑1 (MQ=255)

          cAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCaa                                                                             >  1:1230179/1‑93 (MQ=255)

          cAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCaa                                                                             <  1:3185563/93‑1 (MQ=255)

           aGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGCTCCAt                                                                            >  1:324963/1‑93 (MQ=255)

            gCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATc                                                                           <  1:2875337/93‑1 (MQ=255)

              tGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                                         >  2:5635330/1‑93 (MQ=255)

              tGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                                         <  2:1065716/93‑1 (MQ=255)

               gCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAg                                                                        >  2:3393213/1‑93 (MQ=255)

                 ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                                      <  2:1580172/93‑1 (MQ=255)

                 ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                                      >  2:5501740/1‑93 (MQ=255)

                  cGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGccc                                                                     >  2:4258424/1‑93 (MQ=255)

                   gTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcg                                                                    <  1:3126830/93‑1 (MQ=255)

                   gTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcg                                                                    >  2:5271890/1‑93 (MQ=255)

                    tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcgc                                                                   <  1:5665119/93‑1 (MQ=255)

                       gcgcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGtt                                                                <  1:910055/93‑1 (MQ=255)

                        cgcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGttt                                                               <  1:5995580/93‑1 (MQ=255)

                        cgcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGttt                                                               <  1:6307341/93‑1 (MQ=255)

                         gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                                              >  1:6532362/1‑93 (MQ=255)

                          cTCAATGTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTaa                                                             <  1:2791160/93‑1 (MQ=255)

                           tCAATGTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTACAATCGAGCCCGCGTTTAAc                                                            <  1:5154682/93‑1 (MQ=255)

                            cAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAct                                                           >  1:5961215/1‑93 (MQ=255)

                              aTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcc                                                         <  1:1774902/93‑1 (MQ=255)

                              aTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcc                                                         <  2:3070425/93‑1 (MQ=255)

                              aTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcc                                                         <  1:21283/93‑1 (MQ=255)

                               ttttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcct                                                        <  2:2543429/93‑1 (MQ=255)

                               ttttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcct                                                        <  2:5892188/93‑1 (MQ=255)

                               ttttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcct                                                        <  1:106297/93‑1 (MQ=255)

                               ttttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcct                                                        >  2:3506588/1‑93 (MQ=255)

                               ttttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcct                                                        <  2:371890/93‑1 (MQ=255)

                               ttttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCCCGTTTAActcct                                                        >  2:868442/1‑93 (MQ=255)

                                tttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcctc                                                       <  2:1038436/93‑1 (MQ=255)

                                tttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcctc                                                       >  2:5409303/1‑93 (MQ=255)

                                 ttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATCCCCAGGCGTAACCCACCTGTTCAATCGAGCCCGCGTTTAACTCctct                                                      <  2:2174931/93‑1 (MQ=255)

                                   aCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGc                                                    >  2:6285908/1‑93 (MQ=255)

                                    ccTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTCAACCACCTGTTCAATCGAGCCCGCGTTTAGCTCCTCTGcc                                                   >  2:1757100/1‑93 (MQ=255)

                                     cTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGcca                                                  >  2:6121333/1‑93 (MQ=255)

                                     cTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGcca                                                  <  2:4712753/93‑1 (MQ=255)

                                      tGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccac                                                 >  2:1969632/1‑93 (MQ=255)

                                        aaCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacca                                               <  2:1218936/93‑1 (MQ=255)

                                          ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACcaca                                             <  2:3266766/93‑1 (MQ=255)

                                          ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACcaca                                             <  2:1840722/93‑1 (MQ=255)

                                            gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCCCAGc                                           >  2:2314966/1‑93 (MQ=255)

                                            gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGc                                           >  1:4654660/1‑93 (MQ=255)

                                            gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGc                                           <  2:1442965/93‑1 (MQ=255)

                                             aTAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGcc                                          >  2:4161853/1‑93 (MQ=255)

                                             aTAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGcc                                          <  1:5350227/93‑1 (MQ=255)

                                               aaTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAg                                        <  1:1647167/93‑1 (MQ=255)

                                               aaTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAg                                        <  1:2079477/93‑1 (MQ=255)

                                                aTTGGCAGCGGCTGCGATGGTAAACGTAATCCCCAGGCGTAAACCCCCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGt                                       <  1:3360889/93‑1 (MQ=255)

                                                aTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGt                                       <  2:6706767/93‑1 (MQ=255)

                                                 ttGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGtt                                      >  2:213427/1‑93 (MQ=255)

                                                 ttGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGtt                                      <  1:5509388/93‑1 (MQ=255)

                                                   ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtca                                    <  1:6405063/93‑1 (MQ=255)

                                                   ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtca                                    <  1:4170047/93‑1 (MQ=255)

                                                     cAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatc                                  <  2:4416844/93‑1 (MQ=255)

                                                      aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatca                                 >  2:4558686/1‑93 (MQ=255)

                                                      aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatca                                 <  2:2547708/93‑1 (MQ=255)

                                                        cGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACg                               <  1:3915438/93‑1 (MQ=255)

                                                        cGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACg                               >  1:3573355/1‑93 (MQ=255)

                                                         ggCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGc                              <  2:237953/93‑1 (MQ=255)

                                                          gCTGCGATGTTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCa                             <  1:4568435/93‑1 (MQ=255)

                                                          gCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCa                             <  2:2641950/93‑1 (MQ=255)

                                                          gCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCa                             >  1:4755723/1‑93 (MQ=255)

                                                            tGCGATGGTAAACGTAATCCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGc                           <  2:1712548/93‑1 (MQ=255)

                                                            tGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGc                           <  1:4010242/93‑1 (MQ=255)

                                                             gCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCg                          >  2:4504846/1‑93 (MQ=255)

                                                               gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGc                        >  1:3712649/1‑93 (MQ=255)

                                                                aTGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                       <  2:2843148/93‑1 (MQ=255)

                                                                 tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                      <  1:3042021/93‑1 (MQ=255)

                                                                 tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                      <  1:1986656/93‑1 (MQ=255)

                                                                 tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                      >  2:6397895/1‑93 (MQ=255)

                                                                        cGTAATGCCCAGGCGTAAACCCCCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataa               <  1:2552649/93‑1 (MQ=255)

                                                                         ctAATACCCAGGCTTAAACCACCTGTTCAATCGAGCCCGCGTTTACCTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataat              <  1:2811964/92‑1 (MQ=255)

                                                                           aaTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatca            <  2:5954846/93‑1 (MQ=255)

                                                                           aaTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatca            <  1:5807122/93‑1 (MQ=255)

                                                                            aTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcat           >  2:2707316/1‑93 (MQ=255)

                                                                            aTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcat           <  1:6190709/93‑1 (MQ=255)

                                                                              tcccAGGCGTAACCCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAAtcatca         <  2:6432015/92‑1 (MQ=255)

                                                                              aCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAAtcatca         >  1:4231808/1‑93 (MQ=255)

                                                                               cccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATcatcat        <  1:1660908/93‑1 (MQ=255)

                                                                                ccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTACCTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCAtt       <  1:1540451/93‑1 (MQ=255)

                                                                                 cAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTa      >  1:4600717/1‑93 (MQ=255)

                                                                                 cAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTa      >  1:4944973/1‑93 (MQ=255)

                                                                                  aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaa     >  1:2308222/1‑93 (MQ=255)

                                                                                   ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa    <  1:2802446/93‑1 (MQ=255)

                                                                                   ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa    >  2:3425959/1‑93 (MQ=255)

                                                                                   ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa    >  1:1948233/1‑93 (MQ=255)

                                                                                    gCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAAc   >  2:2773181/1‑93 (MQ=255)

                                                                                    gCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAAc   >  1:877575/1‑93 (MQ=255)

TCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACA  >  NC_000913/1722830‑1723007

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