BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,919,873	G→A	D24N (GAT→AAT)	b1836 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,919,873	0	G	A	98.3%	616.6 / ‑2.2	179	D24N (GAT→AAT)	b1836	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (2/1); new base A (86/90); total (88/91)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

Only 100 of 180 total aligned reads displayed.

AGCTAAATGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTGATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCAGTATCCCTGCCCTG  >  NC_000913/1919783‑1919962

aGCTAAATGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAAt                                                                                         >  2:2132006/1‑93 (MQ=255)

aGCTAAATGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAAt                                                                                         <  1:4295685/93‑1 (MQ=255)

aGCTAAATGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAAt                                                                                         <  2:3063498/93‑1 (MQ=255)

 gCTAAATGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTGATc                                                                                        <  1:3319973/93‑1 (MQ=255)

 gCTAAATGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATc                                                                                        >  1:2091258/1‑93 (MQ=255)

  cTAAATGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCg                                                                                       <  1:4391128/93‑1 (MQ=255)

   tAAATGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGa                                                                                      >  2:3820058/1‑93 (MQ=255)

   tAAATGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGCGAATCGCCAGGTAATCGa                                                                                      >  2:2934373/1‑93 (MQ=255)

    aaaTGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGaa                                                                                     <  1:3087097/93‑1 (MQ=255)

    aaaTGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCCaa                                                                                     >  2:1888497/1‑93 (MQ=255)

     aaTGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAAc                                                                                    >  2:4483082/1‑93 (MQ=255)

      aTGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACg                                                                                   <  1:837207/93‑1 (MQ=255)

      aTGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACg                                                                                   >  1:2933003/1‑93 (MQ=255)

       tGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAActt                                                                                  <  2:2056643/93‑3 (MQ=255)

       tGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGt                                                                                  <  1:2266649/93‑1 (MQ=255)

        ggAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGtt                                                                                 <  2:1216984/93‑1 (MQ=255)

        ggAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGtt                                                                                 <  1:3269659/93‑1 (MQ=255)

          ggaacACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTaa                                                                               <  1:1044296/89‑1 (MQ=255)

           gagCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAAc                                                                              >  1:2604249/1‑93 (MQ=255)

             gCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCa                                                                            >  1:3980575/1‑93 (MQ=255)

             gCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCa                                                                            >  2:2243256/1‑93 (MQ=255)

              gaCGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCAt                                                                           <  1:4250354/92‑1 (MQ=255)

               aCGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCAtt                                                                          >  2:4604544/1‑93 (MQ=255)

                cGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTc                                                                         <  2:2517584/93‑1 (MQ=255)

                 gAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTcc                                                                        <  1:371418/93‑1 (MQ=255)

                 gAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTcc                                                                        <  1:864321/93‑1 (MQ=255)

                  aaGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCt                                                                       >  2:1178708/1‑93 (MQ=255)

                   aGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCtt                                                                      >  1:2873990/1‑93 (MQ=255)

                   aGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCtt                                                                      <  2:3589185/93‑1 (MQ=255)

                    gATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTg                                                                     >  2:4078273/1‑93 (MQ=255)

                    gATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTg                                                                     >  2:4078272/1‑93 (MQ=255)

                     aTGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCtttt                                                                    >  2:864546/1‑91 (MQ=255)

                     aTGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGt                                                                    >  2:4663502/1‑93 (MQ=255)

                      tGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTa                                                                   >  1:4631875/1‑93 (MQ=255)

                       gAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTaa                                                                  <  1:1029160/93‑1 (MQ=255)

                        aaaaCCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTaaa                                                                 <  1:2068657/93‑1 (MQ=255)

                        aaaaCCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTaaa                                                                 >  2:3150696/1‑93 (MQ=255)

                        aaaaCCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTaaa                                                                 <  2:3605456/93‑1 (MQ=255)

                         aaaCCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTTTAAAt                                                                >  2:3494001/1‑93 (MQ=255)

                         aaaCCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAAt                                                                >  2:4202688/1‑93 (MQ=255)

                         aaaCCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAAt                                                                >  1:4670206/1‑93 (MQ=255)

                          aaCCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATc                                                               <  1:4049994/93‑1 (MQ=255)

                            ccAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTg                                                             >  1:2236266/1‑93 (MQ=255)

                            ccAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTg                                                             >  2:1721627/1‑93 (MQ=255)

                              aGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGAc                                                           >  2:4428457/1‑93 (MQ=255)

                               gtgtGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGAcc                                                          <  1:999587/93‑1 (MQ=255)

                                 gtgCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGAcccc                                                        <  2:897708/93‑1 (MQ=255)

                                 gtgCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGAcccc                                                        >  2:1624610/1‑93 (MQ=255)

                                      cATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGAttt                                                   <  2:3917891/93‑1 (MQ=255)

                                       aTAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTta                                                  >  1:4537304/1‑93 (MQ=255)

                                       aTAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTta                                                  >  1:1912019/1‑93 (MQ=255)

                                        tAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTtat                                                 <  2:2824237/93‑1 (MQ=255)

                                        tAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTtat                                                 >  2:3394411/1‑93 (MQ=255)

                                         aGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTtatg                                                <  1:4604544/93‑1 (MQ=255)

                                          ggCGCTTTTGAAATCGGCGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTtatgt                                               <  1:4586164/93‑1 (MQ=255)

                                            cgcTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTtatgtat                                             <  2:2091258/93‑1 (MQ=255)

                                              cTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGc                                           >  1:4557240/1‑93 (MQ=255)

                                              cTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGc                                           >  2:4311358/1‑93 (MQ=255)

                                              cTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGc                                           >  2:1949441/1‑93 (MQ=255)

                                              cTTTTGAAATCGACGACGGCGAATTACACGGTCAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGc                                           <  2:542995/93‑1 (MQ=255)

                                               ttttGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCa                                          <  1:3152049/93‑1 (MQ=255)

                                                tttGAAATCGACGACGGCGAATTACACGGTGAATCGCCCGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCac                                         >  1:2005289/1‑92 (MQ=255)

                                                tttGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCaa                                         <  1:650762/93‑1 (MQ=255)

                                                tttGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCaa                                         >  2:1702619/1‑93 (MQ=255)

                                                   gAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTg                                      <  2:2283662/93‑1 (MQ=255)

                                                    aaaTCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTgg                                     <  1:3631447/93‑1 (MQ=255)

                                                     aaTCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGa                                    >  2:2843623/1‑93 (MQ=255)

                                                      aTCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGAt                                   <  1:4410345/93‑1 (MQ=255)

                                                       tCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATg                                  <  1:1203175/93‑1 (MQ=255)

                                                       tCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATg                                  >  2:3216142/1‑93 (MQ=255)

                                                       tCGACGACGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCACGATTTATGTATGCAACTGGATg                                  >  2:3216143/1‑93 (MQ=255)

                                                        cgacgaCGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGc                                 >  2:3402114/1‑93 (MQ=255)

                                                        cgacgaCGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTGATGTATGCAACTGGATGc                                 <  2:4558734/93‑1 (MQ=255)

                                                         gacgacGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGcc                                <  1:146178/93‑1 (MQ=255)

                                                         gacgacGGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGcc                                <  2:1419786/93‑1 (MQ=255)

                                                          acgacgGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGTCCCCGATTTATGTATGCAACTGGATGCCt                               >  2:2708372/1‑93 (MQ=255)

                                                           cgacgGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTg                              <  2:3947744/93‑1 (MQ=255)

                                                           cgacgGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTg                              <  1:1107622/93‑1 (MQ=255)

                                                           cgacgGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTg                              <  2:2157642/93‑1 (MQ=255)

                                                            gacgGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTgg                             >  2:1033811/1‑93 (MQ=255)

                                                            gacgGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTgg                             <  1:2176980/93‑1 (MQ=255)

                                                              cgGCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGa                           <  1:4497525/93‑1 (MQ=255)

                                                               ggCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGAt                          >  1:4553914/1‑93 (MQ=255)

                                                               ggCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGAt                          <  1:1228756/93‑1 (MQ=255)

                                                               ggCGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGAt                          <  2:2063230/93‑1 (MQ=255)

                                                                 cGAATTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGc                        >  1:3000930/1‑93 (MQ=255)

                                                                    aTTACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGa                     >  2:688461/1‑93 (MQ=255)

                                                                     ttACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGaa                    >  1:4686873/1‑93 (MQ=255)

                                                                      tACACGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGaaa                   >  1:67311/1‑93 (MQ=255)

                                                                         acGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCa                >  1:3579055/1‑93 (MQ=255)

                                                                         acGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCa                >  2:713256/1‑93 (MQ=255)

                                                                         acGGTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGCATGCTGAAACCa                <  1:903886/93‑1 (MQ=255)

                                                                           ggTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCAGt              <  1:2809299/93‑1 (MQ=255)

                                                                            gTGAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCAGTa             <  1:2055363/93‑1 (MQ=255)

                                                                              gAATCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCAGTATc           >  1:1872603/1‑93 (MQ=255)

                                                                                aTCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCAGTATccc         <  2:288790/93‑1 (MQ=255)

                                                                                aTCGCCAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCAGTATccc         >  1:1434919/1‑93 (MQ=255)

                                                                                    ccAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCAGTATccctgcc     <  2:804450/93‑1 (MQ=255)

                                                                                    ccAGGTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCAGTATccctgcc     >  1:2378857/1‑93 (MQ=255)

                                                                                       ggTAATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCAGTATccctgccctg  >  1:2645569/1‑93 (MQ=255)

AGCTAAATGGAGAGCACGAAGATGAAAACCAGTGTGCGCATAGGCGCTTTTGAAATCGACGACGGCGAATTACACGGTGAATCGCCAGGTGATCGAACGTTAACCATTCCTTGTAAATCTGACCCCGATTTATGTATGCAACTGGATGCCTGGGATGCTGAAACCAGTATCCCTGCCCTG  >  NC_000913/1919783‑1919962

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