BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,470,958	0	G	T	21.2%	457.3 / 105.5	199	intergenic (‑45/+28)	ridA/pyrI	enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand): ref base G (78/78); new base T (18/24); total (96/103)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 208 total aligned reads displayed.

GTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGG  >  NC_000913/4470859‑4471056

gTACGGAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAgg                                                                                                    <  2:452846/100‑1 (MQ=255)

 tACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAggg                                                                                                   <  2:1545774/100‑1 (MQ=255)

     tAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTGTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGcc                                                                                               >  1:3105406/1‑100 (MQ=255)

      aaGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCa                                                                                              >  2:3218358/1‑100 (MQ=255)

       aGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAg                                                                                             <  1:1696932/100‑1 (MQ=255)

        ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGa                                                                                            <  1:1447877/100‑1 (MQ=38)

        ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGa                                                                                            <  1:161801/100‑1 (MQ=38)

        ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTCTAGGGAGCCAGa                                                                                            >  2:4095715/1‑100 (MQ=255)

        ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGa                                                                                            >  1:1073778/1‑100 (MQ=255)

        ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGa                                                                                            <  1:2761581/100‑1 (MQ=255)

        ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGa                                                                                            >  1:2025842/1‑100 (MQ=255)

         gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  1:852444/100‑1 (MQ=38)

         gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  2:710660/100‑1 (MQ=38)

         gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  1:2719091/100‑1 (MQ=38)

         gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                           <  1:1919001/100‑1 (MQ=255)

           ncGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGTAGCCAGACtt                                                                                         >  2:3492113/2‑100 (MQ=255)

           ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAATAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  1:2524453/100‑1 (MQ=25)

           ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                         <  1:1837672/100‑1 (MQ=38)

           ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                         >  1:3221901/1‑100 (MQ=255)

           ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                         <  1:728356/100‑1 (MQ=255)

           ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTAATTATAGGGAGCCAGACtt                                                                                         <  2:2650619/100‑1 (MQ=255)

            cGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGACttt                                                                                        <  1:1251639/100‑1 (MQ=38)

             gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                       >  2:2221097/1‑100 (MQ=255)

              aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGACttttt                                                                                      <  1:3606761/100‑2 (MQ=37)

              aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                      <  2:2087483/100‑1 (MQ=255)

              aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                      >  2:1908494/1‑100 (MQ=255)

               tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                                     <  2:1632910/100‑1 (MQ=255)

                aGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGACtttttac                                                                                    <  2:3160161/100‑4 (MQ=37)

                 gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                                   >  1:2233174/1‑100 (MQ=255)

                 gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                                   <  1:2181236/100‑1 (MQ=255)

                     ccGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCaa                                                                               <  2:819877/100‑1 (MQ=255)

                      cGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAc                                                                              <  1:4150764/100‑1 (MQ=255)

                        gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                            <  1:1786710/100‑1 (MQ=255)

                         tGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGc                                                                           >  2:2395909/1‑100 (MQ=255)

                            aTTTTCCGTCGCGATCGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAACTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaat                                                                        >  2:3752875/1‑100 (MQ=255)

                             ttttCCGTCGCGATAGTTTTGCTCATGATGTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatt                                                                       <  1:4178221/100‑1 (MQ=255)

                              tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatta                                                                      >  1:4231326/1‑100 (MQ=255)

                              tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatta                                                                      >  2:4312925/1‑100 (MQ=255)

                              tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatta                                                                      >  1:2159830/1‑100 (MQ=255)

                                  cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTg                                                                  <  2:2016135/100‑1 (MQ=255)

                                   gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                 >  1:2381026/1‑100 (MQ=255)

                                        gATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGc                                                            >  2:3649732/1‑100 (MQ=255)

                                         aTAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGca                                                           <  2:1883114/100‑1 (MQ=255)

                                          tAGTTTTGCTCATGATTTCTCCGTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                          <  1:1536280/100‑1 (MQ=255)

                                              tttGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACcaca                                                      <  2:3974057/100‑1 (MQ=255)

                                               ttGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                                     <  2:2415551/100‑1 (MQ=255)

                                                 gCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATCGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTa                                                   <  2:4065114/100‑1 (MQ=255)

                                                   tCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                 >  1:4416856/1‑100 (MQ=255)

                                                   tCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                 >  1:139138/1‑100 (MQ=255)

                                                   tCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                 <  1:2129722/100‑1 (MQ=255)

                                                    cATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                >  2:561753/1‑100 (MQ=255)

                                                    cATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                >  1:404287/1‑100 (MQ=255)

                                                     aTGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATggg                                               <  2:2971070/100‑1 (MQ=255)

                                                       gATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                             >  1:2917203/1‑100 (MQ=255)

                                                        atgtctcCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                            >  2:3397559/4‑100 (MQ=255)

                                                         tttCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaaa                                           >  2:2484764/1‑100 (MQ=255)

                                                               cTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                     <  1:2348609/100‑1 (MQ=255)

                                                               cTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                     <  2:615839/100‑1 (MQ=255)

                                                               cTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                     >  1:652431/1‑100 (MQ=255)

                                                                tttATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                                    <  1:119498/100‑1 (MQ=255)

                                                                tttATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                                    <  1:164486/100‑1 (MQ=255)

                                                                 ttattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttttt                                   <  1:2645435/100‑1 (MQ=255)

                                                                 ttattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttttt                                   >  2:1432009/1‑100 (MQ=255)

                                                                  tattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                  <  1:1489580/100‑1 (MQ=255)

                                                                    ttaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                <  1:367154/100‑1 (MQ=255)

                                                                    ttaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                <  1:4529789/100‑1 (MQ=255)

                                                                       cAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGt                             >  2:2791857/1‑100 (MQ=255)

                                                                         gCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGCCTTTTATTACCAACCGCAATTAATTGGCCAACACAACATTATGTGAAAACTCTGTTTCAAAGTAt                           >  2:3532551/1‑100 (MQ=255)

                                                                           ggTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAttt                         >  2:4434740/1‑100 (MQ=255)

                                                                            gtaataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                          >  1:4519794/6‑100 (MQ=25)

                                                                            gTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                        <  1:3764251/100‑1 (MQ=255)

                                                                            gTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                        >  1:1847178/1‑100 (MQ=255)

                                                                             taataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAttt                         >  2:3574772/5‑100 (MQ=25)

                                                                             tA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                       <  2:3070032/100‑1 (MQ=255)

                                                                             tA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                       >  2:1160046/1‑100 (MQ=255)

                                                                             tA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                       <  1:3971062/100‑1 (MQ=255)

                                                                              aataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                        >  1:1312454/4‑100 (MQ=25)

                                                                              aataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                        >  2:725931/4‑100 (MQ=25)

                                                                               ataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                       >  1:2499847/3‑100 (MQ=37)

                                                                                taaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtttaca                      >  2:3057152/2‑97 (MQ=37)

                                                                                taaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCa                      >  2:2918754/2‑100 (MQ=37)

                                                                                 aaaaGTCTGGCTCATTATAGGGAGCCCGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATTGGAAAACTCTTTTTCACAGTATTTGCAt                     <  1:3440281/100‑1 (MQ=255)

                                                                                 aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAt                     <  1:4301965/100‑1 (MQ=255)

                                                                                 aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAt                     <  2:1055268/100‑1 (MQ=255)

                                                                                  aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAAGTGGCCAGCACCAGATTATGGGAAAACTCTTTTTCACAGTAGTTGCAtt                    >  2:933691/1‑100 (MQ=37)

                                                                                  aaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt                    <  1:3060855/100‑1 (MQ=255)

                                                                                  aaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATAAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt                    <  1:2816706/100‑1 (MQ=255)

                                                                                   aaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt                   <  1:2449105/100‑1 (MQ=255)

                                                                                     gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGa                 >  1:2103837/1‑100 (MQ=38)

                                                                                     gTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGa                 >  2:1915434/1‑100 (MQ=255)

                                                                                      tCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg                <  1:4349207/100‑1 (MQ=255)

                                                                                       cTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc               >  1:3204359/1‑100 (MQ=38)

                                                                                         ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc             >  2:2267413/1‑100 (MQ=255)

                                                                                         ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc             >  2:1257057/1‑100 (MQ=255)

                                                                                          gCTCATTATAGTGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg            <  2:2611405/100‑1 (MQ=255)

                                                                                           cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa           >  1:3826318/1‑100 (MQ=255)

                                                                                           cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa           >  1:1068269/1‑100 (MQ=255)

                                                                                           cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGAGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa           >  2:2343640/1‑100 (MQ=255)

                                                                                                    gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg  >  1:3186167/1‑100 (MQ=255)

                                                                                                    gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg  >  1:3299083/1‑100 (MQ=255)

GTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGG  >  NC_000913/4470859‑4471056

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