BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,927,743	T→C	Y366C (TAC→TGC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,927,743	0	T	C	99.1%	379.0 / ‑4.9	107	Y366C (TAC→TGC)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (47/59); total (47/60)
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.94e-01

Only 100 of 107 total aligned reads displayed.

GACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCC  >  NC_000913/1927645‑1927834

gacgCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTa                                                                                            <  2:1597113/100‑1 (MQ=255)

 acgCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAg                                                                                           <  2:2447947/100‑1 (MQ=255)

 acgCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAg                                                                                           <  1:908625/100‑1 (MQ=255)

  cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAgg                                                                                          <  1:3090496/100‑1 (MQ=255)

   gCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGc                                                                                         <  2:1079529/100‑1 (MQ=255)

   gCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGc                                                                                         <  2:661876/100‑1 (MQ=255)

    cTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCa                                                                                        >  2:2709769/1‑100 (MQ=255)

      cACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAAt                                                                                      >  2:1261290/1‑100 (MQ=255)

      cACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAAt                                                                                      >  2:271331/1‑100 (MQ=255)

       aCCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATc                                                                                     >  1:1360613/1‑100 (MQ=255)

          ggTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAg                                                                                  <  1:1867728/100‑1 (MQ=255)

           gTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCCTGCAATCCAga                                                                                 <  1:2812853/100‑2 (MQ=255)

           gTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAgg                                                                                 <  1:2194356/100‑1 (MQ=255)

            tatccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCCTGCAATCCAgaa                                                                                <  1:2372512/100‑3 (MQ=255)

             atccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCCTGCAATCCAgaac                                                                               <  1:616884/100‑4 (MQ=255)

             atccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTc                                                                               <  2:376561/100‑1 (MQ=255)

             atccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTc                                                                               <  2:1731824/100‑1 (MQ=255)

              tccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCCTGCAATCCAgaacc                                                                              <  1:520046/100‑5 (MQ=255)

              tccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCCTGCAATCCAgaacc                                                                              <  2:2367090/100‑5 (MQ=255)

               ccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTcac                                                                             >  2:2132929/1‑100 (MQ=255)

                catatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTcaca                                                                            >  1:1909670/1‑100 (MQ=255)

                catatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTcaca                                                                            >  1:925574/1‑100 (MQ=255)

                catatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTATCAGGTTCTGGATTGCAGGCAATCCAGGTcaca                                                                            >  2:45777/1‑100 (MQ=255)

                  tatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATa                                                                          >  2:242472/1‑100 (MQ=255)

                  tatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATa                                                                          <  2:1533958/100‑1 (MQ=255)

                   atccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAg                                                                         >  2:820609/1‑100 (MQ=255)

                   atccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAg                                                                         >  2:2618767/1‑100 (MQ=255)

                    tccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAgg                                                                        <  2:1661324/100‑1 (MQ=255)

                    tacatttCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAgg                                                                        <  1:327589/95‑1 (MQ=255)

                     ccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGc                                                                       >  2:537657/1‑100 (MQ=255)

                      caGTTCAAACAAAGTGTGTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGcc                                                                      <  1:1829187/100‑1 (MQ=255)

                      caGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGcc                                                                      <  1:2190105/100‑1 (MQ=255)

                       aGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCg                                                                     >  2:1181510/1‑100 (MQ=255)

                        gTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCgg                                                                    <  2:1293650/100‑1 (MQ=255)

                         ttCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGa                                                                   <  1:1812415/100‑1 (MQ=255)

                          tCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGat                                                                  <  2:242170/100‑1 (MQ=255)

                          tCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGat                                                                  <  2:1618897/100‑1 (MQ=255)

                          tCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGat                                                                  >  2:317194/1‑100 (MQ=255)

                           caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGatc                                                                 >  2:2564224/1‑100 (MQ=255)

                           caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGatc                                                                 >  1:2725889/1‑100 (MQ=255)

                            aaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGatca                                                                <  1:1420849/100‑1 (MQ=255)

                            aaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGatca                                                                <  1:492224/100‑1 (MQ=255)

                            aaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGatca                                                                <  1:2021048/100‑1 (MQ=255)

                             aacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCCCATAGGCCGGatcat                                                               >  1:3033028/1‑100 (MQ=255)

                             aacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGatcat                                                               <  1:3082919/100‑1 (MQ=255)

                             aacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGatcat                                                               >  1:1560465/1‑100 (MQ=255)

                             aacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGatcat                                                               <  1:2770475/100‑1 (MQ=255)

                              acaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTTCAGGCAATCCAGGTCACATAGGCCGGatcatc                                                              >  2:1842180/1‑100 (MQ=255)

                              acaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGatcatc                                                              <  2:242297/100‑1 (MQ=255)

                              acaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGatcatc                                                              <  2:564795/100‑1 (MQ=255)

                                  aGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAg                                                          <  2:2907329/100‑1 (MQ=255)

                                  aGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAg                                                          <  2:114387/100‑1 (MQ=255)

                                  aGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAg                                                          <  2:1712163/100‑1 (MQ=255)

                                   gtgtCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAgg                                                         >  1:384340/1‑100 (MQ=255)

                                    tgtCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGc                                                        <  1:242473/100‑1 (MQ=255)

                                      tCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCaa                                                      <  1:2573060/100‑1 (MQ=255)

                                       cTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAAt                                                     >  1:1722933/1‑100 (MQ=255)

                                       cTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAAt                                                     <  2:1256319/100‑1 (MQ=255)

                                          gtgtAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATAcc                                                  >  2:150584/1‑100 (MQ=255)

                                              aGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATg                                              <  2:2370844/100‑1 (MQ=255)

                                              aGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATg                                              >  2:696008/1‑100 (MQ=255)

                                               gTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGa                                             >  2:2040150/1‑100 (MQ=255)

                                               gTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGa                                             >  1:335582/1‑100 (MQ=255)

                                                tCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGAc                                            <  1:820614/100‑1 (MQ=255)

                                                tCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGAc                                            <  1:2618790/100‑1 (MQ=255)

                                                  aTGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACtt                                          >  2:2588229/1‑100 (MQ=255)

                                                   tGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTc                                         >  1:391422/1‑100 (MQ=255)

                                                   nggAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTc                                         >  1:394015/2‑100 (MQ=255)

                                                     gaagaaTAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTccc                                       >  1:3083141/1‑100 (MQ=255)

                                                      aagaaTAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCg                                      <  2:925565/100‑1 (MQ=255)

                                                        gaaTAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCggg                                    >  2:2274914/1‑100 (MQ=255)

                                                           tAACCATAACGCAATCGCGCAGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCt                                 >  1:562228/1‑100 (MQ=255)

                                                            aaCCATAACGCAATCGCGCGGATTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCtt                                >  2:2445698/1‑100 (MQ=255)

                                                             aCCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTg                               <  1:2564247/100‑1 (MQ=255)

                                                              ccATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGc                              <  1:1351337/100‑1 (MQ=255)

                                                                aTAACGCAATCGCCCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCgg                            >  1:2097668/1‑100 (MQ=255)

                                                                  aaCGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGCTCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCggct                          >  1:2521148/1‑98 (MQ=255)

                                                                  aaCGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGtt                          <  2:629876/100‑1 (MQ=255)

                                                                  aaCGCAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGtt                          >  1:2810133/1‑100 (MQ=255)

                                                                     gCAATCGCGCGGTTGCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAAt                       <  2:230150/100‑1 (MQ=255)

                                                                      cAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAAtt                      >  2:145301/1‑100 (MQ=255)

                                                                      cAATCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAAtt                      <  2:1917504/100‑1 (MQ=255)

                                                                         tCGCGCGGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTgg                   <  2:1360601/100‑1 (MQ=255)

                                                                          cgcgcgGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGgc                  <  1:1657797/100‑1 (MQ=255)

                                                                          cgcgcgGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGgc                  <  1:317196/100‑1 (MQ=255)

                                                                           gcgcgGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGgcg                 <  1:1067094/100‑1 (MQ=255)

                                                                           gcgcgGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGgcg                 >  1:2139143/1‑100 (MQ=255)

                                                                           gcgcgGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGgcg                 >  1:1302386/1‑100 (MQ=255)

                                                                              cgGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCaa              <  1:260372/100‑1 (MQ=255)

                                                                              cgGTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCaa              >  2:1921397/1‑100 (MQ=255)

                                                                               ggTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCaaa             >  1:2219410/1‑100 (MQ=255)

                                                                                gTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAAc            >  1:1695487/1‑100 (MQ=255)

                                                                                gTTTCAGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAAc            <  1:944758/100‑1 (MQ=255)

                                                                                     aGGTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGtt       <  2:1886023/100‑1 (MQ=255)

                                                                                      ggTTCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGtt       >  1:686631/1‑99 (MQ=255)

                                                                                      ggTTCTGGATTGCAGACAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTa      <  1:3044605/100‑1 (MQ=255)

                                                                                       gTTCCGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTaa     <  1:2116293/100‑1 (MQ=255)

                                                                                        ttCTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAAc    <  1:1961599/100‑1 (MQ=255)

                                                                                        accTGGATTGCCTGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAAc    >  2:2131155/3‑100 (MQ=255)

                                                                                          cTGGATTGCAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAAccc  <  1:937087/100‑1 (MQ=255)

GACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCC  >  NC_000913/1927645‑1927834

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