BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,998,290	+G	coding (379/486 nt)	b4566 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,998,289	1	.	G	99.8%	1862.9 / ‑4.1	538	coding (380/486 nt)	b4566	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (279/258); total (279/259)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

Only 100 of 592 total aligned reads displayed.

GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑GACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376

gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCa                                                                                          <  2:6410289/93‑1 (MQ=255)

gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCa                                                                                          <  2:2516155/93‑1 (MQ=255)

gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCAACATAAAACACGCCGCCa                                                                                          >  1:5119404/1‑93 (MQ=255)

 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑g                                                                                        >  1:965327/1‑93 (MQ=255)

 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑g                                                                                        >  1:286176/1‑93 (MQ=255)

   tGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCGCCCGCGACATAAAACACGCCGCCA‑gca                                                                                      <  1:6794816/93‑3 (MQ=255)

    gCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑ggac                                                                                     <  2:4068706/93‑4 (MQ=255)

       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGa                                                                                   >  2:686805/1‑93 (MQ=255)

       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGa                                                                                   <  2:5335721/93‑1 (MQ=255)

       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGa                                                                                   <  2:7734875/93‑1 (MQ=255)

       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGa                                                                                   >  1:3359082/1‑93 (MQ=255)

        cGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGAt                                                                                  >  2:7044743/1‑93 (MQ=255)

           aaaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGACCGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATaaa                                                                               <  2:4861266/93‑1 (MQ=255)

           aaaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACTATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATaaa                                                                               >  1:7739086/1‑93 (MQ=255)

           aaaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATaaa                                                                               <  1:1850838/93‑1 (MQ=255)

            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAg                                                                              >  2:6113844/1‑93 (MQ=255)

            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAg                                                                              >  2:3599757/1‑93 (MQ=255)

            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAg                                                                              >  2:3119896/1‑93 (MQ=255)

             aCTGAATCTCTTCACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGa                                                                             >  1:1981451/1‑93 (MQ=255)

              cTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGaa                                                                            >  2:8387096/1‑93 (MQ=255)

              cTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGaa                                                                            <  1:3859399/93‑1 (MQ=255)

               tGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGaaa                                                                           >  1:1245775/1‑93 (MQ=255)

               tGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGaaa                                                                           >  2:2930201/1‑93 (MQ=255)

                gAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAAt                                                                          >  1:4159417/1‑93 (MQ=255)

                gAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAAt                                                                          >  2:5492757/1‑93 (MQ=255)

                 aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAAtt                                                                         >  1:1340166/1‑93 (MQ=255)

                    ctctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTccc                                                                      >  1:7643244/1‑93 (MQ=255)

                    ctctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTccc                                                                      <  1:2417081/93‑1 (MQ=255)

                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCt                                                                     >  2:7107558/1‑93 (MQ=255)

                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCt                                                                     >  1:5673677/1‑93 (MQ=255)

                      ctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCtt                                                                    >  2:7909680/1‑93 (MQ=255)

                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTc                                                                   <  2:542949/93‑1 (MQ=255)

                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTc                                                                   >  2:2210851/1‑93 (MQ=255)

                        ttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCa                                                                  >  2:6623179/1‑93 (MQ=255)

                          aCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATTTCCCGCGACATAAAACACGCCGCCAGGACGGTTAAAGACATTCCCTTAAcc                                                                >  2:4549470/1‑93 (MQ=255)

                           cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCt                                                               >  2:7201104/1‑93 (MQ=255)

                            tGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTg                                                              <  2:827941/93‑1 (MQ=255)

                            tGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTg                                                              >  1:8007936/1‑93 (MQ=255)

                              gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGaa                                                            >  1:6555414/1‑93 (MQ=255)

                               ttttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAg                                                           <  2:5959177/93‑1 (MQ=255)

                                tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGc                                                          <  1:5266449/93‑1 (MQ=255)

                                 ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCg                                                         <  2:3742553/93‑1 (MQ=255)

                                 ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCg                                                         <  2:895284/93‑1 (MQ=255)

                                 ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCg                                                         >  2:8010729/1‑93 (MQ=255)

                                  tGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCga                                                        <  1:7112012/93‑1 (MQ=255)

                                   gAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCgag                                                       <  2:3837943/93‑1 (MQ=255)

                                    aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCgaga                                                      >  2:6219820/1‑93 (MQ=255)

                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 >  2:4255630/1‑93 (MQ=255)

                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  2:2718761/1‑93 (MQ=255)

                                           ccAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACg                                               >  1:8363/1‑93 (MQ=255)

                                             aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGTGAGATTTAACGcc                                             >  2:3517321/1‑93 (MQ=255)

                                             aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGcc                                             <  1:4002784/93‑1 (MQ=255)

                                              gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                                            >  2:5036657/1‑93 (MQ=255)

                                               acgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCGCCTGAAGCGAGATTTAACGCCAc                                           <  2:1137661/93‑1 (MQ=255)

                                                cgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAcc                                          <  2:577727/93‑1 (MQ=255)

                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        <  1:4649008/93‑1 (MQ=255)

                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  2:5805764/1‑93 (MQ=255)

                                                   acgaTGCCGACTTTATCTCCCGCGACATAAAACTCGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTc                                       <  1:8515259/93‑1 (MQ=255)

                                                    cgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTct                                      <  1:2776851/93‑1 (MQ=255)

                                                     gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGCAGCGAGATTTAACGCCACCTTctg                                     >  2:7332957/1‑93 (MQ=255)

                                                     gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctg                                     >  2:4629051/1‑93 (MQ=255)

                                                     gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctg                                     >  2:4441893/1‑93 (MQ=255)

                                                     gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctg                                     >  1:7393565/1‑93 (MQ=255)

                                                      aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgc                                    >  1:1755167/1‑93 (MQ=255)

                                                       tGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgct                                   >  1:3669262/1‑93 (MQ=255)

                                                         ccGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGa                                 <  1:62011/93‑1 (MQ=255)

                                                           gACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGatc                               >  2:175953/1‑91 (MQ=255)

                                                           gACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAcc                               >  1:5970099/1‑93 (MQ=255)

                                                           gACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAcc                               >  1:754626/1‑93 (MQ=255)

                                                             cTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGt                             >  1:100364/1‑93 (MQ=255)

                                                             cTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGt                             <  1:3477106/93‑1 (MQ=255)

                                                              tttATCTCCCGCGCCCTCAGACACACCGCCAGGACGGATAAAGAAATTCCCTTCACCTGATGCAAGATTTAACGCCACCTTCTGCTGACCGtt                            <  1:7828892/93‑1 (MQ=255)

                                                               ttATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGttt                           <  2:346889/93‑1 (MQ=255)

                                                               ttATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGttt                           <  1:2252262/93‑1 (MQ=255)

                                                                tATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          <  2:7463229/93‑1 (MQ=255)

                                                                tATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          <  2:570208/93‑1 (MQ=255)

                                                                 aTCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTc                         <  2:1198178/93‑1 (MQ=255)

                                                                 aTCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTc                         <  1:5854290/93‑1 (MQ=255)

                                                                  tctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCg                        >  2:6990632/1‑93 (MQ=255)

                                                                  tctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCg                        >  1:8658357/1‑93 (MQ=255)

                                                                   ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGa                       <  1:4319671/93‑1 (MQ=255)

                                                                   ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGa                       >  1:7454534/1‑93 (MQ=255)

                                                                    tcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt                      <  1:3074968/93‑1 (MQ=255)

                                                                     cccGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATg                     >  2:1590752/1‑93 (MQ=255)

                                                                     cccGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATg                     >  1:1552830/1‑93 (MQ=255)

                                                                      ccGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGt                    <  2:5156547/93‑1 (MQ=255)

                                                                      ccGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGt                    >  2:2116867/1‑93 (MQ=255)

                                                                        gcgACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCCCCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTAc                  <  2:1392007/93‑1 (MQ=255)

                                                                        gcgACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTAc                  >  1:5384065/1‑93 (MQ=255)

                                                                           aCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa               >  1:6969397/1‑93 (MQ=255)

                                                                            cATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc              <  2:8594215/93‑1 (MQ=255)

                                                                            cATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc              >  1:4814348/1‑93 (MQ=255)

                                                                             aTAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGCTTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc             >  1:4475477/1‑93 (MQ=255)

                                                                             aTAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc             <  2:3342013/93‑1 (MQ=255)

                                                                              tAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            >  1:1078276/1‑93 (MQ=255)

                                                                              tAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            >  1:5546671/1‑93 (MQ=255)

                                                                                 aaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgc         >  1:2118155/1‑93 (MQ=255)

                                                                                    acGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc      <  2:1881994/93‑1 (MQ=255)

                                                                                       ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa   >  1:1463261/1‑93 (MQ=255)

                                                                                        cgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAt  <  1:4959347/93‑1 (MQ=255)

GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑GACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376

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