BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,897,800	A→G	Y166H (TAT→CAT)	b2774 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,897,800	0	A	G	95.8%	878.0 / 18.8	259	Y166H (TAT→CAT)	b2774	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (5/6); new base G (138/110); total (143/116)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

Only 100 of 265 total aligned reads displayed.

GGGCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAG  >  NC_000913/2897710‑2897889

gggCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAAtgtg                                                                                         <  1:4809705/93‑1 (MQ=255)

 ggCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGc                                                                                        <  2:636698/93‑1 (MQ=255)

  gCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGTCAAGAGCATGTTTAGTGGCAGAATGTGCa                                                                                       >  2:956615/1‑93 (MQ=255)

   cGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAg                                                                                      >  1:6119675/1‑93 (MQ=255)

    gATACCATTTACCTGAATATTATATTGACGTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAgg                                                                                     >  1:3884851/1‑93 (MQ=255)

    gATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAgg                                                                                     <  2:3766258/93‑1 (MQ=255)

    gATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAgg                                                                                     >  1:4851036/1‑93 (MQ=255)

    gATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAgg                                                                                     <  1:4051715/93‑1 (MQ=255)

    gATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAgg                                                                                     <  2:2235749/93‑1 (MQ=255)

     aTACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGt                                                                                    <  2:5071038/93‑1 (MQ=255)

      tACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTg                                                                                   >  2:4334581/1‑93 (MQ=255)

      tACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTg                                                                                   >  2:1002848/1‑93 (MQ=255)

       aCCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGa                                                                                  >  2:3887640/1‑93 (MQ=255)

        ccATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAACAGCAGGATTAAGGGCAGAATATGCAGGGGAc                                                                                 >  1:5893315/1‑93 (MQ=255)

         cATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGAcc                                                                                >  2:2555856/1‑93 (MQ=255)

          aTTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCa                                                                               >  2:1364245/1‑93 (MQ=255)

           tttACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCAt                                                                              >  1:839377/1‑93 (MQ=255)

           tttACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCAt                                                                              >  1:3559835/1‑93 (MQ=255)

              aCCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGt                                                                           <  2:1057424/93‑1 (MQ=255)

              aCCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGt                                                                           >  2:560291/1‑93 (MQ=255)

              aCCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGt                                                                           >  2:3561548/1‑93 (MQ=255)

              aCCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTAACCATTGt                                                                           >  2:2010563/1‑93 (MQ=255)

                cTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTcc                                                                         >  1:5887164/1‑93 (MQ=255)

                cTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTcc                                                                         <  1:546002/93‑1 (MQ=255)

                 tGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCa                                                                        >  1:1309787/1‑93 (MQ=255)

                  gAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCAc                                                                       <  1:1184180/93‑1 (MQ=255)

                  gAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCAc                                                                       <  2:3321699/93‑1 (MQ=255)

                   aaTATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCAcc                                                                      >  1:719613/1‑93 (MQ=255)

                   aaTATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCAcc                                                                      >  2:2194061/1‑93 (MQ=255)

                   aaTATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCAcc                                                                      >  2:5564498/1‑93 (MQ=255)

                    atattatattGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCt                                                                     <  1:4361131/93‑1 (MQ=255)

                    atattatattGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCt                                                                     >  1:2371819/1‑93 (MQ=255)

                     tattatattGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCta                                                                    >  1:4208895/1‑93 (MQ=255)

                      attatattGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCtaa                                                                   <  1:3669263/93‑1 (MQ=255)

                      attatattGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCtaa                                                                   >  2:5101298/1‑93 (MQ=255)

                          tattGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACgtaagtaa                                                               <  1:448159/93‑1 (MQ=255)

                          tattGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCtaagtaa                                                               >  2:1886948/1‑93 (MQ=255)

                             tGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAagag                                                            <  2:610499/93‑1 (MQ=255)

                             tGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAagag                                                            <  1:5299429/93‑1 (MQ=255)

                             tGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAagag                                                            <  1:3021129/93‑1 (MQ=255)

                              gACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAgaga                                                           <  1:4459262/93‑1 (MQ=255)

                               aCCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGaa                                                          <  1:2199419/93‑1 (MQ=255)

                               aCCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGaa                                                          >  2:4385068/1‑93 (MQ=255)

                                ccTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAAc                                                         <  1:5373311/93‑1 (MQ=255)

                                    gTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATg                                                     <  1:2642706/93‑1 (MQ=255)

                                     ttCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGa                                                    >  1:3342836/1‑93 (MQ=255)

                                     ttCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGa                                                    >  2:2991302/1‑93 (MQ=255)

                                     ttCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGa                                                    <  2:3141671/93‑1 (MQ=255)

                                     ttCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGa                                                    <  2:350911/93‑1 (MQ=255)

                                     ttCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGa                                                    <  1:1455612/93‑1 (MQ=255)

                                      tcatcaCAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGaa                                                   >  1:3787602/1‑93 (MQ=255)

                                      tcatcaCAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGaa                                                   <  2:1295994/93‑1 (MQ=255)

                                      tcatcaCAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGaa                                                   >  1:4878160/1‑93 (MQ=255)

                                        atcaCAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACa                                                 <  1:1005241/93‑1 (MQ=255)

                                        atcaCAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACa                                                 <  1:3950510/93‑1 (MQ=255)

                                        atcaCAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACa                                                 >  1:2291794/1‑93 (MQ=255)

                                          cacaATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGa                                               >  1:1536957/1‑93 (MQ=255)

                                           acaATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGat                                              <  2:3629647/93‑1 (MQ=255)

                                            caATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGata                                             <  1:1069305/93‑1 (MQ=255)

                                            caATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGata                                             <  1:1532888/93‑1 (MQ=255)

                                             aaTAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGatat                                            <  1:4579668/93‑1 (MQ=255)

                                              aTAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATAtt                                           >  2:5326051/1‑93 (MQ=255)

                                              aTAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATAtt                                           <  2:371641/93‑1 (MQ=255)

                                                aaGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATtaa                                         <  1:674850/93‑1 (MQ=255)

                                                aaGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATtaa                                         <  2:2107965/93‑1 (MQ=255)

                                                 aGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATtaat                                        <  1:1909956/93‑1 (MQ=255)

                                                  gCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATtaata                                       <  1:3604521/93‑1 (MQ=255)

                                                  gCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATtaata                                       <  1:2408296/93‑1 (MQ=255)

                                                   cTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATtaataa                                      <  2:4583796/93‑1 (MQ=255)

                                                   cTTTGGTGAACCCGGCAAGAGCACGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATtaataa                                      <  1:4473566/93‑1 (MQ=255)

                                                    tttGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTaataat                                     >  2:5214390/1‑93 (MQ=255)

                                                    tttGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTaataat                                     <  2:5241500/93‑1 (MQ=255)

                                                    tttGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTaataat                                     <  2:5632598/93‑1 (MQ=255)

                                                     ttGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAAtt                                    >  2:4971278/1‑93 (MQ=255)

                                                          gAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGcc                               <  2:1538755/93‑1 (MQ=255)

                                                          gAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGcc                               <  2:4672427/93‑1 (MQ=255)

                                                          gAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGcc                               >  2:417830/1‑93 (MQ=255)

                                                          gAACCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGcc                               <  1:5714929/93‑1 (MQ=255)

                                                            aCCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCgt                             >  2:593687/1‑92 (MQ=255)

                                                            aCCCGGCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGc                             <  2:4513695/93‑1 (MQ=255)

                                                                ggCAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCtttt                         >  2:1112863/1‑93 (MQ=255)

                                                                 gCAAGAGCATGGTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTc                        >  2:2423333/1‑93 (MQ=255)

                                                                  cAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCt                       <  1:3852457/93‑1 (MQ=255)

                                                                  cAAGAGCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCt                       >  1:1971317/1‑93 (MQ=255)

                                                                    agagCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGc                     >  1:1804384/1‑93 (MQ=255)

                                                                    agagCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGc                     <  2:341959/93‑1 (MQ=255)

                                                                    agagCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGc                     >  1:1690607/1‑93 (MQ=255)

                                                                      agCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCGCCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCgg                   <  1:4729123/93‑1 (MQ=255)

                                                                      agCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCgg                   >  1:2957877/1‑93 (MQ=255)

                                                                       gCATGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCggg                  <  1:559827/93‑1 (MQ=255)

                                                                         aTGTTTAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGAt                <  2:406779/93‑1 (MQ=255)

                                                                            tttAGTGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCAt             >  1:1673942/1‑93 (MQ=255)

                                                                            tttAGTGGCAAAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCAt             <  1:5161115/93‑1 (MQ=255)

                                                                                 tGGCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCCGCGGGATCATAAttt        >  1:4983417/1‑93 (MQ=255)

                                                                                  ggCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAAtttt       >  2:807358/1‑93 (MQ=255)

                                                                                   gCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAAttttt      <  2:2709624/93‑1 (MQ=255)

                                                                                   gCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAAttttt      <  2:485082/93‑1 (MQ=255)

                                                                                   gCAGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAAttttt      >  2:1758866/1‑93 (MQ=255)

                                                                                     aGAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGc    <  2:5947895/93‑1 (MQ=255)

                                                                                      gAATGTGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCa   >  2:3757558/1‑93 (MQ=255)

GGGCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAG  >  NC_000913/2897710‑2897889

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