BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,210,743	A→G	Q19R (CAA→CGA)	rpoD →	RNA polymerase sigma(70)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,210,743	0	A	G	99.6%	812.5 / ‑4.9	257	Q19R (CAA→CGA)	rpoD	RNA polymerase sigma(70)
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (124/132); total (124/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 262 total aligned reads displayed.

CACCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCA  >  NC_000913/3210653‑3210832

 acCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAgg                                                                                        >  2:1088944/1‑93 (MQ=255)

 acCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAgg                                                                                        >  1:1474465/1‑93 (MQ=255)

  ccACCCTCATGAAATAAGTGTGTATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGc                                                                                       <  2:1765757/93‑1 (MQ=255)

  ccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGc                                                                                       <  1:4522153/93‑1 (MQ=255)

  ccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGc                                                                                       >  1:4095213/1‑93 (MQ=255)

  ccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGACACAGCTGAAACCTCTTGTCACTCGTGGTAAGGAGCGAGGc                                                                                       >  2:3230429/1‑93 (MQ=255)

   cAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCt                                                                                      >  1:159040/1‑93 (MQ=255)

   cAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCt                                                                                      >  1:3687325/1‑93 (MQ=255)

    aaCCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTa                                                                                     <  2:3973072/93‑1 (MQ=255)

    aaCCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTa                                                                                     >  1:1275285/1‑93 (MQ=255)

     aCCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTAt                                                                                    >  1:2291610/1‑93 (MQ=255)

      ccTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATc                                                                                   <  2:870884/93‑1 (MQ=255)

       cTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCt                                                                                  <  2:2174743/93‑1 (MQ=255)

        tCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTg                                                                                 >  1:4282638/1‑93 (MQ=255)

        tCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTg                                                                                 <  2:1174346/93‑1 (MQ=255)

         cATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGa                                                                                <  2:835064/93‑1 (MQ=255)

         cATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGa                                                                                <  2:2511039/93‑1 (MQ=255)

           tGAAATAAGTGTGGATACCGTCTTATGGAGCAGAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGAcc                                                                              >  1:4550163/1‑93 (MQ=255)

           tGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGAcc                                                                              <  1:2864383/93‑1 (MQ=255)

              aaTAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTAt                                                                           <  1:4010047/93‑1 (MQ=255)

               aTAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATg                                                                          >  2:3510055/1‑93 (MQ=255)

               aTAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATg                                                                          <  2:928539/93‑1 (MQ=255)

                tAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGc                                                                         <  1:4595729/93‑1 (MQ=255)

                tAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGc                                                                         >  2:3668671/1‑93 (MQ=255)

                tAAGAGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGc                                                                         <  1:575688/93‑1 (MQ=255)

                   gtgtGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGa                                                                      <  1:2025405/93‑1 (MQ=255)

                       ggATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTc                                                                  <  1:348015/93‑1 (MQ=255)

                        gATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCa                                                                 >  2:3492283/1‑93 (MQ=255)

                         aTACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCaa                                                                >  1:1826600/1‑93 (MQ=255)

                          tACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAAt                                                               >  1:1437018/1‑93 (MQ=255)

                            ccGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGa                                                             <  2:4488171/93‑1 (MQ=255)

                             cGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGAc                                                            >  2:1156442/1‑93 (MQ=255)

                               tCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCa                                                          >  2:2532087/1‑93 (MQ=255)

                               tCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCa                                                          >  1:178658/1‑93 (MQ=255)

                                 ttATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATc                                                        <  1:2401978/93‑1 (MQ=255)

                                 ttATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATc                                                        >  1:3605099/1‑93 (MQ=255)

                                 ttATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATc                                                        >  1:1636762/1‑93 (MQ=255)

                                  tATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCt                                                       <  1:3740185/93‑1 (MQ=255)

                                    tGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGc                                                     >  2:289338/1‑93 (MQ=255)

                                      gAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCg                                                   >  1:1834397/1‑93 (MQ=255)

                                      gAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCg                                                   >  2:4552121/1‑93 (MQ=255)

                                          aaaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAg                                               >  2:1933434/1‑93 (MQ=255)

                                          aaaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAg                                               >  2:1854116/1‑93 (MQ=255)

                                          aaaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAg                                               >  2:1637733/1‑93 (MQ=255)

                                          aaaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAg                                               >  2:1567080/1‑93 (MQ=255)

                                           aaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGa                                              <  1:2355838/93‑1 (MQ=255)

                                           aaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGa                                              <  1:2633712/93‑1 (MQ=255)

                                           aaaCCCGCAGTCACAGCTGAAACTTCTTGGCACGCGTGGGAAGGAGCGAGGCTATCTGACCTATTCCGAGGACAATGACCATCTTCCGGAAGa                                              >  2:2030576/1‑93 (MQ=255)

                                            aaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGat                                             <  1:729156/93‑1 (MQ=255)

                                             aCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGata                                            <  2:514281/93‑1 (MQ=255)

                                             aCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGata                                            <  1:2262517/93‑1 (MQ=255)

                                             aCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGata                                            <  1:4460518/93‑1 (MQ=255)

                                              cccGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGatat                                           <  1:3007181/93‑1 (MQ=255)

                                               ccGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtc                                          >  1:2245802/1‑93 (MQ=255)

                                               ccGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtc                                          <  1:4372871/93‑1 (MQ=255)

                                               ccGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtc                                          >  1:1024816/1‑93 (MQ=255)

                                                 gCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtcgt                                        <  2:760217/93‑1 (MQ=255)

                                                 gCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtcgt                                        >  2:4659086/1‑93 (MQ=255)

                                                    gTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGa                                     >  2:4085896/1‑93 (MQ=255)

                                                     tCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGAt                                    <  1:4034053/93‑1 (MQ=255)

                                                      cacaGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGAtt                                   >  1:2133808/1‑93 (MQ=255)

                                                        caGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtca                                 <  2:3369248/93‑1 (MQ=255)

                                                          gCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcaga                               >  2:2330342/1‑93 (MQ=255)

                                                          gCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcaga                               <  2:1283651/93‑1 (MQ=255)

                                                            tGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatc                             >  2:3934270/1‑93 (MQ=255)

                                                            tGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatc                             <  2:3341705/93‑1 (MQ=255)

                                                             gAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatca                            >  1:3295021/1‑93 (MQ=255)

                                                             gAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatca                            >  1:656117/1‑93 (MQ=255)

                                                             gAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatca                            >  1:2181643/1‑93 (MQ=255)

                                                             gAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatca                            <  2:2697932/93‑1 (MQ=255)

                                                               aaCTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatcaga                          <  2:1905205/93‑1 (MQ=255)

                                                                aCTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTcagatcagat                         <  2:4274412/93‑1 (MQ=255)

                                                                aCTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTcagatcagat                         <  1:947228/93‑1 (MQ=255)

                                                                aCTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTcagatcagat                         <  2:704356/93‑1 (MQ=255)

                                                                 cttcttGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCagatcagatc                        >  1:3055736/1‑93 (MQ=255)

                                                                 cttcttGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCagatcagatc                        <  1:397384/93‑1 (MQ=255)

                                                                  ttcttGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCagatcagatca                       >  1:977687/1‑92 (MQ=255)

                                                                  ttcttGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCg                       >  2:1431320/1‑93 (MQ=255)

                                                                    tttGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGaa                     <  2:2415548/92‑1 (MQ=255)

                                                                     ttGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAg                    <  2:1242615/93‑1 (MQ=255)

                                                                     ttGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCCTCGATTCAGATCAGATCGAAg                    <  1:4647641/93‑1 (MQ=255)

                                                                     ttGTCACCCGTGGTAAGGAGCGAGGATATCTGACCTATGCCGAGGCCAATGACCATCTGCCGGAAGATATCGTCGATTGAGATCAGATCAAAg                    >  2:2079144/1‑93 (MQ=255)

                                                                       gTCACCCGTGGTAAGTAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAc                  <  2:3141009/93‑1 (MQ=255)

                                                                        tCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAca                 <  1:888379/93‑1 (MQ=255)

                                                                         cACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAcat                <  2:4424573/93‑1 (MQ=255)

                                                                           cccGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAcatca              <  2:3764024/93‑1 (MQ=255)

                                                                           cccGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAcatca              <  2:1776633/93‑1 (MQ=255)

                                                                            ccGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAcatcat             >  1:1326880/1‑93 (MQ=255)

                                                                             cGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACatcatc            <  2:2837084/93‑1 (MQ=255)

                                                                             cGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACatcatc            <  1:3580337/93‑1 (MQ=255)

                                                                              gTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATcc           >  1:1037019/1‑93 (MQ=255)

                                                                                 gTAAGGAGCGAGGGTATCTGACCTATGCGGAGGCCAATGACCATCTGCCGGAAGCTATCGTCGATTCACATCAGATCAAAGACATCATCCaat        >  2:661216/1‑92 (MQ=255)

                                                                                 gTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCaaa        >  1:1669964/1‑93 (MQ=255)

                                                                                 gTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCaaa        <  1:3154617/93‑1 (MQ=255)

                                                                                 gTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCaaa        <  2:2294635/93‑1 (MQ=255)

                                                                                    aGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGa     <  2:544906/93‑1 (MQ=255)

                                                                                     ggAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGAt    >  1:2517128/1‑93 (MQ=255)

                                                                                     ggAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGAt    <  2:2995192/93‑1 (MQ=255)

                                                                                      gAGCGAGGCTATCTGACGTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATc   <  2:4365095/93‑1 (MQ=255)

                                                                                      gAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATc   <  2:4134457/93‑1 (MQ=255)

CACCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCA  >  NC_000913/3210653‑3210832

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