BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,486,485	T→C	Y248C (TAC→TGC)	yfdE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,486,485	0	T	C	91.3%	968.8 / 62.8	357	Y248C (TAC→TGC)	yfdE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (9/22); new base C (151/175); total (160/197)
Fisher's exact test for biased strand distribution p-value = 8.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 362 total aligned reads displayed.

GCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAA  >  NC_000913/2486395‑2486576

  cTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAggg                                                                                         <  1:5848962/93‑1 (MQ=255)

  cTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAggg                                                                                         <  2:1712150/93‑1 (MQ=255)

   tGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGa                                                                                        >  1:9273153/1‑93 (MQ=255)

   tGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGa                                                                                        <  1:2897880/93‑1 (MQ=255)

    ggCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGCGCCATGCAGGGAt                                                                                       >  1:5781684/1‑93 (MQ=255)

     gCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATg                                                                                      <  1:2014814/93‑1 (MQ=255)

     gCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATg                                                                                      <  2:5970226/93‑1 (MQ=255)

      cATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATgg                                                                                     <  2:2004748/93‑1 (MQ=255)

       aTAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGc                                                                                    >  1:5682916/1‑93 (MQ=255)

        tAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCg                                                                                   >  1:1232505/1‑93 (MQ=255)

        tAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCg                                                                                   >  1:2082571/1‑93 (MQ=255)

         aaCGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGa                                                                                  >  1:3717684/1‑93 (MQ=255)

          aCGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGAt                                                                                 <  2:4319302/93‑1 (MQ=255)

            gCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGAttt                                                                               <  1:6680794/93‑1 (MQ=255)

              aGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTcc                                                                             <  1:6015381/93‑1 (MQ=255)

              aGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTcc                                                                             >  1:175022/1‑93 (MQ=255)

               gAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTccc                                                                            >  2:2213123/1‑93 (MQ=255)

               gAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTccc                                                                            <  2:5192343/93‑1 (MQ=255)

                  aaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGa                                                                         <  2:2456338/93‑1 (MQ=255)

                   aaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGAc                                                                        >  1:7137116/1‑93 (MQ=255)

                     aGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGt                                                                      <  2:9583659/93‑1 (MQ=255)

                        ttGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTtgg                                                                   >  2:8658455/1‑92 (MQ=255)

                         tGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTtgtg                                                                  <  1:5983574/93‑1 (MQ=255)

                         tGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTtgtg                                                                  >  2:8289582/1‑93 (MQ=255)

                         tGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTtgtg                                                                  >  1:2350229/1‑93 (MQ=255)

                          ttCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTgtgg                                                                 <  1:8551053/92‑1 (MQ=255)

                          gTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGCTGGCGATTTCCCAGACGTTgtgg                                                                 >  2:5085071/1‑93 (MQ=255)

                          gTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTgtgg                                                                 >  1:8714533/1‑93 (MQ=255)

                           tCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTgtggt                                                                <  1:4356112/93‑1 (MQ=255)

                           tCATGACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTACCCAGACGTTgtggt                                                                <  2:386107/93‑1 (MQ=255)

                             aTTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGa                                                              <  2:7382743/93‑1 (MQ=255)

                              ttACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGAc                                                             >  2:6741010/1‑93 (MQ=255)

                              ttACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGAc                                                             >  2:500538/1‑93 (MQ=255)

                              ttACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGAc                                                             >  1:6866453/1‑93 (MQ=255)

                               tACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACt                                                            >  2:5744239/1‑93 (MQ=255)

                               tACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACt                                                            <  2:8847906/93‑1 (MQ=255)

                                aCCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACtt                                                           >  2:6210437/1‑93 (MQ=255)

                                   acaacaAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTccc                                                        >  2:3758437/1‑93 (MQ=255)

                                    caacaaATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCa                                                       <  1:3919068/93‑1 (MQ=255)

                                     aacaaATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAg                                                      >  1:5492323/1‑93 (MQ=255)

                                      acaaATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGt                                                     <  1:3579579/93‑1 (MQ=255)

                                      acaaATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGt                                                     <  2:6653084/93‑1 (MQ=255)

                                        aaaTCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTg                                                   >  2:7115159/1‑93 (MQ=255)

                                          aTCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCg                                                 >  1:8054459/1‑93 (MQ=255)

                                           tCGTAATCGGCTTATCCTGAGGATTGAAAACATCAAAAGGTGCAATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTgaga                                                >  1:7599055/1‑93 (MQ=255)

                                            cGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGat                                               <  1:4204668/93‑1 (MQ=255)

                                            cGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGat                                               <  2:6866560/93‑1 (MQ=255)

                                             gTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGata                                              >  1:676495/1‑93 (MQ=255)

                                             gTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGata                                              >  1:1087297/1‑93 (MQ=255)

                                              tAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGatat                                             <  1:5300444/93‑1 (MQ=255)

                                               aaTCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGatata                                            <  2:8084664/93‑1 (MQ=255)

                                               aaTCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGatata                                            <  1:9171980/93‑1 (MQ=255)

                                               aaTCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGatata                                            <  1:6395157/93‑1 (MQ=255)

                                               aaTCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGatata                                            <  1:7040088/93‑1 (MQ=255)

                                                aTCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGatatat                                           >  2:1151150/1‑93 (MQ=255)

                                                 tCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATg                                          <  2:36457/93‑1 (MQ=255)

                                                 tCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATg                                          >  1:8391274/1‑93 (MQ=255)

                                                  cGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGc                                         <  2:9752420/93‑1 (MQ=255)

                                                   ggCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGcc                                        >  2:798116/1‑93 (MQ=255)

                                                   ggCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGcc                                        <  1:60896/93‑1 (MQ=255)

                                                    gCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCa                                       <  1:993140/93‑1 (MQ=255)

                                                    gCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCa                                       >  1:7437545/1‑93 (MQ=255)

                                                     cTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCAt                                      <  1:4690690/93‑1 (MQ=255)

                                                     cTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCAt                                      <  1:626565/93‑1 (MQ=255)

                                                      ttATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATc                                     <  1:5127683/93‑1 (MQ=255)

                                                      ttATCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGAGGTGACTTCCCAGTTGCGATATATGCCATc                                     <  2:5440741/93‑1 (MQ=255)

                                                         tCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGa                                  <  2:7833895/93‑1 (MQ=255)

                                                         tCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGa                                  <  1:7884952/93‑1 (MQ=255)

                                                         tCCTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGa                                  <  1:6167783/93‑1 (MQ=255)

                                                          ccTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCTTCCGAc                                 >  1:1369748/1‑93 (MQ=255)

                                                          ccTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGAc                                 <  1:9070791/93‑1 (MQ=255)

                                                           cTGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGTCATCAGAcc                                <  1:6512612/93‑1 (MQ=255)

                                                            tGAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCa                               <  1:9459962/93‑1 (MQ=255)

                                                             gAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCAt                              <  1:2538666/93‑1 (MQ=255)

                                                             gAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCAt                              <  1:1404561/93‑1 (MQ=255)

                                                             gAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCAt                              >  2:1652005/1‑93 (MQ=255)

                                                             gAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCAt                              >  1:3120299/1‑93 (MQ=255)

                                                             gAGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCAt                              <  1:4049308/93‑1 (MQ=255)

                                                              aGTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATg                             <  1:1799466/93‑1 (MQ=255)

                                                               gTATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGc                            <  1:1556887/93‑1 (MQ=255)

                                                                tATTGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCt                           >  2:4695599/1‑93 (MQ=255)

                                                                tATTGAAAACATCAAAAGGCGCCATGCAGGGATGGCGCTTTCCCAGACGCTGAGTTGACTTCCCAGTTGCGATATATTCCATCACACCATGCt                           >  2:2870458/1‑93 (MQ=255)

                                                                 aTTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTc                          >  2:7656532/1‑93 (MQ=255)

                                                                  ttGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTcc                         <  1:5175442/93‑1 (MQ=255)

                                                                  ttGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTcc                         <  2:1231399/93‑1 (MQ=255)

                                                                   tGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCa                        >  2:2635758/1‑93 (MQ=255)

                                                                   tGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCa                        <  1:490596/93‑1 (MQ=255)

                                                                   tGAAAACATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCa                        >  1:8255597/1‑93 (MQ=255)

                                                                     aaaaCATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGa                      >  1:9591443/1‑93 (MQ=255)

                                                                        aCATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGaaaa                   <  2:10034699/93‑1 (MQ=255)

                                                                        aCATCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGaaaa                   <  2:1342651/93‑1 (MQ=255)

                                                                           tCAAAAGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTc                >  2:9313309/1‑93 (MQ=255)

                                                                            cAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCa               <  1:9760745/93‑1 (MQ=255)

                                                                             aaaaGGTGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCATTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAg              >  1:6448904/1‑93 (MQ=255)

                                                                                   tGCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGc        <  2:6132318/93‑1 (MQ=255)

                                                                                   ggCCATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGc        <  2:1795693/92‑1 (MQ=255)

                                                                                      cATGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATc     <  2:9400428/93‑1 (MQ=255)

                                                                                       aTGCAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCa    <  2:4138514/93‑1 (MQ=255)

                                                                                         gTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCaaa  <  2:3968335/93‑1 (MQ=255)

                                                                                         atAGGGATCAAGATTGCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCaaa  <  2:8215984/92‑1 (MQ=255)

GCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAA  >  NC_000913/2486395‑2486576

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