BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,327,579	T→C	Y80C (TAC→TGC)	b2226 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,327,579	0	T	C	99.7%	1077.0 / ‑5.1	313	Y80C (TAC→TGC)	b2226	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (153/159); total (154/159)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

Only 100 of 327 total aligned reads displayed.

CCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGT  >  NC_000913/2327489‑2327668

ccACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAg                                                                                          >  1:6303819/1‑93 (MQ=255)

 cacaGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAgg                                                                                         <  2:7363376/93‑1 (MQ=255)

  acaGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGCTTCCGCGACACTGGCAggg                                                                                        >  2:517010/1‑93 (MQ=255)

   caGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgc                                                                                       <  1:5781985/93‑1 (MQ=255)

    aGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgct                                                                                      <  2:8401417/93‑1 (MQ=255)

    aGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGCGAGTTTCCGCGACACTAGCAGGgct                                                                                      >  2:4500507/1‑93 (MQ=255)

     ggAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgctg                                                                                     >  1:2272949/1‑93 (MQ=255)

      gAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgctgc                                                                                    >  2:6458551/1‑93 (MQ=255)

      gAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgctgc                                                                                    <  2:5998331/93‑1 (MQ=255)

      gAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgctgc                                                                                    <  1:7868766/93‑1 (MQ=255)

       aaCGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgctgct                                                                                   >  2:7907357/1‑93 (MQ=255)

         cGTCCCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑cc                                                                                <  1:8148267/93‑1 (MQ=255)

         cGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑cc                                                                                <  2:7495542/93‑1 (MQ=255)

             ccAGCGCGGTATCGCGCGTTATGCCCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGtt                                                                            <  1:6366271/93‑1 (MQ=255)

               aGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGa                                                                          >  2:1767032/1‑93 (MQ=255)

               aGCGCGGTAACGCGCGTTATGCTCACCCCAGCCTGACAACAATAAACCTCTTCCGGTGAGTTTCCGCGGCACTGGCAGGGCTGCT‑CCGGTTGa                                                                          <  2:1634834/93‑1 (MQ=255)

                gcgcGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGAT‑CCGGTTGAt                                                                         >  2:8309747/1‑93 (MQ=255)

                 cgcgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATg                                                                        <  2:2729277/93‑1 (MQ=255)

                 cgcgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATg                                                                        <  2:3447415/93‑1 (MQ=255)

                 cgcgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATg                                                                        <  1:912603/93‑1 (MQ=255)

                  gcgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATgg                                                                       >  1:7562671/1‑93 (MQ=255)

                  gcgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATgg                                                                       <  2:4302865/93‑1 (MQ=255)

                  gcgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATgg                                                                       >  2:43297/1‑93 (MQ=255)

                   cgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATggg                                                                      >  1:7090414/1‑93 (MQ=255)

                     gTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTa                                                                    <  2:4612215/93‑1 (MQ=255)

                     gTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTa                                                                    <  1:80319/93‑1 (MQ=255)

                         cgcgcgTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTAtcct                                                                >  2:3951292/1‑93 (MQ=255)

                          gcgcgTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTAtcctc                                                               >  2:1233049/1‑93 (MQ=255)

                          gcgcgTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTAtcctc                                                               <  2:5836011/93‑1 (MQ=255)

                          gcgcgTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTAtcctc                                                               <  2:5923563/93‑1 (MQ=255)

                            gcgTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTcca                                                             >  2:3376938/1‑93 (MQ=255)

                             cgTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTccag                                                            <  1:2738704/93‑1 (MQ=255)

                              gTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGGATCCTccagc                                                           >  2:3319350/1‑93 (MQ=255)

                               ttATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTccagcc                                                          >  1:423571/1‑93 (MQ=255)

                                  tGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAgg                                                       >  2:6664407/1‑93 (MQ=255)

                                   gCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGc                                                      <  1:3560023/93‑1 (MQ=255)

                                   gCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGc                                                      <  1:6476491/93‑1 (MQ=255)

                                    ctctCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCg                                                     >  2:1449573/1‑93 (MQ=255)

                                     tctcCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGt                                                    >  2:3574827/1‑93 (MQ=255)

                                     tctcCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGt                                                    >  1:1267077/1‑93 (MQ=255)

                                      ctcCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTa                                                   <  1:218592/93‑1 (MQ=255)

                                      ctcCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTa                                                   >  1:2374316/1‑93 (MQ=255)

                                      ctcCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGGTGATGGGTATCGTCCAGCCAGGCGTa                                                   >  1:5215007/1‑93 (MQ=255)

                                        ccccAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTaaa                                                 <  2:4895478/93‑1 (MQ=255)

                                         cccAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAAc                                                >  1:6164468/1‑93 (MQ=255)

                                         cccAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAAc                                                >  1:3058424/1‑93 (MQ=255)

                                          ccAGCCTGACGACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACa                                               >  2:5655009/1‑93 (MQ=255)

                                           cAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGGCAGGCGTAAACAt                                              >  1:3220737/1‑93 (MQ=255)

                                           cAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACAt                                              <  2:1852555/93‑1 (MQ=255)

                                             gCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATaa                                            >  1:1829792/1‑92 (MQ=255)

                                             gCCTGACAACAATAAACCTCTTCCGGGGAGGTTCCGCCACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATAc                                            >  2:4149951/1‑93 (MQ=255)

                                              ccTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACa                                           <  1:4623678/93‑1 (MQ=255)

                                               cTGACAACAATAAGCCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAg                                          <  1:1080542/93‑1 (MQ=255)

                                               cTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCCGGCGTAAACATACAg                                          >  2:3157216/1‑93 (MQ=255)

                                               cTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAg                                          <  1:8283422/93‑1 (MQ=255)

                                               cTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAg                                          >  1:5737158/1‑93 (MQ=255)

                                                tGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTTCT‑CCGGTTGATGGGTATTCTCCAGCCAGGTCTAAACATACAGt                                         >  2:5489435/1‑93 (MQ=255)

                                                  caatcaaTATTCCTCGTCCGGATAGTTCCCGCTACACTGGCAGGCCTGCT‑CCGGTTTATGGGTATCCTCCAGCCAGGCGTAAACATACAGttt                                       <  2:3697698/89‑1 (MQ=255)

                                                  acaacaATAAACCTCTTCCGGGGAGTTTCCGCAACACTGGCAGGGCTTCTCCCGGTTG‑TCGGCATCCTCCACACAGGCATACACATACAGttt                                       >  2:4947775/1‑93 (MQ=255)

                                                      caaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCcac                                   >  2:1308409/1‑93 (MQ=255)

                                                       aaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCcaca                                  >  1:7409353/1‑93 (MQ=255)

                                                        aTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACaa                                 <  1:5568731/93‑1 (MQ=255)

                                                            aCCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTc                             <  2:6560072/93‑1 (MQ=255)

                                                            aCCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTc                             <  1:1510830/93‑1 (MQ=255)

                                                            aCCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTc                             <  2:8230593/93‑1 (MQ=255)

                                                            aCCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTc                             <  1:5076893/93‑1 (MQ=255)

                                                             ccTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTcc                            >  1:2912997/1‑93 (MQ=255)

                                                              ctctTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACGGTTTCCACAAACTccc                           <  2:7819180/93‑1 (MQ=255)

                                                              ctctTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTccc                           <  2:8662943/93‑1 (MQ=255)

                                                               tctTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCGGGCGGAACCATACAGGTTCCACACACTACCt                          >  2:3398468/1‑93 (MQ=255)

                                                               tctTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCt                          >  1:6585972/1‑93 (MQ=255)

                                                               tctTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCt                          >  1:3738080/1‑93 (MQ=255)

                                                                 ttCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTaa                        >  1:4456603/1‑93 (MQ=255)

                                                                  tCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAAc                       <  1:7521033/93‑1 (MQ=255)

                                                                   ccGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACg                      >  2:7033864/1‑93 (MQ=255)

                                                                    cGGGGAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACgg                     <  1:2421054/93‑1 (MQ=255)

                                                                     ggggAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGc                    >  1:451324/1‑93 (MQ=255)

                                                                      gggAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCg                   >  1:5394947/1‑93 (MQ=255)

                                                                       ggggTTTCCGCGACACTGGGAGGGCTGCG‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgt                  <  1:2606234/90‑1 (MQ=255)

                                                                       ggAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgt                  <  2:5256963/93‑1 (MQ=255)

                                                                       ggAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgt                  <  1:8460085/93‑1 (MQ=255)

                                                                       ggAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgt                  <  1:6772123/93‑1 (MQ=255)

                                                                       ggAGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgt                  <  1:1950093/93‑1 (MQ=255)

                                                                         aGTTTCCGCGCCACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgtgg                <  2:6930967/93‑1 (MQ=255)

                                                                         aGTTTCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgtgg                <  1:4685416/93‑1 (MQ=255)

                                                                           tttCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgtggtg              <  1:2453995/93‑1 (MQ=255)

                                                                           tttCCGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgtggtg              <  1:1708591/93‑1 (MQ=255)

                                                                              ccGCGACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTGTCCACCAACTCCCTAACGGCgtggtcacc           >  1:695991/1‑89 (MQ=255)

                                                                               cgcgACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCg          <  2:3001218/93‑1 (MQ=255)

                                                                                gcgACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGa         >  1:327991/1‑93 (MQ=255)

                                                                                gcgACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGa         >  2:230526/1‑93 (MQ=255)

                                                                                  gACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGt       >  2:2926356/1‑93 (MQ=255)

                                                                                  gACACTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGt       >  1:773387/1‑93 (MQ=255)

                                                                                   acacTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTc      <  1:1440339/93‑1 (MQ=255)

                                                                                   acacTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTc      <  1:4979902/93‑1 (MQ=255)

                                                                                    cacTGGCAGGTCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAACCTCCCTAACGGCGTGGTGAGCGAGTCa     <  2:6136501/93‑1 (MQ=255)

                                                                                    cacTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCa     <  1:96086/93‑1 (MQ=255)

                                                                                    cacTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCa     <  1:3841209/93‑1 (MQ=255)

                                                                                      cTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAgg   >  2:5597117/1‑93 (MQ=255)

                                                                                      cTGGCAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGAGAGCGAGTCAgg   <  2:743746/93‑1 (MQ=255)

CCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCT‑CCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGT  >  NC_000913/2327489‑2327668

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