BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	654,812	G→A	A127V (GCC→GTC)	dcuC ←	DcuC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	654,812	0	G	A	97.6%	438.7 / ‑1.7	127	A127V (GCC→GTC)	dcuC	DcuC
Reads supporting (aligned to +/- strand): ref base G (1/2); new base A (61/62); total (62/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

Only 100 of 135 total aligned reads displayed.

TACCAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGGCGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCATTCGCG  >  NC_000913/654723‑654901

tACCAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGa                                                                                        >  2:2484725/1‑93 (MQ=255)

 aCCAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAgacgac                                                                                       <  2:1469871/93‑1 (MQ=255)

     aCGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACaaaa                                                                                   >  1:2833010/1‑93 (MQ=255)

     aCGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACaaaa                                                                                   >  2:1570670/1‑93 (MQ=255)

       gTTTACCATCACCGGAAATAGGGTTGCCATCAGCAACACACCCATACCCGTTGCGTTACAGACGGCCAGAGACATCAGACAGACGACAAAATa                                                                                 >  2:4252061/1‑93 (MQ=255)

        tttACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAg                                                                                >  1:1446061/1‑93 (MQ=255)

         ttACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAgg                                                                               <  2:2456751/93‑1 (MQ=255)

         ttACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAgg                                                                               >  1:3965077/1‑93 (MQ=255)

         ttACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAgg                                                                               <  1:1321305/93‑1 (MQ=255)

         ttACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCACACAGACGACAAAATAgg                                                                               >  1:2393791/1‑93 (MQ=255)

            ccATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGGCGACAAAATAggcgg                                                                            <  2:2128280/93‑1 (MQ=255)

             cATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAggcggc                                                                           >  2:3404943/1‑93 (MQ=255)

              aTCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGTCGGCCAGAGACATCAGACAGACGACAAAATATGCGGCa                                                                          >  2:2233498/1‑93 (MQ=255)

               tCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCaa                                                                         <  1:2455241/93‑1 (MQ=255)

                cACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAat                                                                        >  1:1776730/1‑93 (MQ=255)

                 aCCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatc                                                                       <  1:4228618/93‑1 (MQ=255)

                  ccGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatca                                                                      <  1:2868629/93‑1 (MQ=255)

                   cGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatcat                                                                     <  2:4385812/93‑1 (MQ=255)

                    ggAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatcatc                                                                    >  1:3341771/1‑93 (MQ=255)

                    ggAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatcatc                                                                    >  1:3773517/1‑93 (MQ=255)

                    ggAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatcatc                                                                    <  2:3714900/93‑1 (MQ=255)

                    ggAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatcatc                                                                    <  2:3916259/93‑1 (MQ=255)

                     gAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAAtcatca                                                                   >  2:344422/1‑93 (MQ=255)

                     gAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAAtcatca                                                                   >  2:1080747/1‑93 (MQ=255)

                        aTAGTGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATcagca                                                                <  2:252771/93‑1 (MQ=255)

                        aTAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGGCGACAAAATAGGCGGCAATCATcagca                                                                >  1:2153713/1‑93 (MQ=255)

                        aTAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATcagca                                                                <  2:1410752/93‑1 (MQ=255)

                         tAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATcagcag                                                               >  1:903600/1‑93 (MQ=255)

                         tAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATcagaag                                                               >  1:2208732/1‑93 (MQ=255)

                          aGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAgg                                                              >  2:4467054/1‑93 (MQ=255)

                          aGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAgg                                                              >  1:3547688/1‑93 (MQ=255)

                           gggTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGt                                                             <  1:245214/93‑1 (MQ=255)

                            ggTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTa                                                            >  1:2283094/1‑93 (MQ=255)

                            ggTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTa                                                            <  1:1343222/93‑1 (MQ=255)

                             gTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTaa                                                           <  1:888820/93‑1 (MQ=255)

                              ttGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAg                                                          <  1:929832/93‑1 (MQ=255)

                              ttGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAg                                                          <  1:2511829/93‑1 (MQ=255)

                               tGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAgg                                                         >  1:2086296/1‑93 (MQ=255)

                                tccATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAggg                                                        <  1:1000437/92‑1 (MQ=255)

                                gCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAggg                                                        <  2:3069198/93‑1 (MQ=255)

                                gCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAggg                                                        <  2:399485/93‑1 (MQ=255)

                                gCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAggg                                                        <  1:1000436/93‑1 (MQ=255)

                                 ccATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAgggg                                                       >  2:2559792/1‑93 (MQ=255)

                                  cATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGa                                                      >  2:3308874/1‑93 (MQ=255)

                                  cATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGa                                                      <  1:951634/93‑1 (MQ=255)

                                  cATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGa                                                      >  2:1130283/1‑93 (MQ=255)

                                   aTCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAg                                                     >  2:2397139/1‑93 (MQ=255)

                                    tCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGt                                                    >  2:1950312/1‑93 (MQ=255)

                                    acAGCAAAACACCCAGACTGGTTGCGTAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGt                                                    <  2:3599631/92‑1 (MQ=255)

                                     cAGCAAAAGACCCAGACCGGGTGCGGAGGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGtt                                                   <  1:4706744/93‑1 (MQ=255)

                                     cAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGtt                                                   <  1:2270848/93‑1 (MQ=255)

                                      aGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGCTa                                                  >  2:85748/1‑93 (MQ=255)

                                       gCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGGCGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTaa                                                 <  1:4001979/93‑1 (MQ=255)

                                       gCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTaa                                                 >  2:3690326/1‑93 (MQ=255)

                                       gCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTaa                                                 >  1:371973/1‑93 (MQ=255)

                                         aaaaCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAata                                               >  1:1894076/1‑93 (MQ=255)

                                         aaaaCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAata                                               >  2:1015132/1‑93 (MQ=255)

                                         aaaaCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAata                                               <  2:566016/93‑1 (MQ=255)

                                          aaaCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAatat                                              >  2:58600/1‑93 (MQ=255)

                                           aaCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAAtata                                             <  1:1175369/93‑1 (MQ=255)

                                            acacCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGAGGCAATCATCAGCAGGTAAGGGGATTTAATATAc                                            <  2:2522409/93‑1 (MQ=255)

                                             cacCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACt                                           <  1:1494951/93‑1 (MQ=255)

                                             cacCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACt                                           <  1:740762/93‑1 (MQ=255)

                                               cccAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGc                                         >  2:3851662/1‑93 (MQ=255)

                                               cccAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGc                                         >  1:4423323/1‑93 (MQ=255)

                                                ccAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCa                                        <  1:4485104/93‑1 (MQ=255)

                                                 cAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCaa                                       <  2:616026/93‑1 (MQ=255)

                                                  aGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAAt                                      <  1:1574866/93‑1 (MQ=255)

                                                   gACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATg                                     >  2:4533140/1‑93 (MQ=255)

                                                     ccGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGt                                   <  1:4543539/93‑1 (MQ=255)

                                                        gTTGCGGAAGAGACGGCCAGAGANATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGtttt                                <  1:2189082/93‑1 (MQ=255)

                                                          tGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGa                              >  2:322323/1‑93 (MQ=255)

                                                           gCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGAc                             <  2:1941506/93‑1 (MQ=255)

                                                           gCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGAc                             <  2:1865242/93‑1 (MQ=255)

                                                              gAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGcc                          <  2:2696361/93‑1 (MQ=255)

                                                              gAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGcc                          <  2:3699128/93‑1 (MQ=255)

                                                              gAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGcc                          >  2:1431919/1‑93 (MQ=255)

                                                               aaGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCa                         <  2:1916806/93‑1 (MQ=255)

                                                               aaGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCa                         <  2:1916804/93‑1 (MQ=255)

                                                                agagACGGCCAGAGACATCAGACAGACGACACAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGATTTGACGCCAg                        >  1:1323385/1‑93 (MQ=255)

                                                                agagACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAg                        >  1:4530821/1‑93 (MQ=255)

                                                                  agaCGGCCAGAGACATCAGACAGTCGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCt                      <  2:2595970/93‑1 (MQ=255)

                                                                  agaCGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCt                      <  2:3093137/93‑1 (MQ=255)

                                                                   aaCGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCtt                     <  1:214881/92‑1 (MQ=255)

                                                                    aCGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTg                    <  1:4358727/93‑1 (MQ=255)

                                                                      ggCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGac                  <  2:722268/93‑1 (MQ=255)

                                                                      ggCAAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGac                  <  2:530818/93‑1 (MQ=255)

                                                                       gCCAGAGACATCAGACAGACGACAAACTAGGCGGCAATCGTCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGacc                 <  2:4210623/93‑1 (MQ=255)

                                                                       gCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGacc                 >  2:31215/1‑93 (MQ=255)

                                                                          agagACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGaccacc              >  1:3147153/1‑93 (MQ=255)

                                                                           gagaCATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGAccacca             >  2:2894570/1‑93 (MQ=255)

                                                                               cATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATc         <  1:2447080/93‑1 (MQ=255)

                                                                               cATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATc         >  1:1182359/1‑93 (MQ=255)

                                                                                aTCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCa        <  1:1096241/93‑1 (MQ=255)

                                                                                aTCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCa        >  1:4726767/1‑93 (MQ=255)

                                                                                aTCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCa        <  1:1180603/93‑1 (MQ=255)

                                                                                 tCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCAt       <  1:316563/93‑1 (MQ=255)

                                                                                  cagacagaCGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCAtt      <  1:1771293/93‑1 (MQ=255)

                                                                                   agacagacGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCATTc     >  2:2111319/1‑93 (MQ=255)

                                                                                      cAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCATTcgcg  >  2:1034089/1‑93 (MQ=255)

TACCAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGGCGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCATTCGCG  >  NC_000913/654723‑654901

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