BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	27,903	T→C	L204P (CTG→CCG)	yaaF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	27,903	0	T	C	99.4%	1070.6 / ‑3.4	340	L204P (CTG→CCG)	yaaF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (179/159); total (180/160)
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.03e-01

Only 100 of 354 total aligned reads displayed.

CCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACG  >  NC_000913/27813‑27992

ccTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGc                                                                                         >  2:4951731/1‑93 (MQ=255)

ccTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGc                                                                                         <  1:2878996/93‑1 (MQ=255)

 cTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGcc                                                                                        >  1:6202836/1‑93 (MQ=255)

    tgtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCa                                                                                     <  2:4051859/93‑1 (MQ=255)

     gtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAg                                                                                    <  1:3511414/93‑1 (MQ=255)

      tCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGt                                                                                   >  1:4206572/1‑93 (MQ=255)

       cTTCCGCAGTGGTATTGAAATCGTCATGTTCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGtt                                                                                  <  1:7390289/93‑1 (MQ=255)

        ttCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTa                                                                                 >  2:1382418/1‑93 (MQ=255)

          cctcAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTaaa                                                                               <  1:1653745/90‑1 (MQ=255)

          ccGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTATCTCCTGACTATCGCTCTACACCGCCGCAGTTaaa                                                                               <  2:7329916/93‑1 (MQ=255)

           cGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAAc                                                                              <  2:1111909/93‑1 (MQ=255)

           cGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAAc                                                                              <  1:7639694/93‑1 (MQ=255)

            gCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAAcc                                                                             >  2:6823475/1‑93 (MQ=255)

            gCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAAcc                                                                             <  1:4256464/93‑1 (MQ=255)

              aGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGt                                                                           >  1:5549613/1‑93 (MQ=255)

               gTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTa                                                                          <  1:2790681/93‑1 (MQ=255)

                  gTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCg                                                                       >  1:1560538/1‑93 (MQ=255)

                  gTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCg                                                                       >  1:4525922/1‑93 (MQ=255)

                   tATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCgg                                                                      >  1:4434975/1‑93 (MQ=255)

                   tATTGAAATCGTCATGTGCGGTGTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCgg                                                                      >  2:7001126/1‑93 (MQ=255)

                     ttGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGa                                                                    >  1:2673693/1‑93 (MQ=255)

                      tGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGaa                                                                   <  1:6416667/93‑1 (MQ=255)

                        aaaTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGaaaa                                                                 <  1:5390995/93‑1 (MQ=255)

                        aaaTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGaaaa                                                                 >  2:4311332/1‑93 (MQ=255)

                        aaaTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGaaaa                                                                 <  1:2239162/93‑1 (MQ=255)

                         aaTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAAt                                                                <  2:1644699/93‑1 (MQ=255)

                           tCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGc                                                              >  2:4690963/1‑93 (MQ=255)

                           tCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGc                                                              <  2:476720/93‑1 (MQ=255)

                            cGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCt                                                             <  1:172401/93‑1 (MQ=255)

                             gTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCtt                                                            <  1:975240/93‑1 (MQ=255)

                               cATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCa                                                          >  2:2809727/1‑93 (MQ=255)

                               cATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCa                                                          <  1:7152699/93‑1 (MQ=255)

                               cATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCa                                                          <  1:7061381/93‑1 (MQ=255)

                                aTGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCAc                                                         <  1:771079/93‑1 (MQ=255)

                                  gtgCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGc                                                       <  1:1228090/93‑1 (MQ=255)

                                   tgCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGcc                                                      >  1:7241372/1‑93 (MQ=255)

                                   tgCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGcc                                                      <  1:1908233/93‑1 (MQ=255)

                                    gCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGccc                                                     <  2:7108838/93‑1 (MQ=255)

                                     cGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCt                                                    <  1:7825057/93‑1 (MQ=255)

                                       gTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGt                                                  >  2:3805289/1‑93 (MQ=255)

                                        tttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGtt                                                 >  2:6047887/1‑93 (MQ=255)

                                         ttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGttt                                                <  2:3166973/93‑1 (MQ=255)

                                          tGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTa                                               >  1:4140006/1‑93 (MQ=255)

                                          tGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTa                                               >  2:4167734/1‑93 (MQ=255)

                                           ggATGTCACCAATCATGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAg                                              <  1:341307/93‑1 (MQ=255)

                                           ggATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAg                                              <  2:3339495/93‑1 (MQ=255)

                                             aTGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGcc                                            >  1:2398265/1‑93 (MQ=255)

                                             aTGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGcc                                            >  2:6685387/1‑93 (MQ=255)

                                              ttgcACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCa                                           <  2:7418335/90‑1 (MQ=255)

                                               gTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACCCCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCAc                                          <  1:1391305/93‑1 (MQ=255)

                                                tCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACt                                         >  2:257918/1‑93 (MQ=255)

                                                 cACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTa                                        >  1:1205863/1‑93 (MQ=255)

                                                  aCCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAc                                       >  1:3896209/1‑93 (MQ=255)

                                                  aCCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAc                                       <  2:5054563/93‑1 (MQ=255)

                                                  aCCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAc                                       >  2:6281217/1‑93 (MQ=255)

                                                   ccAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAcc                                      >  2:1552307/1‑93 (MQ=255)

                                                   ccAATCAGGCAATATTAACTCCTCACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAcc                                      <  2:5659053/93‑1 (MQ=255)

                                                         aGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCg                                >  2:864924/1‑93 (MQ=255)

                                                          ggCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCgg                               >  1:6156177/1‑93 (MQ=255)

                                                              atatTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGt                           >  2:235311/1‑93 (MQ=255)

                                                              atatTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGt                           >  2:2895430/1‑93 (MQ=255)

                                                              atatTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGt                           >  2:7665905/1‑93 (MQ=255)

                                                                atTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTAt                         <  2:639988/93‑1 (MQ=255)

                                                                 ttAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATg                        >  1:3172506/1‑93 (MQ=255)

                                                                   aaCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCa                      >  1:6860614/1‑93 (MQ=255)

                                                                     cTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaaa                    <  1:7013615/93‑1 (MQ=255)

                                                                      tCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAg                   >  2:1876781/1‑93 (MQ=255)

                                                                      tCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAg                   <  2:6446107/93‑1 (MQ=255)

                                                                       ccTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGc                  <  2:6218375/93‑1 (MQ=255)

                                                                       ccTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGc                  <  2:5521244/93‑1 (MQ=255)

                                                                       ccTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGc                  <  1:2550371/93‑1 (MQ=255)

                                                                       ccTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGc                  <  1:909671/93‑1 (MQ=255)

                                                                        cTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCg                 <  2:7543299/93‑1 (MQ=255)

                                                                        cTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCg                 >  2:2173854/1‑93 (MQ=255)

                                                                         tGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCgg                >  2:6860163/1‑93 (MQ=255)

                                                                         tGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCgg                >  1:4895453/1‑93 (MQ=255)

                                                                         tGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCgg                <  2:5592042/93‑1 (MQ=255)

                                                                          gACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGc               <  2:325484/93‑1 (MQ=255)

                                                                          gACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGc               >  1:5637413/1‑93 (MQ=255)

                                                                            cTATCTCTCTACCCCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt             <  2:523745/93‑1 (MQ=255)

                                                                            cTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt             <  1:5867092/93‑1 (MQ=255)

                                                                            cTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt             <  1:4366050/93‑1 (MQ=255)

                                                                             tATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCttt            >  1:2059400/1‑92 (MQ=255)

                                                                             tATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTg            >  1:3291496/1‑93 (MQ=255)

                                                                             tATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTg            >  1:1032238/1‑93 (MQ=255)

                                                                             tATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTg            >  2:2038283/1‑93 (MQ=255)

                                                                             tATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTg            <  2:5991853/93‑1 (MQ=255)

                                                                              aTCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtttg           >  1:6530354/1‑91 (MQ=255)

                                                                              aTCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGc           <  1:3979256/93‑1 (MQ=255)

                                                                               tctctcTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCg          >  1:6631681/1‑93 (MQ=255)

                                                                               tctctcTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCg          >  1:1061789/1‑93 (MQ=255)

                                                                                 tctctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGTGaa        >  2:4506784/1‑93 (MQ=255)

                                                                                 tctctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGaa        <  2:3608577/93‑1 (MQ=255)

                                                                                 tctctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGaa        >  1:7213525/1‑93 (MQ=255)

                                                                                 cctctACACCGCCGCAGTTAAACCGTCCCGGTAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGaa        <  1:4021440/92‑1 (MQ=255)

                                                                                    ctACCCCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGc     <  1:6576441/93‑1 (MQ=255)

                                                                                    ctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGc     <  1:4827280/93‑1 (MQ=255)

                                                                                     tACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTTCGGATGCa    >  1:3032118/1‑93 (MQ=255)

                                                                                      acacCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGGATGCAc   >  2:2955/1‑93 (MQ=255)

                                                                                       cacCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACg  <  2:3662071/93‑1 (MQ=255)

CCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACG  >  NC_000913/27813‑27992

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