BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,493	T→C	C124C (TGT→TGC)	yidX →	putative replicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,493	0	T	C	99.9%	2502.0 / ‑5.7	755	C124C (TGT→TGC)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (375/379); total (376/379)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

Only 100 of 768 total aligned reads displayed.

CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATT  >  NC_000913/3873403‑3873582

cGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCAc                                                                                         >  2:7248707/1‑93 (MQ=255)

cGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCAc                                                                                         >  2:4672201/1‑93 (MQ=255)

  gCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCa                                                                                       <  2:3145978/93‑1 (MQ=255)

   cacaATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAg                                                                                      <  1:690791/93‑1 (MQ=255)

   cacaATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAg                                                                                      <  1:1383946/93‑1 (MQ=255)

   cacaATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAg                                                                                      <  2:8761169/93‑1 (MQ=255)

    acaATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAgg                                                                                     >  2:4082946/1‑93 (MQ=255)

    acaATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAgg                                                                                     >  2:7668286/1‑93 (MQ=255)

         atCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACt                                                                                <  1:9133845/93‑1 (MQ=255)

         atCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACt                                                                                <  1:3710811/93‑1 (MQ=255)

          tCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTg                                                                               <  2:9342976/93‑1 (MQ=255)

            gCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGa                                                                             >  1:3312834/1‑93 (MQ=255)

            gCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGa                                                                             <  2:6758481/93‑1 (MQ=255)

             cTTCACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaa                                                                            <  2:5435661/93‑1 (MQ=255)

               ttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaaa                                                                          <  1:2728098/93‑1 (MQ=255)

               ttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaaa                                                                          >  1:4045042/1‑93 (MQ=255)

                  ctactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAAAACTGCACCAGGTCACTGAAAAACCt                                                                       >  2:7910527/1‑93 (MQ=255)

                   taccactGGAATAACGTAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCATGTCACTGGAAAACCTa                                                                      >  2:8920421/1‑93 (MQ=255)

                    actactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaa                                                                     >  2:7002003/1‑93 (MQ=255)

                     ctactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaa                                                                    >  2:10026769/1‑93 (MQ=255)

                     ctactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaa                                                                    <  2:2417122/93‑1 (MQ=255)

                      tactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaaa                                                                   >  1:840978/1‑93 (MQ=255)

                        ctGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAgg                                                                 >  1:6602427/1‑93 (MQ=255)

                         tGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                                <  1:1506161/93‑1 (MQ=255)

                          ggAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGca                                                               <  2:5922293/93‑1 (MQ=255)

                          ggAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGca                                                               >  1:4885846/1‑93 (MQ=255)

                           gAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGcac                                                              >  2:1720819/1‑93 (MQ=255)

                            aaTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGcaca                                                             >  1:7657051/1‑93 (MQ=255)

                             aTAACGAAGAAGAAAAGGCACTGGCAACGTATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                            <  2:1473076/93‑1 (MQ=255)

                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                            <  2:654572/93‑1 (MQ=255)

                              tAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAAt                                                           >  2:9118380/1‑93 (MQ=255)

                              tAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAAt                                                           >  2:1324499/1‑93 (MQ=255)

                              tAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAAt                                                           <  1:9224415/93‑1 (MQ=255)

                               aaCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATc                                                          <  1:5503328/93‑1 (MQ=255)

                                aCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATcc                                                         >  1:8878494/1‑93 (MQ=255)

                                aCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATcc                                                         <  2:339427/93‑1 (MQ=255)

                                aCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATcc                                                         <  2:6024382/93‑1 (MQ=255)

                                    agaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATga                                                     >  1:6425778/1‑93 (MQ=255)

                                    agaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATga                                                     <  1:4426791/93‑1 (MQ=255)

                                     gaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgag                                                    <  1:4579539/93‑1 (MQ=255)

                                      aagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgaga                                                   >  1:7879819/1‑93 (MQ=255)

                                      aagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgaga                                                   <  1:123368/93‑1 (MQ=255)

                                      aagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgaga                                                   >  1:7172031/1‑93 (MQ=255)

                                        gaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaa                                                 <  1:583815/93‑1 (MQ=255)

                                         aaaaGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaaa                                                >  2:3909957/1‑93 (MQ=255)

                                         aaaaGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaaa                                                >  1:2519388/1‑93 (MQ=255)

                                         aaaaGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaaa                                                >  2:8505191/1‑93 (MQ=255)

                                           aaGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAAt                                              <  2:55873/93‑1 (MQ=255)

                                           aaGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAAt                                              >  1:6680783/1‑93 (MQ=255)

                                            aGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATa                                             >  2:4842787/1‑93 (MQ=255)

                                              gCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaa                                           <  2:6431060/93‑1 (MQ=255)

                                               cACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaaa                                          <  2:8057988/93‑1 (MQ=255)

                                               cACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaaa                                          >  2:1061124/1‑93 (MQ=255)

                                               cACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaaa                                          <  1:5256227/93‑1 (MQ=255)

                                                aCTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaaaa                                         >  1:5032299/1‑93 (MQ=255)

                                                aCTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaaaa                                         <  1:2993390/93‑1 (MQ=255)

                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                        >  2:4490021/1‑93 (MQ=255)

                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                        <  1:5447550/93‑1 (MQ=255)

                                                  tGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACa                                       <  1:4742818/93‑1 (MQ=255)

                                                     cAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATgg                                    <  1:2202309/93‑1 (MQ=255)

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                   >  1:747985/1‑93 (MQ=255)

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                   >  2:6297166/1‑93 (MQ=255)

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                   <  2:290278/93‑1 (MQ=255)

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                   <  1:5920799/93‑1 (MQ=255)

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                   <  1:1982981/93‑1 (MQ=255)

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                   <  2:9261195/93‑1 (MQ=255)

                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                 <  2:2078395/93‑1 (MQ=255)

                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                 <  2:9466435/93‑1 (MQ=255)

                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                 <  1:3919950/93‑1 (MQ=255)

                                                         gCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTa                                <  1:8047510/93‑1 (MQ=255)

                                                           aTTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACt                              >  2:4277840/1‑93 (MQ=255)

                                                           aTTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACt                              >  1:9812802/1‑93 (MQ=255)

                                                           aTTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACt                              >  1:2395850/1‑93 (MQ=255)

                                                            ttATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTc                             >  2:3703299/1‑93 (MQ=255)

                                                            ttATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTc                             >  1:4326562/1‑93 (MQ=255)

                                                            ttATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTc                             >  1:10011480/1‑93 (MQ=255)

                                                              aTGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCaa                           <  2:8334825/93‑1 (MQ=255)

                                                                ggTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAg                         <  1:9730196/93‑1 (MQ=255)

                                                                  ttttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAggg                       >  1:8277563/1‑92 (MQ=255)

                                                                  ttttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGt                       >  2:5262256/1‑93 (MQ=255)

                                                                  ttttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGt                       >  1:6123499/1‑93 (MQ=255)

                                                                     tGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGAt                    <  1:1090466/93‑1 (MQ=255)

                                                                     tGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGAt                    <  1:958262/93‑1 (MQ=255)

                                                                       ccTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATgg                  <  2:8145771/93‑1 (MQ=255)

                                                                         tgtgGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCg                >  1:5721878/1‑93 (MQ=255)

                                                                         tgtgGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCg                >  1:7044374/1‑93 (MQ=255)

                                                                         tgtgGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCg                <  1:9345885/93‑1 (MQ=255)

                                                                            ggtgTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTc             >  2:9893855/1‑93 (MQ=255)

                                                                            ggtgTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTc             <  1:1665848/93‑1 (MQ=255)

                                                                            ggtgTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTc             >  1:8692263/1‑93 (MQ=255)

                                                                              tgtTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTa           >  1:9006671/1‑93 (MQ=255)

                                                                                ttCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAcc         >  2:2610935/1‑93 (MQ=255)

                                                                                gtCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAcc         <  1:7338772/92‑1 (MQ=255)

                                                                                  cAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccat       <  1:3531163/93‑1 (MQ=255)

                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatc      <  2:1585354/93‑1 (MQ=255)

                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatc      <  2:3507291/93‑1 (MQ=255)

                                                                                      acTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatccat   >  2:1832894/1‑93 (MQ=255)

                                                                                       cTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAtt  <  2:568313/93‑1 (MQ=255)

                                                                                       cTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAtt  >  1:2630595/1‑93 (MQ=255)

                                                                                       cTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAtt  >  2:5097201/1‑93 (MQ=255)

CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATT  >  NC_000913/3873403‑3873582

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