BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,369,637	T→C	R345R (CGT→CGC)	mopA →	GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,369,637	0	T	C	99.6%	823.6 / NA	256	R345R (CGT→CGC)	mopA	GroEL
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (121/134); minor base A (0/1); total (121/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

AGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGG  >  NC_000913/4369547‑4369726

 gACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGtt                                                                                        >  2:3625287/1‑93 (MQ=255)

 gACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGtt                                                                                        >  1:5211203/1‑93 (MQ=255)

   ccTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGc                                                                                      >  2:5249354/1‑93 (MQ=255)

   ccTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGc                                                                                      <  1:3407781/93‑1 (MQ=255)

    cTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCt                                                                                     >  1:6659434/1‑93 (MQ=255)

    cTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCt                                                                                     >  2:4305519/1‑93 (MQ=255)

      gggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGGCGCGTTGCTCa                                                                                   >  1:2336077/1‑93 (MQ=255)

      gggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCa                                                                                   >  2:878543/1‑93 (MQ=255)

      gggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCa                                                                                   >  1:3366403/1‑93 (MQ=255)

      gggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGCTGCTCa                                                                                   >  1:6285501/1‑93 (MQ=255)

       ggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAg                                                                                  >  2:946831/1‑93 (MQ=255)

       ggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAg                                                                                  <  1:4220350/93‑1 (MQ=255)

        gTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGa                                                                                 >  1:2863024/1‑93 (MQ=255)

        gTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGa                                                                                 <  2:6345024/93‑1 (MQ=255)

        gTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGa                                                                                 >  2:4872404/1‑93 (MQ=255)

        gTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGa                                                                                 <  1:5634079/93‑1 (MQ=255)

        gTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGa                                                                                 <  2:4026771/93‑1 (MQ=255)

        gTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGa                                                                                 >  2:5322071/1‑93 (MQ=255)

         tCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGAt                                                                                <  2:6089504/93‑1 (MQ=255)

         tCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGAt                                                                                <  2:3446644/93‑1 (MQ=255)

         tCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGAt                                                                                <  2:1915520/93‑1 (MQ=255)

          cAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATc                                                                               <  2:6619024/93‑1 (MQ=255)

          cAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATc                                                                               <  2:2066208/93‑1 (MQ=255)

          cAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATc                                                                               <  1:2979364/93‑1 (MQ=255)

           aGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATcc                                                                              >  2:6005168/1‑93 (MQ=255)

            ggCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCg                                                                             >  2:1103357/1‑93 (MQ=255)

               tAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTca                                                                          >  1:2627084/1‑93 (MQ=255)

               tAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTca                                                                          <  1:3129476/93‑1 (MQ=255)

                    gtgtTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGa                                                                     <  2:6233076/93‑1 (MQ=255)

                    gtgtTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGa                                                                     <  1:6176501/93‑1 (MQ=255)

                     tgttgtGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGAt                                                                    <  1:2122890/93‑1 (MQ=255)

                         gtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaa                                                                <  2:2789101/93‑1 (MQ=255)

                         gtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaa                                                                <  2:4998026/93‑1 (MQ=255)

                          tgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaag                                                               <  2:699587/93‑1 (MQ=255)

                           gATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaaga                                                              <  1:4663534/93‑1 (MQ=255)

                             tCAACAAAGACACCACCACTATCATCGATGGCGTAGGTGAAGAAGCTGCAATCCAGGCCCGCGTTGCTCAGATCCGTCAGCAGATTGaagaag                                                            <  1:3963827/93‑1 (MQ=255)

                              caacaaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGc                                                           <  2:5875429/93‑1 (MQ=255)

                               aacaaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCa                                                          >  2:2515737/1‑93 (MQ=255)

                               aacaaAGACACCACCACTATCATCGATGGCGTGAGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCa                                                          >  2:6547049/1‑93 (MQ=255)

                                 caaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAAc                                                        <  2:4946610/93‑1 (MQ=255)

                                 caaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAAc                                                        <  2:2605393/93‑1 (MQ=255)

                                 caaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAAc                                                        <  1:3601482/93‑1 (MQ=255)

                                   aaGACACCACCACTATCATCGATGGCGTGGGTGACGAAGCTGCCATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACtt                                                      <  2:1180007/93‑1 (MQ=255)

                                    aGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTc                                                     >  1:1505888/1‑93 (MQ=255)

                                    aGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTc                                                     >  2:3872709/1‑93 (MQ=255)

                                    aGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTc                                                     >  1:1568864/1‑93 (MQ=255)

                                    aGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGCCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTc                                                     <  2:1631969/93‑1 (MQ=255)

                                     gccaccaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCt                                                    <  2:1786431/91‑1 (MQ=255)

                                     gACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCt                                                    <  2:2894502/93‑1 (MQ=255)

                                      acacCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTg                                                   <  1:3181410/93‑1 (MQ=255)

                                      acacCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTg                                                   <  2:3140002/93‑1 (MQ=255)

                                       caccaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGa                                                  <  2:3225486/93‑1 (MQ=255)

                                       caccaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGa                                                  >  1:2726792/1‑93 (MQ=255)

                                        tccaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGAGCCGTCAGCAGATTGAAGAAGCAACTTCTGAc                                                 <  2:4108983/92‑1 (MQ=255)

                                        cccaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGAc                                                 <  1:4551298/92‑1 (MQ=255)

                                        cccaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGAc                                                 <  2:3409712/92‑1 (MQ=255)

                                        accaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGAc                                                 >  1:3296778/1‑93 (MQ=255)

                                          caccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTa                                               >  1:220142/1‑93 (MQ=255)

                                           accacTATCATCGATGGCGTGGGTGAAGAATCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTAc                                              >  2:2664575/1‑93 (MQ=255)

                                           accacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTag                                              >  1:1104792/1‑92 (MQ=255)

                                             cacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGa                                            >  2:2266713/1‑93 (MQ=255)

                                             cacTATCATCGATGGCGTGGATGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGa                                            <  2:1499489/93‑1 (MQ=255)

                                              acTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGAc                                           >  2:2098429/1‑93 (MQ=255)

                                                  tcatcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTg                                       >  2:825660/1‑93 (MQ=255)

                                                  tcatcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTg                                       <  2:2026315/93‑1 (MQ=255)

                                                  tcatcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTg                                       >  2:6364874/1‑93 (MQ=255)

                                                   catcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGa                                      <  2:4387103/93‑1 (MQ=255)

                                                   catcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGa                                      <  1:830235/93‑1 (MQ=255)

                                                     tcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaa                                    <  1:440536/93‑1 (MQ=255)

                                                       gATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaaaa                                  <  1:904503/93‑1 (MQ=255)

                                                       gATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaaaa                                  <  1:1027503/93‑1 (MQ=255)

                                                        aTGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAAc                                 >  1:5348818/1‑93 (MQ=255)

                                                         tGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACt                                <  2:1120151/93‑1 (MQ=255)

                                                             gTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAg                            <  1:4792665/93‑1 (MQ=255)

                                                               gggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGa                          <  1:2886401/93‑1 (MQ=255)

                                                                  tGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAcg                       >  2:6239469/1‑93 (MQ=255)

                                                                  tGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAcg                       <  1:6013415/93‑1 (MQ=255)

                                                                   gaagaaGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAcgc                      <  1:208668/93‑1 (MQ=255)

                                                                   gaagaaGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAcgc                      <  2:4722051/93‑1 (MQ=255)

                                                                    aagaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAcgcg                     >  1:4638901/1‑93 (MQ=255)

                                                                    aagaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAcgcg                     >  2:5478401/1‑93 (MQ=255)

                                                                    aagaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAcgcg                     <  1:1344290/93‑1 (MQ=255)

                                                                      gaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTa                   >  2:145257/1‑93 (MQ=255)

                                                                      gaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTa                   <  1:5111954/93‑1 (MQ=255)

                                                                      gaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTa                   <  1:6305474/93‑1 (MQ=255)

                                                                      gaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTa                   <  1:4037663/93‑1 (MQ=255)

                                                                       aagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAg                  <  2:3738914/93‑1 (MQ=255)

                                                                         gCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCg                >  2:5637323/1‑93 (MQ=255)

                                                                          cTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa               >  2:3846158/1‑93 (MQ=255)

                                                                          cTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa               <  1:983800/93‑1 (MQ=255)

                                                                            gCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGaaa             <  1:3683604/93‑1 (MQ=255)

                                                                            gCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGaaa             >  1:3312575/1‑93 (MQ=255)

                                                                              aaTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACt           <  1:4298807/93‑1 (MQ=255)

                                                                               aTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTg          <  2:2460499/93‑1 (MQ=255)

                                                                               aTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTg          >  1:2498881/1‑93 (MQ=255)

                                                                               aTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTg          >  1:5495600/1‑93 (MQ=255)

                                                                                tCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTgg         >  2:3189151/1‑93 (MQ=255)

                                                                                  cAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCa       <  2:2557660/93‑1 (MQ=255)

                                                                                    gggCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAgg     >  1:2797862/1‑93 (MQ=255)

                                                                                     ggCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAggc    <  1:528425/93‑1 (MQ=255)

AGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGG  >  NC_000913/4369547‑4369726

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