BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,071,901	G→A	A197T (GCC→ACC)	yihV →	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,071,901	0	G	A	97.4%	514.3 / ‑0.6	152	A197T (GCC→ACC)	yihV	putative kinase
Reads supporting (aligned to +/- strand): ref base G (2/2); new base A (77/71); total (79/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 155 total aligned reads displayed.

CCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGGCCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAATGG  >  NC_000913/4071811‑4071989

ccgcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGAcc                                                                                        <  1:4617415/93‑1 (MQ=255)

 cgcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCa                                                                                       >  2:4134902/1‑93 (MQ=255)

  gcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAg                                                                                      <  1:1302250/93‑1 (MQ=255)

  gcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAg                                                                                      <  2:52836/93‑1 (MQ=255)

    aGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTg                                                                                    >  1:1147108/1‑93 (MQ=255)

    aGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTg                                                                                    >  2:2736926/1‑93 (MQ=255)

     ggATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTgc                                                                                   >  2:4384203/1‑93 (MQ=255)

      gATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTgcg                                                                                  >  2:508711/1‑93 (MQ=255)

      gATATCAGCGGGCTGGTGTCATTAAGCGATCCCGCGGACTTTTAAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTgcg                                                                                  >  2:3761518/1‑93 (MQ=255)

        tatCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCt                                                                                >  1:3814880/1‑93 (MQ=255)

          tCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTaa                                                                              >  1:3621855/1‑93 (MQ=255)

          tCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTaa                                                                              >  1:2725060/1‑93 (MQ=255)

             gTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAAc                                                                           >  2:1793338/1‑93 (MQ=255)

              tGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACa                                                                          <  1:617791/93‑1 (MQ=255)

               gagaTGGTGGAATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAg                                                                         <  1:115544/93‑1 (MQ=255)

               gAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGGCCAGTGCGCTAAAACAg                                                                         >  2:1997961/1‑93 (MQ=255)

               gAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAg                                                                         <  2:2300990/93‑1 (MQ=255)

                 gCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTATGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGc                                                                       <  1:333919/93‑1 (MQ=255)

                  cTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGca                                                                      <  1:1792848/93‑1 (MQ=255)

                   cggtggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGcac                                                                     <  2:3697006/92‑1 (MQ=255)

                    ggtggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGcaca                                                                    >  1:2514471/1‑93 (MQ=255)

                    gggggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGcaca                                                                    >  2:1065434/4‑93 (MQ=255)

                      tggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACaaa                                                                  <  1:4160419/93‑1 (MQ=255)

                      tggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACaaa                                                                  >  2:3121668/1‑93 (MQ=255)

                        gCATTAAGCGCTCCTGCGACCTCTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAATCAGGCACAAACg                                                                <  2:4165054/93‑1 (MQ=255)

                        gCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACg                                                                >  1:4103966/1‑93 (MQ=255)

                          aTTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCt                                                              >  2:3294584/1‑93 (MQ=255)

                          aTTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCt                                                              >  1:839112/1‑93 (MQ=255)

                           ttAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTc                                                             >  1:3091749/1‑93 (MQ=255)

                           ttAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGAACAAACGCTc                                                             >  2:3646178/1‑93 (MQ=255)

                            tAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTca                                                            <  2:3452704/93‑1 (MQ=255)

                                cGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACaaa                                                        >  1:3559110/1‑93 (MQ=255)

                                cGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACaaa                                                        <  2:2677204/93‑1 (MQ=255)

                                 gATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCCGTGCGCTAAAACAGGCACAAACGCTCACAAAt                                                       >  1:3896108/1‑93 (MQ=255)

                                 gATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAAt                                                       >  1:3896107/1‑93 (MQ=255)

                                 gATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAAt                                                       >  1:2485448/1‑93 (MQ=255)

                                  aTCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATg                                                      <  1:4520083/93‑1 (MQ=255)

                                    cACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTAACAAATggc                                                    >  1:46852/1‑92 (MQ=255)

                                      cgcgGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCa                                                  >  2:235229/1‑93 (MQ=255)

                                      cgcgGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCa                                                  <  1:1528188/93‑1 (MQ=255)

                                       gcgGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCAt                                                 <  2:1555163/93‑1 (MQ=255)

                                       gcgGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCAt                                                 <  1:3094115/93‑1 (MQ=255)

                                        cgGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATg                                                <  1:2422696/93‑1 (MQ=255)

                                         ggCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGt                                               <  2:296420/93‑1 (MQ=255)

                                         ggCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGt                                               >  1:4201375/1‑93 (MQ=255)

                                         ggCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGt                                               >  2:4417904/1‑93 (MQ=255)

                                          gCCTTTTCAGCACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTc                                              >  2:2950739/1‑93 (MQ=255)

                                          gCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTc                                              >  1:3847294/1‑93 (MQ=255)

                                          gCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTc                                              >  2:1865925/1‑93 (MQ=255)

                                           ccTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCt                                             >  2:1984907/1‑93 (MQ=255)

                                            cTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTa                                            <  2:4503175/93‑1 (MQ=255)

                                            cTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTa                                            <  2:1897356/93‑1 (MQ=255)

                                             ttttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAt                                           <  2:2565509/93‑1 (MQ=255)

                                             ttttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAt                                           >  2:221749/1‑93 (MQ=255)

                                              tttCGGCACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCCGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAtg                                          <  2:4220015/93‑1 (MQ=255)

                                              tttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAtg                                          <  1:4428995/93‑1 (MQ=255)

                                               ttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAtgt                                         <  1:483059/93‑1 (MQ=255)

                                                 cAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGGCCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGa                                       <  1:2004675/93‑1 (MQ=255)

                                                  aGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAc                                      <  1:4456490/93‑1 (MQ=255)

                                                  aGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAc                                      <  2:3362673/93‑1 (MQ=255)

                                                  aGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAc                                      >  2:352439/1‑93 (MQ=255)

                                                   gAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAcc                                     >  1:2512938/1‑93 (MQ=255)

                                                    aaCCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAccc                                    >  1:905913/1‑93 (MQ=255)

                                                     aCCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCa                                   <  1:1780064/93‑1 (MQ=255)

                                                      ccGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAg                                  <  2:3370433/93‑1 (MQ=255)

                                                      ccGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAg                                  >  1:2343675/1‑93 (MQ=255)

                                                       cGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAgg                                 <  2:2042861/93‑1 (MQ=255)

                                                       cGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAgg                                 <  2:1286902/93‑1 (MQ=255)

                                                       cGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAgg                                 >  1:2039802/1‑93 (MQ=255)

                                                        tggTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCTCAATTGGTCATGTCTATGTGACCCAggg                                >  1:4222888/2‑93 (MQ=255)

                                                        gggTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAggg                                >  1:98826/1‑93 (MQ=255)

                                                         ggTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAggga                               >  1:4202480/1‑92 (MQ=255)

                                                          gTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTa                              <  2:1704964/93‑1 (MQ=255)

                                                             tGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAgcg                           >  2:1483448/1‑93 (MQ=255)

                                                             tGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCAAAAATGGTCATGTCTATGTGACCCAGGGTAgcg                           >  2:3334800/1‑93 (MQ=255)

                                                              ggCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAgcgc                          >  2:4343458/1‑93 (MQ=255)

                                                              ggCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAgcgc                          <  2:603006/93‑1 (MQ=255)

                                                               gcgcgcTTAACGAGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCa                         <  2:3078980/93‑1 (MQ=255)

                                                                 gcgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAgg                       >  1:251584/1‑93 (MQ=255)

                                                                 gcgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAgg                       >  2:1628101/1‑93 (MQ=255)

                                                                  cgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGc                      <  1:1589605/93‑1 (MQ=255)

                                                                   gcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCt                     >  2:4593975/1‑93 (MQ=255)

                                                                     ttAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGc                   <  2:217342/93‑1 (MQ=255)

                                                                         cGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGAct               <  2:2056272/93‑1 (MQ=255)

                                                                           ggCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActgg             >  1:802176/1‑93 (MQ=255)

                                                                           ggCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActgg             >  1:802177/1‑93 (MQ=255)

                                                                             cGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggct           <  1:2198234/93‑1 (MQ=255)

                                                                              gTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggctg          >  2:1112494/1‑93 (MQ=255)

                                                                              gTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggctg          >  1:2601871/1‑93 (MQ=255)

                                                                              gTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggctg          >  2:189859/1‑93 (MQ=255)

                                                                              gTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggctg          >  2:2299246/1‑93 (MQ=255)

                                                                                 aaaGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaa       >  1:40213/1‑93 (MQ=255)

                                                                                  aaGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaaa      <  1:4261652/93‑1 (MQ=255)

                                                                                  aaGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaaa      >  1:3351561/1‑93 (MQ=255)

                                                                                   agagATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaaaa     >  2:1721225/1‑93 (MQ=255)

                                                                                    gagaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAt    >  1:2933010/1‑93 (MQ=255)

                                                                                     agaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAtg   <  2:2620568/93‑1 (MQ=255)

                                                                                      gaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAtgg  <  2:1766349/93‑1 (MQ=255)

                                                                                      gaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAtgg  >  2:4663266/1‑93 (MQ=255)

                                                                                      gaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAtgg  >  1:1672931/1‑93 (MQ=255)

CCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGGCCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAATGG  >  NC_000913/4071811‑4071989

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