BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,570,594	G→A	A878V (GCG→GTG)	pqqL ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,570,594	0	G	A	99.5%	1358.3 / ‑5.9	375	A878V (GCG→GTG)	pqqL	putative peptidase
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (210/162); total (211/164)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 379 total aligned reads displayed.

TTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTA  >  NC_000913/1570504‑1570683

ttCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACt                                                                                         >  2:7734635/1‑93 (MQ=255)

 tctcTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTa                                                                                        >  2:8271026/1‑93 (MQ=255)

 tctcTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTa                                                                                        >  2:7144617/1‑93 (MQ=255)

  ctctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTaa                                                                                       <  2:6769136/93‑1 (MQ=255)

  ctctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTaa                                                                                       <  2:1986172/93‑1 (MQ=255)

   tctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAt                                                                                      >  1:1303377/1‑93 (MQ=255)

    ctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAtt                                                                                     <  2:653662/93‑1 (MQ=255)

      aCCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAttgt                                                                                   <  2:6933627/93‑1 (MQ=255)

       ccGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAttgtt                                                                                  <  2:9067757/93‑1 (MQ=255)

       ccGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAttgtt                                                                                  >  1:7779658/1‑93 (MQ=255)

        cGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAttgttg                                                                                 <  1:5033616/93‑1 (MQ=255)

        cGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAttgttg                                                                                 >  2:1075171/1‑93 (MQ=255)

         gTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGa                                                                                >  1:1072371/1‑93 (MQ=255)

         gTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGa                                                                                <  1:6726956/93‑1 (MQ=255)

          tCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGaa                                                                               >  2:3281456/1‑93 (MQ=255)

          tCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGaa                                                                               >  2:5699722/1‑93 (MQ=255)

          tCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGaa                                                                               >  2:4717718/1‑93 (MQ=255)

           cATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAac                                                                              >  1:4890716/1‑93 (MQ=255)

            aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacg                                                                             >  2:4057294/1‑93 (MQ=255)

            aTTTGTTTCAACAATTGATCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCCCTAATTGTTGaaag                                                                             >  1:7447439/1‑91 (MQ=255)

              ttGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACACTAGAGTTCACTAATTGTTGAacgct                                                                           >  1:8686757/1‑93 (MQ=255)

              ttGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgct                                                                           <  1:1149459/93‑1 (MQ=255)

               tGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgcta                                                                          >  1:761160/1‑93 (MQ=255)

                  ttCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgctacgc                                                                       <  2:8728595/93‑1 (MQ=255)

                   tCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgctacgct                                                                      >  2:4344220/1‑93 (MQ=255)

                   tCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgctacgct                                                                      >  1:678654/1‑93 (MQ=255)

                   tCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgctacgct                                                                      >  2:6679169/1‑93 (MQ=255)

                    caacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCAATAATTGTTGAACGCTACGCTg                                                                     >  1:3011662/1‑93 (MQ=255)

                     aacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGt                                                                    <  2:5046417/93‑1 (MQ=255)

                     aacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGt                                                                    >  1:9775018/1‑93 (MQ=255)

                     aacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGt                                                                    >  2:8444164/1‑93 (MQ=255)

                      acaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGtt                                                                   >  1:9332124/1‑93 (MQ=255)

                       caaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTg                                                                  <  2:6492530/93‑1 (MQ=255)

                          ttGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGat                                                               >  2:3740216/1‑93 (MQ=255)

                           tGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGata                                                              <  1:3177180/93‑1 (MQ=255)

                           tGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGata                                                              <  2:2229435/93‑1 (MQ=255)

                            gCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGatat                                                             <  1:417496/93‑1 (MQ=255)

                             cTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATc                                                            <  1:6144530/93‑1 (MQ=255)

                              tCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCg                                                           >  2:586936/1‑93 (MQ=255)

                               ccTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGa                                                          <  2:4614216/93‑1 (MQ=255)

                               ccTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGa                                                          <  2:1555960/93‑1 (MQ=255)

                               ccTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGa                                                          >  2:5251558/1‑93 (MQ=255)

                               ccTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGa                                                          >  1:3381837/1‑93 (MQ=255)

                               ccTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGa                                                          <  1:7954379/93‑1 (MQ=255)

                                cTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAg                                                         >  2:9493085/1‑93 (MQ=255)

                                cTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAg                                                         >  2:7666251/1‑93 (MQ=255)

                                 tGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAgg                                                        <  2:3954805/93‑1 (MQ=255)

                                 tGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAgg                                                        >  2:8541169/1‑93 (MQ=255)

                                  gCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGc                                                       >  1:3573628/1‑93 (MQ=255)

                                  gCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGc                                                       <  2:2964409/93‑1 (MQ=255)

                                   cTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCt                                                      >  2:9685351/1‑93 (MQ=255)

                                   cTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCt                                                      >  1:621361/1‑93 (MQ=255)

                                    tCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTg                                                     >  2:6125467/1‑93 (MQ=255)

                                    tCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTg                                                     >  1:8290067/1‑93 (MQ=255)

                                    tCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTg                                                     <  2:4278251/93‑1 (MQ=255)

                                      aGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTgcg                                                   >  1:5891060/1‑93 (MQ=255)

                                       gTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTgcgc                                                  >  2:7433363/1‑93 (MQ=255)

                                         ccATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTg                                                >  1:7016215/1‑93 (MQ=255)

                                            tgctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAAc                                             >  1:2753620/1‑93 (MQ=255)

                                            tgctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAAc                                             >  1:1498019/1‑93 (MQ=255)

                                              ctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGt                                           <  1:8980528/93‑1 (MQ=255)

                                               tgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACggt                                          >  1:7144959/1‑91 (MQ=255)

                                               tgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGtt                                          >  1:5378855/1‑93 (MQ=255)

                                                gcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGctt                                         >  2:9439967/1‑93 (MQ=255)

                                                    gATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAAGGTTTTGCt                                     >  1:2154844/1‑93 (MQ=255)

                                                     aTCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTg                                    <  2:3227989/93‑1 (MQ=255)

                                                      tCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTgg                                   >  2:1684476/1‑93 (MQ=255)

                                                         tCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTat                                >  1:9470658/1‑93 (MQ=255)

                                                            tatTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattca                             >  2:1458085/1‑93 (MQ=255)

                                                             atTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcat                            >  1:4741821/1‑93 (MQ=255)

                                                               tGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcattc                          <  2:4908485/93‑1 (MQ=255)

                                                                gAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGTTGAACGTTTTGCTGGTattcattcg                         >  2:7959679/1‑92 (MQ=255)

                                                                 aaTAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAg                        >  1:6743911/1‑93 (MQ=255)

                                                                   tAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGtt                      <  2:1227138/93‑1 (MQ=255)

                                                                   tAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGtt                      >  1:3958042/1‑93 (MQ=255)

                                                                     aGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCt                    >  1:1656656/1‑93 (MQ=255)

                                                                     aGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCt                    <  2:4867869/93‑1 (MQ=255)

                                                                     aGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCt                    >  1:4415522/1‑93 (MQ=255)

                                                                     aGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCt                    <  2:2726335/93‑1 (MQ=255)

                                                                     aGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCt                    <  2:3029223/93‑1 (MQ=255)

                                                                      gACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCtt                   >  2:2338361/1‑93 (MQ=255)

                                                                      gACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCtt                   >  2:3465166/1‑93 (MQ=255)

                                                                       aCTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTg                  <  1:3725996/93‑1 (MQ=255)

                                                                        cTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGc                 <  1:1818810/93‑1 (MQ=255)

                                                                         tATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCt                >  1:3849612/1‑93 (MQ=255)

                                                                         tATTTACAATAGTGTTCACTAAGTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACATCTTGCTGTTGTTCACTCAGCTCTTGcg                >  1:8871680/1‑92 (MQ=255)

                                                                           tttACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCa              >  1:8049526/1‑93 (MQ=255)

                                                                           tttACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCa              >  2:9880845/1‑93 (MQ=255)

                                                                            ttACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCAc             >  2:2759298/1‑93 (MQ=255)

                                                                             tACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTCTCTCACt            >  1:6641756/1‑93 (MQ=255)

                                                                                aaTAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGAt         <  1:4067616/93‑1 (MQ=255)

                                                                                aaTAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGAt         <  2:2107953/93‑1 (MQ=255)

                                                                                aaTAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGAt         >  2:9219670/1‑93 (MQ=255)

                                                                                aaTAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGAt         <  1:9880815/93‑1 (MQ=255)

                                                                                  tAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATcc       >  1:24485/1‑93 (MQ=255)

                                                                                    gtgtTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCt     <  2:9504982/93‑1 (MQ=255)

                                                                                    gtgtTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCt     <  2:4558468/93‑1 (MQ=255)

                                                                                      gtTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCttt   <  1:9615800/93‑1 (MQ=255)

                                                                                      gtTCACAAATTGATGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGCAAGTTTTGCTGGTATTCATTCAGTTCTAGCACACTGATCACttt   <  2:8943460/93‑1 (MQ=255)

                                                                                       ttCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTa  >  1:8547548/1‑93 (MQ=255)

TTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTA  >  NC_000913/1570504‑1570683

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