BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,788,045	A→G	S408S (AGT→AGC)	b1706 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,788,045	0	A	G	96.4%	337.4 / 1.3	112	S408S (AGT→AGC)	b1706	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/3); new base G (60/48); total (61/51)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 119 total aligned reads displayed.

ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGC  >  NC_000913/1787955‑1788135

aCCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGAcgctcgct                                                                                          >  2:239865/1‑93 (MQ=255)

 cccTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTa                                                                                         >  1:2465511/1‑93 (MQ=255)

 cccTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTa                                                                                         <  1:4109750/93‑1 (MQ=255)

  ccTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTAt                                                                                        <  2:50648/93‑1 (MQ=255)

  ccTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTAt                                                                                        <  1:2065187/93‑1 (MQ=255)

    ttttCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCg                                                                                      >  2:3067647/1‑93 (MQ=255)

     tttCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGc                                                                                     >  1:346341/1‑93 (MQ=255)

       tCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCtcgctatcgcta                                                                                   <  2:1100768/93‑1 (MQ=255)

        cTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATGCCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTaa                                                                                  <  1:1989368/93‑1 (MQ=255)

         tGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAAc                                                                                 >  1:3163011/1‑93 (MQ=255)

          gCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAAcc                                                                                <  2:3340778/93‑1 (MQ=255)

          gCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAAcc                                                                                <  2:3478287/93‑1 (MQ=255)

              cTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGt                                                                            >  1:3157780/1‑93 (MQ=255)

              cTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGt                                                                            <  2:3044657/93‑1 (MQ=255)

               ttCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTc                                                                           >  2:66019/1‑93 (MQ=255)

               ttCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTc                                                                           >  1:3806556/1‑93 (MQ=255)

                tCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCt                                                                          <  2:587095/93‑1 (MQ=255)

                 cAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCtt                                                                         <  1:2613252/93‑1 (MQ=255)

                 cAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCtt                                                                         >  1:82417/1‑93 (MQ=255)

                 cAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCCGACGCTCGCTATCGCTAACCTCGTCtt                                                                         >  1:2401013/1‑93 (MQ=255)

                  aaTCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCttg                                                                        <  1:1208104/93‑1 (MQ=255)

                   aTCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCttgt                                                                       >  1:713143/1‑93 (MQ=255)

                    tCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCttgtt                                                                      >  1:2281715/1‑93 (MQ=255)

                    tCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCttgtt                                                                      <  2:39201/93‑1 (MQ=255)

                     cGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCttgttg                                                                     <  1:4652604/93‑1 (MQ=255)

                      gCCCGTTGCGCCAACCAGGTGCGCAATACCAGAGCGGGGCTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGc                                                                    >  2:3648182/1‑93 (MQ=255)

                        ccGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaa                                                                  <  2:4638748/93‑1 (MQ=255)

                        ccGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaa                                                                  >  1:1650108/1‑93 (MQ=255)

                        ccGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaa                                                                  >  1:3945428/1‑93 (MQ=255)

                        ccGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaa                                                                  >  2:2102628/1‑93 (MQ=255)

                         cGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaaa                                                                 <  1:4686356/93‑1 (MQ=255)

                         cGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaaa                                                                 <  1:4686357/93‑1 (MQ=255)

                         cGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCaaa                                                                 >  2:3987089/1‑93 (MQ=255)

                          gtagcgcCAACCAGTTGCGCAATTCCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCTTCTTGTTGCAAAc                                                                <  2:4327266/90‑1 (MQ=255)

                          gTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAAc                                                                <  1:1326850/93‑1 (MQ=255)

                             gcgcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGcc                                                             <  1:4035146/93‑1 (MQ=255)

                              cgcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGccc                                                            <  1:3008649/93‑1 (MQ=255)

                              cgcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCACTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGccc                                                            >  2:1124868/1‑93 (MQ=255)

                               gcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGcccc                                                           <  2:2178404/93‑1 (MQ=255)

                               gcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGcccc                                                           <  1:2519809/93‑1 (MQ=255)

                               gcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGcccc                                                           >  2:2141349/1‑93 (MQ=255)

                               gcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGcccc                                                           >  1:2957916/1‑93 (MQ=255)

                                ccAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGccccc                                                          >  1:3539701/1‑93 (MQ=255)

                                ccAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGccccc                                                          <  1:2877459/93‑1 (MQ=255)

                                  aaCCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGa                                                        >  1:2309654/1‑93 (MQ=255)

                                  aaCCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGa                                                        >  2:2656894/1‑93 (MQ=255)

                                    ccAGTTGCGCAATCCCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTCTCGCTAACCTCGTCTTGTTGCAAACGCCCCCCATa                                                      >  1:810546/1‑93 (MQ=255)

                                    ccAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATa                                                      >  1:810547/1‑93 (MQ=255)

                                    ccAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATa                                                      >  2:1859728/1‑93 (MQ=255)

                                    ccAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATa                                                      >  2:1859725/1‑93 (MQ=255)

                                       gTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACg                                                   <  1:2737013/93‑1 (MQ=255)

                                            gcAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCaa                                              >  1:3928424/1‑93 (MQ=255)

                                            gcAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCaa                                              >  2:3048490/1‑93 (MQ=255)

                                              aaTACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAc                                            >  2:2357736/1‑93 (MQ=255)

                                              aaTACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAc                                            >  1:4569313/1‑93 (MQ=255)

                                              aaTACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAc                                            >  1:1390827/1‑93 (MQ=255)

                                              aaTACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAc                                            >  1:1768913/1‑93 (MQ=255)

                                                tACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCa                                          >  2:1104708/1‑93 (MQ=255)

                                                tACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCa                                          >  2:3854426/1‑93 (MQ=255)

                                                 aCCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAg                                         >  2:690945/1‑93 (MQ=255)

                                                   cAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGtc                                       >  1:1958130/1‑93 (MQ=255)

                                                    agagCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGtca                                      <  2:1293119/93‑1 (MQ=255)

                                                     gagCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGtcat                                     >  1:2821823/1‑93 (MQ=255)

                                                     gagCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGtcat                                     <  2:654343/93‑1 (MQ=255)

                                                       gCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGtcatca                                   >  2:4609622/1‑93 (MQ=255)

                                                         tgggTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCaaa                                 >  2:3617225/2‑93 (MQ=255)

                                                         ggggTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCaaa                                 <  1:3897064/93‑1 (MQ=255)

                                                          gtgTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAAt                                <  2:4285778/91‑1 (MQ=255)

                                                          gggTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAAt                                <  1:1697328/93‑1 (MQ=255)

                                                           ggTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATg                               <  2:1364536/93‑1 (MQ=255)

                                                            gTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGc                              <  1:4605302/93‑1 (MQ=255)

                                                             ttAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGcc                             >  2:4735560/1‑93 (MQ=255)

                                                             ttAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGcc                             >  1:1931638/1‑93 (MQ=255)

                                                                 cGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCAc                         <  2:526407/93‑1 (MQ=255)

                                                                  gCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACg                        >  2:2346575/1‑93 (MQ=255)

                                                                   cTTTGCATCAGTTGCTGGCGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGa                       <  2:220286/93‑1 (MQ=255)

                                                                   cTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGa                       <  2:1344109/93‑1 (MQ=255)

                                                                   cTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGa                       <  1:2973917/93‑1 (MQ=255)

                                                                     ttGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATc                     <  2:2306188/93‑1 (MQ=255)

                                                                      tGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCa                    >  1:1679229/1‑93 (MQ=255)

                                                                      tGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCa                    >  1:2709964/1‑93 (MQ=255)

                                                                        cATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAAt                  <  2:1239960/93‑1 (MQ=255)

                                                                         aTCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATa                 <  1:2455028/93‑1 (MQ=255)

                                                                          tCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATaa                >  1:2990263/1‑93 (MQ=255)

                                                                          tCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATaa                >  2:2118235/1‑93 (MQ=255)

                                                                           cAGTTGCAGACGCTCGCTATCGCTAACCTCGTCGTGTTGCAAACGCCTCCGATAACGGGCAACCCAGTCATCAAATGCCGCACGATCAATaaa               <  2:1371/93‑1 (MQ=255)

                                                                            aGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAAc              >  1:3525419/1‑93 (MQ=255)

                                                                            aGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAAc              <  2:673471/93‑1 (MQ=255)

                                                                              ttGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACtc            <  1:1026765/93‑1 (MQ=255)

                                                                              ttGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACtc            >  2:4471426/1‑93 (MQ=255)

                                                                              ttGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACtc            >  2:2917223/1‑93 (MQ=255)

                                                                                gCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCAGCACGATCAATAAACtcat          >  1:1436927/1‑93 (MQ=255)

                                                                                  tGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATCAACtcataa        >  2:1858437/1‑91 (MQ=255)

                                                                                  tGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACtcatca        <  1:1307553/93‑1 (MQ=255)

                                                                                   gACGCTCGCTCTCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCAc       >  2:4347789/1‑93 (MQ=255)

                                                                                     cgctcgctATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGt     >  1:208304/1‑93 (MQ=255)

                                                                                     cgctcgctATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGt     <  1:2793508/93‑1 (MQ=255)

                                                                                      gCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTa    <  2:3896453/93‑1 (MQ=255)

                                                                                      ccTCGCTATCGCTAACCTCGTCTTTTTGCAAACGCCCCCGATAACGGGCCAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTa    <  1:1095027/92‑1 (MQ=255)

                                                                                        tCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGc  <  1:649534/93‑1 (MQ=255)

ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGC  >  NC_000913/1787955‑1788135

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