BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,130,334	A→G	V30V (GTT→GTC)	b2983 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,130,334	0	A	G	99.7%	988.3 / NA	288	V30V (GTT→GTC)	b2983	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (135/151); total (135/153)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 326 total aligned reads displayed.

CGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACAACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGCCAA  >  NC_000913/3130244‑3130423

 gATTGCAGAATCACCAGTACGCAAAACAGGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCacgaca                                                                                        <  2:702484/93‑1 (MQ=255)

 gATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGGCCCCGGGGGCGCCAACAGACCCAGCAAGGCGCCCacgaca                                                                                        >  1:1947837/1‑93 (MQ=255)

  aTTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACAt                                                                                       >  2:532661/1‑93 (MQ=255)

  aTTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACAt                                                                                       >  2:4578212/1‑93 (MQ=255)

        caATCACCAGTACGCCAAACATGGCGGGCGTCCTCCCCTTACCCGCACACGAGAGCGGCAAACGACCCAGCAACGCGCTCACGAAATTACtag                                                                                 >  2:4601851/1‑91 (MQ=255)

        caATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTgg                                                                                 <  1:8701045/93‑1 (MQ=255)

        caATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTgg                                                                                 <  2:2981876/93‑1 (MQ=255)

          aTCACCAGTACGCCAAACAGGGTGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAg                                                                               <  1:605920/93‑1 (MQ=255)

            cACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGcc                                                                             <  1:8533486/93‑1 (MQ=255)

            cACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTAACGGCACACGAGAGCGCCAACAGACCCAGCAAGGCCCTCACGACATTACTGGAGcc                                                                             >  1:2570637/1‑93 (MQ=255)

             aCCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCCACAGACCCAGCCAGGCGCTCACGACATTACTGGAGCCt                                                                            >  2:3633627/1‑93 (MQ=255)

              ccAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTa                                                                           >  2:7631944/1‑93 (MQ=255)

                     gCCAAACATGGCGGGCGTCATCCCTTTACCGGCCCACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGgc                                                                    >  2:3030587/1‑93 (MQ=255)

                     gCCAAACATGGCGGGCGGCATCCCTTTAACGGGACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCAGGACAATACTGGGGCCCAACCAggg                                                                    >  2:5911047/1‑92 (MQ=255)

                      ccAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGTAGCCTAACCAGgcg                                                                   >  1:7489956/1‑93 (MQ=255)

                      ccAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGgcg                                                                   <  1:1276499/93‑1 (MQ=255)

                      ccAAACATGGCGGGCGTCATCCCTTTACCCGCACACGAGAGCGCCCACAGACCCGGCAAGGCGCTCACGACATTACTGGAGCCTAAAAAGgcg                                                                   >  2:1970643/1‑93 (MQ=255)

                        aaaCATGGCGGGGGTCATCCCTTTACCGGCACAAGAGAGGGCCAAAAGACCCAGCAAGGCGGTCATGACCTATCTGGCGCCTTAACAGGCGCt                                                                 >  2:2846025/1‑93 (MQ=255)

                        aaaCATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTTGAGCCTAACCAGGCGCt                                                                 <  1:8386451/93‑1 (MQ=255)

                         aaCATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTa                                                                >  2:8807528/1‑93 (MQ=255)

                         aaCATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTa                                                                <  1:5980577/93‑1 (MQ=255)

                         aaCATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTa                                                                >  2:3302652/1‑93 (MQ=255)

                         aaCATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTa                                                                <  2:5457506/93‑1 (MQ=255)

                         aaCATGGCGGGCGTCAGCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTa                                                                <  2:1347977/93‑1 (MQ=255)

                          aCATGGCGGGCGTCATCCCTTTCCCGCCACACGAGAGCGCCAACAGCCCCAGCAGGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTAt                                                               <  1:245371/93‑1 (MQ=255)

                          aCATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTAt                                                               >  1:2061438/1‑93 (MQ=255)

                          aCATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTAt                                                               >  1:7964578/1‑93 (MQ=255)

                           cATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTAtt                                                              <  2:6059731/93‑1 (MQ=255)

                           cATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCAAGACATTACTGGAGCCTAACCAGGCGCTAtt                                                              >  1:6171152/1‑93 (MQ=255)

                               gcgggcgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTCGAGCCTAACCAGGCGCTATTGCGa                                                          >  1:1844839/1‑93 (MQ=255)

                                ggggcgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaa                                                         <  1:7574839/92‑1 (MQ=255)

                                ggggcgCCACCCCTTCCCCGGCACACGAGACCCCAAACAGACCCAGCAGGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaa                                                         <  1:6439460/92‑1 (MQ=255)

                                cgggcgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaa                                                         >  1:458037/1‑93 (MQ=255)

                                    cgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataat                                                     <  1:2942464/93‑1 (MQ=255)

                                    cgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataat                                                     >  2:201877/1‑93 (MQ=255)

                                    cctCACCCCTTTACCGGCACACGAAAGCGCCAATAGCCCCAGCAAGGCCCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataat                                                     <  2:8509861/91‑1 (MQ=255)

                                     gTCATCCCTTTACCGGCACACGAGAGCGCCAACAGGCCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATg                                                    <  1:7177276/93‑1 (MQ=255)

                                     gTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATg                                                    <  1:5560494/93‑1 (MQ=255)

                                     gTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATg                                                    >  2:4803645/1‑93 (MQ=255)

                                      tCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATgg                                                   >  2:1687213/1‑93 (MQ=255)

                                       ccccccTTTCCCGGCACACGAGAGCGCCAACAGCCCCAGCAGGGCCCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGt                                                  <  1:3269402/90‑1 (MQ=255)

                                       cATCCCTTTCCGGGCACACGAGAGCCCAAACAGCCCCAGCAGGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGt                                                  <  2:8095946/93‑1 (MQ=255)

                                        aTCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTg                                                 <  2:2601171/93‑1 (MQ=255)

                                         tCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGc                                                >  2:7470241/1‑93 (MQ=255)

                                         tCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGc                                                >  2:3999998/1‑93 (MQ=255)

                                          cccTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCg                                               <  2:1023712/93‑1 (MQ=255)

                                          cccTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCg                                               >  1:5274070/1‑93 (MQ=255)

                                          cccTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCg                                               >  2:3251192/1‑93 (MQ=255)

                                           ccTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGa                                              >  2:2517431/1‑93 (MQ=255)

                                            cTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGaa                                             >  1:5697875/1‑93 (MQ=255)

                                            cTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGaa                                             >  2:267739/1‑93 (MQ=255)

                                             tttACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGaaa                                            >  1:1144729/1‑93 (MQ=255)

                                             tttACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGaaa                                            <  1:2836631/93‑1 (MQ=255)

                                              ttACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAAc                                           >  1:4097536/1‑93 (MQ=255)

                                                acgggCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAAcgc                                         <  1:6548631/90‑1 (MQ=255)

                                                 ccGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAAcgcg                                        >  1:7660945/1‑93 (MQ=255)

                                                 ccGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAAcgcg                                        >  1:5082850/1‑93 (MQ=255)

                                                   ggCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGcc                                      <  1:5469677/93‑1 (MQ=255)

                                                   ggCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGcc                                      >  1:5414655/1‑93 (MQ=255)

                                                   ggCACACCAGACCGTCAACAGACACAGCAACACGCTCACGACATTACTTGAGCCTCAGCAGGCGCTATTGCGAATAATCTTGCGAAACGCGcc                                      >  2:7064173/1‑93 (MQ=255)

                                                    gCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCa                                     <  1:7215429/93‑1 (MQ=255)

                                                     cacaCGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCAt                                    <  2:6946565/93‑1 (MQ=255)

                                                      acacGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATc                                   >  2:7925314/1‑93 (MQ=255)

                                                       cacGAGGCCGCCAACAGCCCCAGGAAGGCCCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCt                                  <  1:2273742/93‑1 (MQ=255)

                                                       cacGAGAGCGCCAACAGCCCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCt                                  <  1:3720333/93‑1 (MQ=255)

                                                       cacGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCt                                  <  2:4194965/93‑1 (MQ=255)

                                                       cacGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCt                                  <  1:8217714/93‑1 (MQ=255)

                                                          gagaGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCa                               <  1:8083392/93‑1 (MQ=255)

                                                             agCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAAc                            <  1:4224232/93‑1 (MQ=255)

                                                              gcgcCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAAcc                           <  1:7043808/93‑1 (MQ=255)

                                                                gcCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCaa                         <  1:786232/93‑1 (MQ=255)

                                                                 ccAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAAt                        <  1:2325793/93‑1 (MQ=255)

                                                                 ccAACAGACCCAGAAAGGCGCGCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAAt                        <  2:2418812/93‑1 (MQ=255)

                                                                  cAACAGCCCCAGCAAGCCCCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATg                       <  2:6826321/93‑1 (MQ=255)

                                                                  cAACAGCCCCAGCAAGACGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATg                       <  2:1314462/93‑1 (MQ=255)

                                                                  cAACAGACCCAGCACGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATg                       <  1:4747281/93‑1 (MQ=255)

                                                                    aCAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGtt                     <  2:5700113/93‑1 (MQ=255)

                                                                    aCAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGtt                     <  2:4648503/93‑1 (MQ=255)

                                                                     cAGCCCCAGCAGGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGttt                    <  1:3395981/93‑1 (MQ=255)

                                                                     cAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGttt                    <  1:6996899/93‑1 (MQ=255)

                                                                     cAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGttt                    <  2:2272990/93‑1 (MQ=255)

                                                                     cAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGttt                    <  2:1149671/93‑1 (MQ=255)

                                                                      aGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTg                   <  1:7343961/93‑1 (MQ=255)

                                                                       gACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGa                  <  1:1514731/93‑1 (MQ=255)

                                                                          ccAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGt               >  1:4550261/1‑93 (MQ=255)

                                                                          caaGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGt               <  1:3503760/91‑1 (MQ=255)

                                                                           cAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGtt              >  2:6189040/1‑93 (MQ=255)

                                                                            aGCAGGCCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTa             <  1:3158633/93‑1 (MQ=255)

                                                                            aGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTa             >  2:3460847/1‑93 (MQ=255)

                                                                             gCACGGCCCCCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaa            <  1:6340992/93‑1 (MQ=255)

                                                                             gCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGTAAACCAATGTTTGATGTTaa            >  2:6513484/1‑93 (MQ=255)

                                                                             gCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaa            <  2:3899762/93‑1 (MQ=255)

                                                                               aaGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaaaa          >  1:3617300/1‑93 (MQ=255)

                                                                                aGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAAc         >  2:1635252/1‑93 (MQ=255)

                                                                                aGGCGATCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAAc         >  2:8744121/1‑93 (MQ=255)

                                                                                 ggCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAAcc        >  1:3941936/1‑93 (MQ=255)

                                                                                  gcccTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACCCGCCATCTGCAAACCAATGTTTGATGTTAAAAccc       <  1:6845509/90‑1 (MQ=255)

                                                                                   cgcTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCg      <  2:8396466/93‑1 (MQ=255)

                                                                                    gcTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGc     >  1:8784257/1‑93 (MQ=255)

                                                                                     cTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGcc    <  1:4368867/93‑1 (MQ=255)

CGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACAACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGCCAA  >  NC_000913/3130244‑3130423

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