BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,786,415	A→G	T230A (ACG→GCG)	b3614 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,786,415	0	A	G	99.6%	823.1 / ‑5.3	258	T230A (ACG→GCG)	b3614	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (148/109); total (149/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

Only 100 of 271 total aligned reads displayed.

CGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAT  >  NC_000913/3786325‑3786504

 gCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTgcgc                                                                                        <  1:305151/93‑1 (MQ=255)

  cAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCt                                                                                       >  1:534511/1‑93 (MQ=255)

   aaCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTg                                                                                      >  2:1157956/1‑93 (MQ=255)

   aaCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTg                                                                                      <  2:2794207/93‑1 (MQ=255)

    aCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTgg                                                                                     >  1:4384868/1‑93 (MQ=255)

     cGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGt                                                                                    >  2:703005/1‑93 (MQ=255)

      ggTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTg                                                                                   <  2:2874024/93‑1 (MQ=255)

        ttCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGaa                                                                                 <  1:3863971/93‑1 (MQ=255)

        ttCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGaa                                                                                 >  1:3888744/1‑93 (MQ=255)

         tCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGaaa                                                                                <  2:2974306/93‑1 (MQ=255)

          cGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAg                                                                               >  2:3729959/1‑93 (MQ=255)

          cGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAg                                                                               <  2:2821599/93‑1 (MQ=255)

           gACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGc                                                                              >  2:2688663/1‑93 (MQ=255)

           gACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGc                                                                              >  2:2601192/1‑93 (MQ=255)

            aCCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGcc                                                                             >  1:2278763/1‑93 (MQ=255)

              cattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcag                                                                           <  2:3550776/93‑1 (MQ=255)

               attgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagc                                                                          <  1:4583454/93‑1 (MQ=255)

               attgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagc                                                                          >  2:1097563/1‑93 (MQ=255)

                ttgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagca                                                                         <  2:1261738/93‑1 (MQ=255)

                   ccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCt                                                                      <  2:3500916/93‑1 (MQ=255)

                    cattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCtt                                                                     >  1:1023143/1‑93 (MQ=255)

                     attgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTg                                                                    >  1:736776/1‑93 (MQ=255)

                      ttgGGCATCCTCCCCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGc                                                                   <  1:1259529/93‑1 (MQ=255)

                       tgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCt                                                                  <  1:3995028/93‑1 (MQ=255)

                       tgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCt                                                                  <  1:4397959/93‑1 (MQ=255)

                          gCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaa                                                               >  2:2293544/1‑93 (MQ=255)

                          gCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaa                                                               <  1:165157/93‑1 (MQ=255)

                            aTCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatg                                                             <  1:3370183/93‑1 (MQ=255)

                            aTCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatg                                                             <  2:479058/93‑1 (MQ=255)

                             tCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatga                                                            >  1:3205362/1‑93 (MQ=255)

                             tCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatga                                                            >  1:1433875/1‑93 (MQ=255)

                             tCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatga                                                            >  1:2375923/1‑93 (MQ=255)

                             tCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatga                                                            <  1:2110646/93‑1 (MQ=255)

                              ccTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatgaa                                                           <  1:3533581/93‑1 (MQ=255)

                              ccTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatgaa                                                           >  1:1266227/1‑93 (MQ=255)

                               cTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAAc                                                          >  1:3599117/1‑93 (MQ=255)

                                tCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAAcc                                                         <  1:1574715/93‑1 (MQ=255)

                                tCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGATGAAAGCCAGCAGCTTGCTGAATGAAcc                                                         >  1:4509784/1‑93 (MQ=255)

                                 cACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCg                                                        >  1:442932/1‑93 (MQ=255)

                                 cACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCg                                                        <  1:2554779/93‑1 (MQ=255)

                                   cccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCa                                                      <  1:1922544/93‑1 (MQ=255)

                                    ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAg                                                     >  2:871591/1‑93 (MQ=255)

                                    ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAg                                                     <  1:2813962/93‑1 (MQ=255)

                                     cTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAgg                                                    >  2:4462211/1‑93 (MQ=255)

                                     cTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAgg                                                    >  2:769558/1‑93 (MQ=255)

                                       gCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGtt                                                  >  2:59860/1‑93 (MQ=255)

                                        cAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTg                                                 >  1:2980838/1‑93 (MQ=255)

                                          aCGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGGTGAt                                               >  1:2849345/1‑93 (MQ=255)

                                            ggTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATAc                                             <  1:3291830/93‑1 (MQ=255)

                                              tACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACtt                                           >  1:2490499/1‑93 (MQ=255)

                                              tACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACtt                                           >  2:2122171/1‑93 (MQ=255)

                                              tACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCCCAGGTTGATACtt                                           >  2:4590797/1‑93 (MQ=255)

                                               aCGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTc                                          <  2:542782/93‑1 (MQ=255)

                                                 gcgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTa                                        <  1:1080202/93‑1 (MQ=255)

                                                  cgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTGTac                                       <  2:3312579/93‑1 (MQ=255)

                                                   gtgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTaca                                      >  2:4227922/1‑93 (MQ=255)

                                                   gtgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTaca                                      >  2:3793267/1‑93 (MQ=255)

                                                     gtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAcc                                    >  2:3936199/1‑93 (MQ=255)

                                                     gtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAcc                                    >  2:1071508/1‑93 (MQ=255)

                                                     attGCAACAGATGGTTTATATCCTGCAGCCAGACATTGCTCTGGTGAAAGCCAGCAGCTTGCTGACTGAACCGCAGGTTGATACTTCTACAcc                                    <  2:4374745/92‑1 (MQ=255)

                                                       tGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgc                                  <  1:2135171/93‑1 (MQ=255)

                                                       tGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgc                                  <  1:3720816/93‑1 (MQ=255)

                                                        gCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgcc                                 >  2:1815314/1‑93 (MQ=255)

                                                          aaCAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGa                               >  1:3432881/1‑93 (MQ=255)

                                                             aGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaa                            <  1:218696/93‑1 (MQ=255)

                                                              gATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaaa                           >  2:3231098/1‑93 (MQ=255)

                                                              gATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaaa                           >  2:132076/1‑93 (MQ=255)

                                                               aTGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAAc                          <  1:994722/93‑1 (MQ=255)

                                                               aTGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAAc                          <  2:1668558/93‑1 (MQ=255)

                                                               aTGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAAc                          <  2:3597388/93‑1 (MQ=255)

                                                               aTGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAAc                          <  1:136937/93‑1 (MQ=255)

                                                                tGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACg                         >  2:4342522/1‑93 (MQ=255)

                                                                tGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACg                         >  2:150161/1‑93 (MQ=255)

                                                                 gttttATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGc                        <  2:4647257/91‑1 (MQ=255)

                                                                 ggTTTATACCCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGc                        <  2:2558800/93‑1 (MQ=255)

                                                                 ggTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGc                        <  2:2015590/93‑1 (MQ=255)

                                                                 ggTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGc                        <  1:4476582/93‑1 (MQ=255)

                                                                  gTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCt                       <  1:831445/93‑1 (MQ=255)

                                                                    ttATACCCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATTCTTCTACACCGCCGAAAAACGCtgt                     <  2:4159990/93‑1 (MQ=255)

                                                                    ttATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgt                     >  2:1768914/1‑93 (MQ=255)

                                                                    ttATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGGTGATACTTCTACACCGCCGAAAAACGCtgt                     >  1:4041870/1‑93 (MQ=255)

                                                                     tataACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgtg                    >  1:266676/1‑93 (MQ=255)

                                                                      ataACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGc                   >  1:2630363/1‑93 (MQ=255)

                                                                       taACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGcc                  >  2:3981075/1‑93 (MQ=255)

                                                                        aaCCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCt                 >  2:1642813/1‑93 (MQ=255)

                                                                           cTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGAt              >  2:832775/1‑93 (MQ=255)

                                                                           cTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGAt              >  2:3157231/1‑93 (MQ=255)

                                                                           cTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGAt              <  1:3753307/93‑1 (MQ=255)

                                                                            tGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATg             <  1:687751/93‑1 (MQ=255)

                                                                                gccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCAcc         <  1:4215111/93‑1 (MQ=255)

                                                                                 ccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcg        >  1:1223137/1‑93 (MQ=255)

                                                                                 ccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcg        >  2:235136/1‑93 (MQ=255)

                                                                                  cAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcgc       >  2:42406/1‑93 (MQ=255)

                                                                                  cAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcgc       >  1:3416793/1‑93 (MQ=255)

                                                                                    gACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGt     >  1:239869/1‑93 (MQ=255)

                                                                                    gACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGt     >  1:2145035/1‑93 (MQ=255)

                                                                                     aCATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTa    >  2:1973798/1‑93 (MQ=255)

                                                                                       aTTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAt  >  2:4270451/1‑93 (MQ=255)

                                                                                       aTTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAt  <  1:195315/93‑1 (MQ=255)

                                                                                       aTTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAt  >  1:1447757/1‑93 (MQ=255)

CGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAT  >  NC_000913/3786325‑3786504

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: `‑`