BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,515,312	T→C	M181V (ATG→GTG)	b2395 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,515,312	0	T	C	99.7%	1029.7 / ‑4.7	315	M181V (ATG→GTG)	b2395	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (144/170); minor base G (0/1); total (144/171)
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.41e-01

Only 100 of 337 total aligned reads displayed.

CTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTA  >  NC_000913/2515222‑2515401

ctTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAcgt                                                                                         <  2:4443126/93‑1 (MQ=255)

ctTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAcgt                                                                                         <  2:49926/93‑1 (MQ=255)

 tttGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAcgtt                                                                                        >  1:9027621/1‑93 (MQ=255)

  ttGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAcgttg                                                                                       >  1:2552301/1‑93 (MQ=255)

  ttGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAcgttg                                                                                       <  2:6076143/93‑1 (MQ=255)

   tGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAcgttgt                                                                                      <  1:4314477/93‑1 (MQ=255)

    gCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTa                                                                                     <  2:9814751/93‑1 (MQ=255)

     cccAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTaa                                                                                    >  2:9058755/1‑93 (MQ=255)

       cAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAAt                                                                                  >  1:1732079/1‑93 (MQ=255)

        aaCGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATc                                                                                 <  2:1868164/93‑1 (MQ=255)

        aaCGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATc                                                                                 <  2:6436979/93‑1 (MQ=255)

        aaCGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATc                                                                                 <  2:8157629/93‑1 (MQ=255)

          cgacgGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAg                                                                               <  2:2287438/93‑1 (MQ=255)

           gacgGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGc                                                                              <  2:6257689/93‑1 (MQ=255)

            acgGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGca                                                                             <  2:8808810/93‑1 (MQ=255)

              ggAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGcaca                                                                           <  2:8281397/93‑1 (MQ=255)

              ggAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGcaca                                                                           <  2:7085784/93‑1 (MQ=255)

                 gCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCg                                                                        >  2:3213804/1‑93 (MQ=255)

                   gATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGt                                                                      <  2:5558515/93‑1 (MQ=255)

                   gATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGt                                                                      >  1:6123587/1‑93 (MQ=255)

                    atatCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTg                                                                     >  2:3624271/1‑93 (MQ=255)

                    atatCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTg                                                                     >  1:9406639/1‑93 (MQ=255)

                    atatCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTg                                                                     >  2:1123412/1‑93 (MQ=255)

                    atatCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTg                                                                     >  1:3464885/1‑93 (MQ=255)

                     tatCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGa                                                                    <  2:1666530/93‑1 (MQ=255)

                     tatCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGa                                                                    >  2:8904975/1‑93 (MQ=255)

                      atCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGaa                                                                   >  1:7885064/1‑93 (MQ=255)

                        cccTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGaaaa                                                                 <  2:8553085/93‑1 (MQ=255)

                        cccTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGaaaa                                                                 <  2:5114474/93‑1 (MQ=255)

                        cccTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGaaaa                                                                 >  2:4635050/1‑93 (MQ=255)

                         ccTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAg                                                                <  1:1759617/93‑1 (MQ=255)

                         ccTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAg                                                                <  2:939062/93‑1 (MQ=255)

                           tGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTGAATCAGCACAGCGGTGAAAAGGc                                                              <  2:2659439/93‑1 (MQ=255)

                           tGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGCGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGc                                                              >  1:7267918/1‑93 (MQ=255)

                            gcgcCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCt                                                             <  1:3052265/93‑1 (MQ=255)

                             cgcCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCtt                                                            <  2:2793340/93‑1 (MQ=255)

                             cgcCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCtt                                                            <  2:4872608/93‑1 (MQ=255)

                              gcCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTa                                                           <  2:3685647/93‑1 (MQ=255)

                              gcCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTa                                                           <  2:4145812/93‑1 (MQ=255)

                              gcCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTa                                                           >  2:481446/1‑93 (MQ=255)

                              gcCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTa                                                           <  2:5436096/93‑1 (MQ=255)

                                cacaATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGc                                                         >  1:8990696/1‑93 (MQ=255)

                                cacaATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGc                                                         <  1:4500757/93‑1 (MQ=255)

                                   aaTATCTGCGCAAAGTGGGGCCTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTACATCAGCACACAGGTGAAAAGGCTTAGCaaa                                                      >  1:4190160/1‑93 (MQ=255)

                                   aaTATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCaaa                                                      >  1:4959400/1‑93 (MQ=255)

                                    atatCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAAt                                                     <  1:9515589/93‑1 (MQ=255)

                                    atatCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAAt                                                     <  2:75633/93‑1 (MQ=255)

                                     tatCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATc                                                    >  1:2701952/1‑93 (MQ=255)

                                       tCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCaa                                                  <  2:3727922/93‑1 (MQ=255)

                                       tCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCaa                                                  >  2:1494148/1‑93 (MQ=255)

                                        cTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAAc                                                 <  1:3715130/93‑1 (MQ=255)

                                        cTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGACAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAAc                                                 <  2:9730179/93‑1 (MQ=255)

                                          gcgcAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGa                                               <  1:6866996/93‑1 (MQ=255)

                                          gcgaAAAGTGGGGCCTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGa                                               <  1:1365900/93‑1 (MQ=255)

                                            gcAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGAcc                                             >  1:8333667/1‑93 (MQ=255)

                                             cAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCg                                            >  1:572781/1‑93 (MQ=255)

                                              aaaGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGt                                           >  2:1374179/1‑93 (MQ=255)

                                               aaGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTg                                          >  1:8137733/1‑93 (MQ=255)

                                                aGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGa                                         <  1:2232176/93‑1 (MQ=255)

                                                 gTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGaa                                        <  2:10020508/93‑1 (MQ=255)

                                                  gggggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGaaa                                       <  2:7281843/92‑1 (MQ=255)

                                                   ggggACTTACAATCATGCGAGTGAGATAGTAGGAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGaaaa                                      <  1:6670758/93‑1 (MQ=255)

                                                   ggggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCGCAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGaaaa                                      <  2:5781623/93‑1 (MQ=255)

                                                   ggggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGaaaa                                      >  2:9541256/1‑93 (MQ=255)

                                                   ggggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGaaaa                                      <  2:3668534/93‑1 (MQ=255)

                                                    gggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAAt                                     >  2:7893247/1‑93 (MQ=255)

                                                     ggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATg                                    >  2:6859735/1‑93 (MQ=255)

                                                      gACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATgg                                   >  2:4286907/1‑93 (MQ=255)

                                                       aCTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATgga                                  >  1:7322003/1‑92 (MQ=255)

                                                        cTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCa                                 <  1:6647512/93‑1 (MQ=255)

                                                          tACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATc                               <  2:3898637/93‑1 (MQ=255)

                                                          tACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATc                               <  2:2057304/93‑1 (MQ=255)

                                                              aTCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCa                           >  1:518198/1‑93 (MQ=255)

                                                              aTCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCa                           <  2:5590263/93‑1 (MQ=255)

                                                               tCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCAt                          >  1:354009/1‑93 (MQ=255)

                                                               tCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCAt                          <  2:8997295/93‑1 (MQ=255)

                                                               tCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCAt                          >  2:1581443/1‑93 (MQ=255)

                                                               tCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCAt                          >  1:9711399/1‑93 (MQ=255)

                                                                cATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATc                         >  1:5941238/1‑93 (MQ=255)

                                                                  tGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGc                       <  1:1657844/93‑1 (MQ=255)

                                                                       gTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaa                  >  1:5464430/1‑93 (MQ=255)

                                                                       gTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaa                  >  1:9779174/1‑93 (MQ=255)

                                                                        tGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaaa                 >  1:7099557/1‑93 (MQ=255)

                                                                         gagaTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaaaa                <  2:5970481/93‑1 (MQ=255)

                                                                         gagaTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaaaa                >  1:2788354/1‑93 (MQ=255)

                                                                           gaTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAgg              <  1:6358061/93‑1 (MQ=255)

                                                                            aTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGt             <  1:799057/93‑1 (MQ=255)

                                                                            aTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGt             <  1:8676691/93‑1 (MQ=255)

                                                                            aTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGt             <  1:532344/93‑1 (MQ=255)

                                                                              agtagAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTgg           <  2:3485661/93‑1 (MQ=255)

                                                                                tagAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGat         <  2:9449864/93‑1 (MQ=255)

                                                                                 agAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGata        >  2:1044513/1‑93 (MQ=255)

                                                                                  gAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGatat       >  1:660930/1‑93 (MQ=255)

                                                                                   aaaaGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATc      >  2:8698197/1‑93 (MQ=255)

                                                                                   aaaaGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATc      <  1:9092603/93‑1 (MQ=255)

                                                                                   aaaaGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATc      >  2:7369621/1‑93 (MQ=255)

                                                                                   aaaaGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATc      >  1:1454605/1‑93 (MQ=255)

                                                                                    aaaGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCt     >  2:7946522/1‑93 (MQ=255)

                                                                                    aaaGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCt     <  1:7372173/93‑1 (MQ=255)

                                                                                    aaaGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCt     >  1:1052481/1‑93 (MQ=255)

CTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTA  >  NC_000913/2515222‑2515401

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