BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,246,351	T→C	S68G (AGT→GGT)	lysP ←	lysine‑specific permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,246,351	0	T	C	97.3%	457.7 / 0.4	149	S68G (AGT→GGT)	lysP	lysine‑specific permease
Reads supporting (aligned to +/- strand): ref base T (2/2); new base C (88/57); total (90/59)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 152 total aligned reads displayed.

GCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGG  >  NC_000913/2246261‑2246440

gcgcGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCgg                                                                                         >  2:4111180/1‑93 (MQ=255)

     aaGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGtcatc                                                                                    >  2:4020337/1‑93 (MQ=255)

     aaGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtcatc                                                                                    >  1:2815254/1‑93 (MQ=255)

     aaGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtcatc                                                                                    >  1:1196744/1‑93 (MQ=255)

     aaGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtcatc                                                                                    <  1:4540850/93‑1 (MQ=255)

      aGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtcatca                                                                                   >  2:1722767/1‑93 (MQ=255)

        ccAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAgg                                                                                 >  2:3743048/1‑93 (MQ=255)

        ccAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAgg                                                                                 <  1:3928308/93‑1 (MQ=255)

        ccAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAgg                                                                                 >  1:1306236/1‑93 (MQ=255)

         cAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGa                                                                                >  1:4238411/1‑93 (MQ=255)

         cAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGa                                                                                >  2:781000/1‑93 (MQ=255)

           aaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAg                                                                              >  2:3477026/1‑93 (MQ=255)

            aGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGt                                                                             <  1:1086990/93‑1 (MQ=255)

            aGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGt                                                                             >  1:1773627/1‑93 (MQ=255)

             gCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTa                                                                            >  2:2427677/1‑93 (MQ=255)

             gCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTa                                                                            >  2:1545327/1‑93 (MQ=255)

              ccTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTAACCGAGACCGGTCATCAGGAAGTaa                                                                           >  2:3881195/1‑93 (MQ=255)

                ttcttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAc                                                                         <  1:1882413/93‑1 (MQ=255)

                ttcttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAc                                                                         >  2:1859947/1‑93 (MQ=255)

                ttcttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAc                                                                         >  2:621493/1‑93 (MQ=255)

                   ttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCAt                                                                      <  2:3779305/93‑1 (MQ=255)

                    tcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATc                                                                     <  1:2378064/93‑1 (MQ=255)

                      aaCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAg                                                                   <  2:4223357/93‑1 (MQ=255)

                       aCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCTAGACCGGTCATCAGGAAGTAAACCATCAgg                                                                  >  2:4742461/1‑93 (MQ=255)

                        cATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGc                                                                 <  1:586157/93‑1 (MQ=255)

                        cATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGc                                                                 <  1:3512537/93‑1 (MQ=255)

                        cATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGc                                                                 <  1:1811509/93‑1 (MQ=255)

                          tAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCa                                                               >  1:3507441/1‑93 (MQ=255)

                            gTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAAt                                                             >  2:134536/1‑93 (MQ=255)

                            gTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAAt                                                             >  2:3012385/1‑93 (MQ=255)

                             ttCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGTAAGTAAACCATCAGGCCAATc                                                            <  1:549065/93‑1 (MQ=255)

                              tCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCa                                                           <  2:830263/93‑1 (MQ=255)

                              tCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCa                                                           <  2:2994977/93‑1 (MQ=255)

                              tCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCa                                                           <  1:1514257/93‑1 (MQ=255)

                               cTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAg                                                          <  1:272386/93‑1 (MQ=255)

                                tGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGc                                                         >  1:308384/1‑93 (MQ=255)

                                  aCCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCat                                                       <  1:321565/93‑1 (MQ=255)

                                  aCCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCat                                                       <  2:4670104/93‑1 (MQ=255)

                                  aCCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCat                                                       >  1:3211762/1‑93 (MQ=255)

                                   ccGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCata                                                      <  1:3796628/93‑1 (MQ=255)

                                   ccGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCata                                                      >  2:4136369/1‑93 (MQ=255)

                                   ccGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCata                                                      >  2:219614/1‑93 (MQ=255)

                                    cGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCatat                                                     <  1:1669353/93‑1 (MQ=255)

                                    cGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCatat                                                     >  1:2238254/1‑93 (MQ=255)

                                    cGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCatat                                                     >  1:4404549/1‑93 (MQ=255)

                                    cGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCatat                                                     <  2:505224/93‑1 (MQ=255)

                                      tAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATAc                                                   <  2:2838538/93‑1 (MQ=255)

                                      tAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATAc                                                   >  1:2885094/1‑93 (MQ=255)

                                        aGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACga                                                 <  1:4648820/93‑1 (MQ=255)

                                        aGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACga                                                 >  1:1691699/1‑93 (MQ=255)

                                         gTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACgag                                                <  1:1076483/93‑1 (MQ=255)

                                         gTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACAGAGACAGGTCATCAGGAAGGAAACCATCAAGCCAATCAGCATATACgag                                                >  2:3328998/1‑93 (MQ=255)

                                          tGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACgaga                                               <  2:4681846/93‑1 (MQ=255)

                                            gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGc                                             >  2:272162/1‑93 (MQ=255)

                                            gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGc                                             <  1:2388086/93‑1 (MQ=255)

                                             cAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCa                                            >  1:934077/1‑93 (MQ=255)

                                              aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCaa                                           <  1:4453111/93‑1 (MQ=255)

                                               aaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCaag                                          >  1:1361067/1‑92 (MQ=255)

                                               aaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAAt                                          >  1:1233485/1‑93 (MQ=255)

                                               aaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAAt                                          >  1:1977697/1‑93 (MQ=255)

                                               aaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAAt                                          >  1:1977698/1‑93 (MQ=255)

                                               aaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAAt                                          <  2:1692355/93‑1 (MQ=255)

                                               aaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAAt                                          <  1:2976663/93‑1 (MQ=255)

                                                   aaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGccc                                      >  2:4364153/1‑93 (MQ=255)

                                                   aaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGccc                                      <  2:801370/93‑1 (MQ=255)

                                                    aCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGcccc                                     <  1:2345068/93‑1 (MQ=255)

                                                       ggAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgcc                                  <  2:1313878/93‑1 (MQ=255)

                                                          aaCCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCgg                               >  2:2452456/1‑93 (MQ=255)

                                                           aCCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGa                              >  2:138651/1‑93 (MQ=255)

                                                            ccGGCATATATGCAGCCAGTTCACCGATACCGGTCATCAGGAAGTAAACCATCAAGCCAATCAGCATATACGAGAGCAATCGCCAGCCCGGAc                             >  2:3570069/1‑93 (MQ=255)

                                                            ccGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAc                             >  2:3432641/1‑93 (MQ=255)

                                                            ccGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAc                             >  1:2259186/1‑93 (MQ=255)

                                                             cGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcc                            <  2:2079918/93‑1 (MQ=255)

                                                               gCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctg                          >  1:368233/1‑93 (MQ=255)

                                                               gCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctg                          >  2:3401234/1‑93 (MQ=255)

                                                                cATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgc                         >  2:1133847/1‑93 (MQ=255)

                                                                cATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgc                         >  1:119138/1‑93 (MQ=255)

                                                                  tatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcct                       >  1:592726/1‑93 (MQ=255)

                                                                    tttGCAGCCAGTTCACCGAGACCGATCACCAGGAAGTAAACCATCAGGCCAGTCAGCATATACGAGAGACACGCCCCGCGCGGACCTGCTTga                     <  2:4193176/91‑1 (MQ=255)

                                                                    tatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTga                     >  1:4649428/1‑93 (MQ=255)

                                                                    tatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTga                     <  1:1396635/93‑1 (MQ=255)

                                                                         aGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAAt                >  1:2021653/1‑93 (MQ=255)

                                                                          gCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATc               <  2:3761016/93‑1 (MQ=255)

                                                                            cAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGt             >  2:2797140/1‑93 (MQ=255)

                                                                              gTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTg           >  2:4661780/1‑93 (MQ=255)

                                                                              gTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTg           <  2:2556113/93‑1 (MQ=255)

                                                                              gTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTg           <  2:4318422/93‑1 (MQ=255)

                                                                                 cACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTgcgc        >  1:3806156/1‑93 (MQ=255)

                                                                                 cACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTgcgc        >  1:1262375/1‑93 (MQ=255)

                                                                                  aCCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGcc       <  2:1926952/93‑1 (MQ=255)

                                                                                  aCCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGcc       >  2:3930268/1‑93 (MQ=255)

                                                                                  aCCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACATGCCTGAGAAATCGTTGCGcc       >  2:1521993/1‑93 (MQ=255)

                                                                                   ccGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCa      <  1:4050814/93‑1 (MQ=255)

                                                                                   ccGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCa      >  1:2561490/1‑93 (MQ=255)

                                                                                    cGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCag     >  1:4362758/1‑93 (MQ=255)

                                                                                     gagaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCaga    >  1:980223/1‑93 (MQ=255)

                                                                                     gagaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCaga    >  1:4146973/1‑93 (MQ=255)

                                                                                      agaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCagag   >  1:3682793/1‑93 (MQ=255)

                                                                                       gaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAgg  <  1:1907393/93‑1 (MQ=255)

                                                                                       gaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAgg  >  1:2091437/1‑93 (MQ=255)

GCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGG  >  NC_000913/2246261‑2246440

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