BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	526,519	T→C	C1345C (TGT→TGC)	rhsD →	RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	526,519	0	T	C	99.2%	387.9 / ‑4.5	125	C1345C (TGT→TGC)	rhsD	RhsD
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (56/68); total (56/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

Only 100 of 126 total aligned reads displayed.

TGTCTGTAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGA  >  NC_000913/526429‑526608

tGTCTGTAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGccc                                                                                         <  1:2303338/93‑1 (MQ=255)

  tCTGTAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCca                                                                                       <  2:520141/93‑1 (MQ=255)

  tCTGTAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCca                                                                                       >  2:3510482/1‑93 (MQ=255)

    tGTAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCcaca                                                                                     <  1:1155462/93‑1 (MQ=255)

    tGTAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCcaca                                                                                     <  2:1113375/93‑1 (MQ=255)

     gTAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAg                                                                                    >  2:3038100/1‑93 (MQ=255)

     gTAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAg                                                                                    >  2:3566195/1‑93 (MQ=255)

     gTAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAg                                                                                    >  1:2198551/1‑93 (MQ=255)

      tAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAgg                                                                                   <  2:3922234/93‑1 (MQ=255)

      tAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAgg                                                                                   >  1:1561517/1‑93 (MQ=255)

         aGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTAc                                                                                <  1:3456440/93‑1 (MQ=255)

          gATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACa                                                                               <  1:4185183/93‑1 (MQ=255)

          gATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACa                                                                               <  2:4309567/93‑1 (MQ=255)

          gATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTAAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACa                                                                               <  2:913290/93‑1 (MQ=255)

            tACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGa                                                                             <  2:945236/93‑1 (MQ=255)

             aCTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGAt                                                                            >  2:3684926/1‑93 (MQ=255)

              cTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTATGAAAGGAGTCATATACTCTGAAAAGAAAATGCCCCACAGGCACAGATa                                                                           >  2:1679106/1‑93 (MQ=255)

              cTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATa                                                                           <  2:2134523/93‑1 (MQ=255)

               tAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATag                                                                          >  2:4659610/1‑93 (MQ=255)

               tAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATag                                                                          <  1:2707848/93‑1 (MQ=255)

               tAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATag                                                                          >  2:36676/1‑93 (MQ=255)

               tAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATag                                                                          >  1:843158/1‑93 (MQ=255)

                aaTGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATaga                                                                         >  2:2741468/1‑93 (MQ=255)

                aaTGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATaga                                                                         <  2:1955068/93‑1 (MQ=255)

                  tGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGt                                                                       <  1:2216944/93‑1 (MQ=255)

                   ggAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGtt                                                                      <  2:1197436/93‑1 (MQ=255)

                   ggAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGtt                                                                      >  2:2359486/1‑93 (MQ=255)

                   ggAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGtt                                                                      <  1:394523/93‑1 (MQ=255)

                   ggAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGtt                                                                      >  2:1264141/1‑93 (MQ=255)

                    gAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTtg                                                                     <  1:3308267/93‑1 (MQ=255)

                     aaaaTATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTtgc                                                                    <  2:2020933/93‑1 (MQ=255)

                       aaTATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTtgctg                                                                  >  2:4309350/1‑93 (MQ=255)

                          atTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCtt                                                               >  2:4472443/1‑93 (MQ=255)

                          atTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCtt                                                               >  1:807369/1‑93 (MQ=255)

                           ttttGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTa                                                              <  2:981279/93‑1 (MQ=255)

                             ttGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATc                                                            <  1:361418/93‑1 (MQ=255)

                              tGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCa                                                           <  1:1401488/93‑1 (MQ=255)

                               gCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCAt                                                          <  1:3627430/93‑1 (MQ=255)

                                cATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATa                                                         <  2:3543279/93‑1 (MQ=255)

                                 aTCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATAc                                                        >  2:1006241/1‑93 (MQ=255)

                                 aTCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCCGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATAc                                                        <  2:545471/93‑1 (MQ=255)

                                  tCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACt                                                       >  1:106997/1‑93 (MQ=255)

                                  tCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACt                                                       >  1:180362/1‑93 (MQ=255)

                                  tCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACt                                                       >  1:3013376/1‑93 (MQ=255)

                                    tAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCa                                                     >  1:536934/1‑93 (MQ=255)

                                    tAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCa                                                     >  1:1128791/1‑93 (MQ=255)

                                    tAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCa                                                     <  1:2822424/93‑1 (MQ=255)

                                       ggCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACgg                                                  <  2:1616094/93‑1 (MQ=255)

                                       ggCAGAAACTGATAATGTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACgg                                                  <  2:4641342/93‑1 (MQ=255)

                                        gCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGt                                                 >  2:617562/1‑93 (MQ=255)

                                          aGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGc                                               >  2:2432472/1‑93 (MQ=255)

                                            aaaCTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAg                                             >  2:817813/1‑93 (MQ=255)

                                             aaCTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGa                                            >  2:4016395/1‑93 (MQ=255)

                                             aaCTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGa                                            >  1:1533949/1‑93 (MQ=255)

                                               cTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATa                                          <  1:1412887/93‑1 (MQ=255)

                                                tGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAg                                         <  2:1185053/93‑1 (MQ=255)

                                                tGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAg                                         >  1:233339/1‑93 (MQ=255)

                                                 gATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGt                                        <  1:3655492/93‑1 (MQ=255)

                                                  aTAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTc                                       <  2:4234404/93‑1 (MQ=255)

                                                  aTAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGGTGCTGCTTATCATACTCACGGTGCAGATAGTc                                       >  1:2858420/1‑93 (MQ=255)

                                                    aaTTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCAt                                     <  1:3948820/93‑1 (MQ=255)

                                                    aaTTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCAt                                     <  2:2675826/93‑1 (MQ=255)

                                                     aTTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATg                                    <  1:3995482/93‑1 (MQ=255)

                                                     aTTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATg                                    >  1:3133762/1‑93 (MQ=255)

                                                       ttAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGc                                  >  1:3225691/1‑93 (MQ=255)

                                                         aaGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATAATACTCACGGTGCAGATAGTCATGGCGa                                <  1:3270263/93‑1 (MQ=255)

                                                          aGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGAt                               >  2:2243579/1‑93 (MQ=255)

                                                          aGAAAGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGAt                               <  1:928045/93‑1 (MQ=255)

                                                             aaGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTAt                            <  1:280651/93‑1 (MQ=255)

                                                              aGGAGTCATATCCTCTGAATAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATg                           >  2:4581335/1‑93 (MQ=255)

                                                              aGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATg                           >  2:3773361/1‑93 (MQ=255)

                                                              aGGAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGCACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATg                           >  1:4160857/1‑93 (MQ=255)

                                                               ggAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGt                          >  1:1471506/1‑93 (MQ=255)

                                                               ggAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGt                          <  2:556614/93‑1 (MQ=255)

                                                                gAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGtt                         >  1:4006839/1‑93 (MQ=255)

                                                                gAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGtt                         <  1:595634/93‑1 (MQ=255)

                                                                gAGTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGtt                         <  2:1018215/93‑1 (MQ=255)

                                                                  gTCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTtga                       <  2:2211250/93‑1 (MQ=255)

                                                                   tCATATCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTtgat                      <  2:848887/93‑1 (MQ=255)

                                                                      tatCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGaa                   >  2:995100/1‑93 (MQ=255)

                                                                       atCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAAt                  >  2:344834/1‑93 (MQ=255)

                                                                       atCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAAt                  <  1:1291884/93‑1 (MQ=255)

                                                                        tCCTCTGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAAtt                 <  2:188781/93‑1 (MQ=255)

                                                                          ctctGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAAtttt               >  2:3945103/1‑93 (MQ=255)

                                                                          ctctGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAAtttt               >  2:189991/1‑93 (MQ=255)

                                                                           tctGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAAttttt              <  2:4091849/93‑1 (MQ=255)

                                                                            ctGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAAtttttt             >  1:1083790/1‑93 (MQ=255)

                                                                            ctGAAAAGAAAATGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAAtttttt             <  1:249855/93‑1 (MQ=255)

                                                                              gaaaagaaaaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTc           >  1:3699169/1‑93 (MQ=255)

                                                                               aaaagaaaaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCa          <  2:301352/93‑1 (MQ=255)

                                                                               aaaagaaaaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCa          <  1:3564522/93‑1 (MQ=255)

                                                                                aaagaaaaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCaa         >  1:1476429/1‑93 (MQ=255)

                                                                                aaagaaaaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCaa         <  1:2069559/93‑1 (MQ=255)

                                                                                  agaaaaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGt       <  2:525055/93‑1 (MQ=255)

                                                                                  agaaaaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGt       <  1:4508664/93‑1 (MQ=255)

                                                                                   gaaaaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTa      <  1:4184853/93‑1 (MQ=255)

                                                                                     aaaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGc    <  1:2603976/93‑1 (MQ=255)

                                                                                     aaaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGACAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGc    <  1:1353794/93‑1 (MQ=255)

                                                                                      aaTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAACTAGCg   >  2:409843/1‑93 (MQ=255)

                                                                                       aTGCCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGa  >  2:2079524/1‑93 (MQ=255)

TGTCTGTAAAGATACTAATGGAAAATATTTTGCATCTAAGGCAGAAACTGATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGA  >  NC_000913/526429‑526608

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