BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	277,917	A→G	L41P (CTG→CCG)	afuC ←	AfuC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	277,917	0	A	G	99.3%	475.2 / NA	152	L41P (CTG→CCG)	afuC	AfuC
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (64/87); minor base T (1/0); total (65/87)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 162 total aligned reads displayed.

CCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTC  >  NC_000913/277827‑278005

ccATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGagcggc                                                                                        >  2:3010892/1‑93 (MQ=255)

ccATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGagcggc                                                                                        >  1:3260042/1‑93 (MQ=255)

ccATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGagcggc                                                                                        <  2:3631186/93‑1 (MQ=255)

  aTCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGTACGGGCCGAGCGGCGt                                                                                      <  2:2519092/93‑1 (MQ=255)

  aTCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGt                                                                                      <  2:1563894/93‑1 (MQ=255)

   tCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTc                                                                                     >  2:2894545/1‑93 (MQ=255)

    cAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCa                                                                                    >  1:371454/1‑93 (MQ=255)

       tGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCa                                                                                 <  2:755366/93‑1 (MQ=255)

       tGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGCGCCGAGCGGCGTCACCa                                                                                 >  2:3496787/1‑93 (MQ=255)

         aaaaTTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCAtt                                                                               <  2:2141395/93‑1 (MQ=255)

         aaaaTTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCAtt                                                                               >  1:338771/1‑93 (MQ=255)

         aaaaTTTCCCCCTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCAtt                                                                               <  1:1344163/93‑1 (MQ=39)

           aaTTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTg                                                                             >  1:3789710/1‑93 (MQ=255)

            aTTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGc                                                                            >  1:4394472/1‑93 (MQ=255)

             tttGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGcc                                                                           >  2:4179728/1‑93 (MQ=255)

             tttGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGcc                                                                           <  1:311942/93‑1 (MQ=255)

                tcccTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCt                                                                        <  1:2666090/92‑1 (MQ=255)

                   cTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCg                                                                     >  1:1691817/1‑93 (MQ=255)

                    ttCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCgg                                                                    <  1:1377868/93‑1 (MQ=255)

                       gCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGag                                                                 >  2:1707822/1‑92 (MQ=255)

                       gCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGAt                                                                 <  2:64663/93‑1 (MQ=255)

                        cTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATg                                                                >  1:2617683/1‑93 (MQ=255)

                         tCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATgg                                                               >  1:4601346/1‑93 (MQ=255)

                          cGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGt                                                              >  1:457546/1‑93 (MQ=255)

                            gTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGa                                                            >  1:3074039/1‑93 (MQ=255)

                             tttttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTTCCGCAGCCGGCCGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAg                                                           <  2:891289/93‑1 (MQ=255)

                             tttttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAg                                                           <  1:1966817/93‑1 (MQ=255)

                              ttttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAgg                                                          <  1:868197/93‑1 (MQ=255)

                              ttttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAgg                                                          >  2:2731813/1‑93 (MQ=255)

                              ttttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAgg                                                          <  2:22411/93‑1 (MQ=255)

                               tttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGt                                                         >  2:1654009/1‑93 (MQ=255)

                               tttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGt                                                         <  1:3207528/93‑1 (MQ=255)

                               tttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGt                                                         <  1:199948/93‑1 (MQ=255)

                               tttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGt                                                         >  1:688873/1‑93 (MQ=255)

                               tttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGt                                                         >  2:3865727/1‑93 (MQ=255)

                               tttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGt                                                         <  2:2480691/93‑1 (MQ=255)

                                 tCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTg                                                       <  2:2503491/93‑1 (MQ=255)

                                  ccAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGa                                                      >  2:1908617/1‑93 (MQ=255)

                                  ccAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGa                                                      <  2:1321490/93‑1 (MQ=255)

                                  ccAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGa                                                      >  1:2379827/1‑93 (MQ=255)

                                    aGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                    <  1:4442698/93‑1 (MQ=255)

                                    aGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                    <  2:844559/93‑1 (MQ=255)

                                    aGCCCGGCAACCAGGCGCAAAATAGTGGTTTCGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                    <  2:1497946/93‑1 (MQ=255)

                                     gCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGatat                                                   <  2:2382121/93‑1 (MQ=255)

                                     gCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGatat                                                   <  1:2297176/93‑1 (MQ=255)

                                     gCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGGTGatat                                                   >  2:3473130/1‑93 (MQ=255)

                                     gCCCGACAACCAGGAGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCGCTGCGGGATGGTGAGGTTGatat                                                   <  2:4216330/93‑1 (MQ=255)

                                      cccGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATAtt                                                  <  1:450302/93‑1 (MQ=255)

                                       ccGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTg                                                 >  2:3042310/1‑93 (MQ=255)

                                        cGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGt                                                >  1:1631706/1‑93 (MQ=255)

                                         ggCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTc                                               >  1:2220332/1‑93 (MQ=255)

                                          gCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCg                                              >  1:1256287/1‑93 (MQ=255)

                                           cAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGa                                             >  2:3900040/1‑93 (MQ=255)

                                              ccAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTa                                          >  1:3583809/1‑93 (MQ=255)

                                              ccAGGCGCAAAATAGTGGTTTTCCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTa                                          <  2:1551325/93‑1 (MQ=255)

                                                aGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTAcc                                        <  2:3667321/93‑1 (MQ=255)

                                                 ggCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCg                                       >  1:4527158/1‑93 (MQ=255)

                                                   cgcAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTa                                     <  1:3367015/93‑1 (MQ=255)

                                                    gcAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTAt                                    <  1:3024135/93‑1 (MQ=255)

                                                    gcAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTAt                                    <  1:3805701/93‑1 (MQ=255)

                                                    gcAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTAt                                    >  2:1794140/1‑93 (MQ=255)

                                                     cAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTAtt                                   <  2:4533707/93‑1 (MQ=255)

                                                     cAAAATAGTGGTTGTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTAtt                                   <  2:578240/93‑1 (MQ=255)

                                                      aaaaTAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTa                                  <  2:3072536/93‑1 (MQ=255)

                                                       aaaTAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTAc                                 >  2:4391617/1‑93 (MQ=255)

                                                        aaTAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTAACGTATTACt                                >  1:4350118/1‑93 (MQ=255)

                                                         aTAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATCTTGTCGATTACCGTATTACTg                               <  1:2274837/93‑1 (MQ=255)

                                                          tAGTGTTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGc                              >  1:4392021/1‑93 (MQ=255)

                                                          tAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGc                              >  2:364812/1‑93 (MQ=255)

                                                             tGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCaa                           >  2:1116795/1‑93 (MQ=255)

                                                             tGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCaa                           >  2:4254791/1‑93 (MQ=255)

                                                             tGGTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCaa                           <  2:4258309/93‑1 (MQ=255)

                                                              ggTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCaaa                          <  1:4156773/93‑1 (MQ=255)

                                                               gTTTTGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAAt                         >  2:3716928/1‑93 (MQ=255)

                                                                ttttGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATc                        <  1:2785350/93‑1 (MQ=255)

                                                                  ttGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGt                      <  1:4295347/93‑1 (MQ=255)

                                                                  ttGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGt                      >  1:2583471/1‑93 (MQ=255)

                                                                   tGCCGCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGtt                     <  2:81442/93‑1 (MQ=255)

                                                                     acGCAGCCGGACGGGCCTAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTCCCGTATTACTGCCAAATCGTTTa                   <  2:3106584/92‑1 (MQ=255)

                                                                       gCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGt                 <  1:2815073/93‑1 (MQ=255)

                                                                       gCAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGCGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGt                 >  2:3102807/1‑93 (MQ=255)

                                                                        cAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTg                >  2:4469497/1‑93 (MQ=255)

                                                                        cAGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTg                <  1:6980/93‑1 (MQ=255)

                                                                         aGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGa               >  1:3614708/1‑93 (MQ=255)

                                                                         aGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGa               <  1:4108415/93‑1 (MQ=255)

                                                                         aGCCGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGa               >  2:4091058/1‑93 (MQ=255)

                                                                           ccGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACg             >  1:3917190/1‑93 (MQ=255)

                                                                            cagACGGGCAGAGCGGCGTCACCATATGCCCCTGCGGGATGGTGAGGTTGATAGTGTCGCTTACCGTCTTCCTGCCAACTCGTTTAGTGACGt            <  1:4682732/91‑1 (MQ=255)

                                                                            cGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGt            <  1:2525283/93‑1 (MQ=255)

                                                                            cGGACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGt            <  1:2943595/93‑1 (MQ=255)

                                                                             ggAGGGGCCGAGCTGCGTCACCATTTGCCCCTGCGGGATGATGAGGTTGATATTTTCGATTAGCGTATTACTGCCACATCGTTTAGTGACgct           >  2:3310138/1‑91 (MQ=255)

                                                                             ggACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGtt           <  2:3760721/93‑1 (MQ=255)

                                                                              gACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGATGTTg          >  1:2493892/1‑93 (MQ=255)

                                                                              gACGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTg          >  2:1816576/1‑93 (MQ=255)

                                                                               aCGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTgc         <  2:294655/93‑1 (MQ=255)

                                                                               aCGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTgc         <  1:4252733/93‑1 (MQ=255)

                                                                                cGGGCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTgcg        <  1:1498494/93‑1 (MQ=255)

                                                                                  ggCCGAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCa      >  1:1100600/1‑93 (MQ=255)

                                                                                    ccGAGGGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGt    <  2:2433337/93‑1 (MQ=255)

                                                                                      gAGCGGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTc  <  2:3960595/93‑1 (MQ=255)

CCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTC  >  NC_000913/277827‑278005

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