BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	127,038	G→A	A448A (GCG→GCA)	aceF →	pyruvate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	127,038	0	G	A	99.0%	1064.7 / ‑3.2	299	A448A (GCG→GCA)	aceF	pyruvate dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (1/2); new base A (160/136); total (161/138)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 315 total aligned reads displayed.

TAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGATGCCTCGCTTCAA  >  NC_000913/126948‑127127

  aCTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGa                                                                                       <  1:958641/93‑1 (MQ=255)

  aCTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGa                                                                                       >  2:2561951/1‑93 (MQ=255)

   cTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGaa                                                                                      >  2:5438626/1‑93 (MQ=255)

    tGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGaaa                                                                                     >  2:3328432/1‑93 (MQ=255)

     gggTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAAc                                                                                    <  2:5742257/93‑1 (MQ=255)

     gggTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAAc                                                                                    <  1:3368037/93‑1 (MQ=255)

     gggTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAAc                                                                                    >  1:2849186/1‑93 (MQ=255)

      ggTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACg                                                                                   >  1:3307700/1‑93 (MQ=255)

       gTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGt                                                                                  >  1:717180/1‑93 (MQ=255)

        tAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTa                                                                                 <  2:2968915/93‑1 (MQ=255)

        tAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTa                                                                                 >  1:6908245/1‑93 (MQ=255)

        tAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTa                                                                                 <  2:488244/93‑1 (MQ=255)

         aaTGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTaa                                                                                <  2:195477/93‑1 (MQ=255)

          aTGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAg                                                                               >  2:375286/1‑93 (MQ=255)

          aTGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAg                                                                               <  2:4835506/93‑1 (MQ=255)

           tGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGc                                                                              >  1:5910283/1‑93 (MQ=255)

           tGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGc                                                                              <  1:4202041/93‑1 (MQ=255)

             aTCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTg                                                                            >  1:1333515/1‑93 (MQ=255)

             aTCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTg                                                                            <  1:5118187/93‑1 (MQ=255)

              tCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGACACGTAAGCTgg                                                                           >  1:3987551/1‑93 (MQ=255)

              tCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTgg                                                                           <  2:1616781/93‑1 (MQ=255)

              tCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTgg                                                                           >  1:515853/1‑93 (MQ=255)

               cccGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGa                                                                          >  1:1540425/1‑93 (MQ=255)

               cccGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGa                                                                          <  2:3995729/93‑1 (MQ=255)

                 cGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGAtg                                                                        >  2:6097923/1‑93 (MQ=255)

                   cATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGAtgtg                                                                      <  2:6780289/93‑1 (MQ=255)

                   cATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGAtgtg                                                                      >  2:2326728/1‑93 (MQ=255)

                    aTGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGa                                                                     <  2:2781804/93‑1 (MQ=255)

                    aTGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGa                                                                     <  1:4460326/93‑1 (MQ=255)

                    aTGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGa                                                                     >  1:2361745/1‑93 (MQ=255)

                    aTGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGa                                                                     >  1:6395616/1‑93 (MQ=255)

                     tGTTCCTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGaa                                                                    <  1:7050414/93‑1 (MQ=255)

                     tGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGaa                                                                    >  1:7400693/1‑93 (MQ=255)

                     tGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGaa                                                                    <  2:4711925/93‑1 (MQ=255)

                     tGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGaa                                                                    >  2:1841412/1‑93 (MQ=255)

                       ttACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGa                                                                  <  2:3008785/93‑1 (MQ=255)

                        tACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGAt                                                                 <  2:842119/93‑1 (MQ=255)

                         aCTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATc                                                                >  1:5862558/1‑93 (MQ=255)

                          cTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCa                                                               <  2:3449087/93‑1 (MQ=255)

                          cTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCa                                                               <  1:4816856/93‑1 (MQ=255)

                          cTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCa                                                               <  2:6447926/93‑1 (MQ=255)

                            cACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCAcc                                                             <  2:7705584/93‑1 (MQ=255)

                            cACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCAcc                                                             >  1:7842564/1‑93 (MQ=255)

                             aCTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCAccc                                                            >  1:7355661/1‑93 (MQ=255)

                             aCTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCAccc                                                            >  1:6947096/1‑93 (MQ=255)

                              cTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCAcccc                                                           <  1:3535553/93‑1 (MQ=255)

                              cTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCAcccc                                                           >  2:6638226/1‑93 (MQ=255)

                              cTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCAcccc                                                           >  1:2636068/1‑93 (MQ=255)

                               ttCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTCAGCTGGATGTGAAGATCACCCCg                                                          >  2:5649837/1‑93 (MQ=255)

                               ttCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCg                                                          >  2:1488263/1‑93 (MQ=255)

                                tCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCgg                                                         <  2:3962882/93‑1 (MQ=255)

                                  gACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGtt                                                       >  1:3222044/1‑93 (MQ=255)

                                  gACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGtt                                                       >  2:3509727/1‑93 (MQ=255)

                                   aCAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTg                                                      <  2:877156/93‑1 (MQ=255)

                                      aaaCCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCt                                                   <  2:7912524/93‑1 (MQ=255)

                                         ccGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTtca                                                >  1:1840135/1‑93 (MQ=255)

                                         ccGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTtca                                                <  1:4542903/93‑1 (MQ=255)

                                            atatCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTtcatca                                             <  2:3878633/93‑1 (MQ=255)

                                             tatCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTcatcat                                            >  1:848204/1‑93 (MQ=255)

                                               tCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGa                                          >  2:2250906/1‑93 (MQ=255)

                                               tCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGa                                          >  2:459704/1‑93 (MQ=255)

                                                cACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGaa                                         <  2:3087490/93‑1 (MQ=255)

                                                  ccGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAg                                       <  1:7782021/93‑1 (MQ=255)

                                                   cGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGc                                      >  2:322878/1‑93 (MQ=255)

                                                    gAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGcc                                     >  1:5711328/1‑93 (MQ=255)

                                                    gAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGcc                                     >  2:6269080/1‑93 (MQ=255)

                                                      gTTGGAAGAGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGACGt                                   >  2:6754034/1‑93 (MQ=255)

                                                       ttGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGtt                                  <  1:824237/93‑1 (MQ=255)

                                                         ggAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTtgc                                >  2:7273354/1‑93 (MQ=255)

                                                         ggAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTtgc                                >  2:6461730/1‑93 (MQ=255)

                                                          gAAGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTtgct                               >  1:123807/1‑93 (MQ=255)

                                                            aGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTtgctgc                             <  1:6825844/93‑1 (MQ=255)

                                                            aGCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACATAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTtgctgc                             >  1:1009087/1‑93 (MQ=255)

                                                             gCGTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCa                            <  1:2239770/93‑1 (MQ=255)

                                                               gTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGc                          <  2:4397571/93‑1 (MQ=255)

                                                               gTTCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGc                          >  1:3833833/1‑93 (MQ=255)

                                                               gTTCCGGAAACAGGCGAACGAAGGAGCAGCGAAACGTAAGGTTGATGTGAAGATCCCCCCAGTTGTCTTCATCATGAAAGCCGTTGCTGCAGc                          <  1:555734/93‑1 (MQ=255)

                                                                ttCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGct                         >  1:7505621/1‑93 (MQ=255)

                                                                 tCCGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGctc                        >  2:4258008/1‑93 (MQ=255)

                                                                   cGTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCtt                      <  1:3718624/93‑1 (MQ=255)

                                                                    gTAAACAGCAGAACGAAGAAGCAGCGTAACGTAAGCTGGATGTGAAGTTCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTg                     >  2:4682154/1‑93 (MQ=255)

                                                                    gTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGNAAGCCGTTGCTGCAGCTCTTg                     >  1:6054444/1‑93 (MQ=255)

                                                                    gTAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTg                     <  2:5215544/93‑1 (MQ=255)

                                                                     tAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGa                    <  1:4492686/93‑1 (MQ=255)

                                                                     tAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGa                    >  2:2041706/1‑93 (MQ=255)

                                                                     tAAACAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGa                    <  2:7001299/93‑1 (MQ=255)

                                                                      aaaCAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAg                   <  2:2180199/93‑1 (MQ=255)

                                                                        aCAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCa                 >  2:133490/1‑93 (MQ=255)

                                                                        aCAGCAGAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCa                 >  1:5822803/1‑93 (MQ=255)

                                                                          agcagAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGa               <  2:2585634/93‑1 (MQ=255)

                                                                           gcagAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGAt              <  2:5943566/93‑1 (MQ=255)

                                                                           gcagAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGAt              <  2:1359739/93‑1 (MQ=255)

                                                                              gAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGATGcc           >  1:2312760/1‑93 (MQ=255)

                                                                              gAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGATGcc           <  2:6474996/93‑1 (MQ=255)

                                                                              gAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGATGcc           >  1:3911006/1‑93 (MQ=255)

                                                                              gAACGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGATGcc           <  2:5139817/93‑1 (MQ=255)

                                                                               aaCGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGATGCCt          <  2:6506771/93‑1 (MQ=255)

                                                                               aaCGAAGAAGCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGATGCCt          <  1:531263/93‑1 (MQ=255)

                                                                                   aagaagCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGACTTCATCATCAAAGCCGTTGCTCCAGCTCTTGAGCAGATGTCCCGCt      >  2:5801520/1‑93 (MQ=255)

                                                                                     gaagCAGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGATGCCTCGCTTc    <  1:2727335/93‑1 (MQ=255)

TAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGATGCCTCGCTTCAA  >  NC_000913/126948‑127127

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