BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	715,849	A→G	L146L (TTA→TTG)	ybfP →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	715,849	0	A	G	98.2%	385.3 / ‑2.6	113	L146L (TTA→TTG)	ybfP	lipoprotein
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (48/63); total (48/65)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 120 total aligned reads displayed.

TGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGT  >  NC_000913/715751‑715945

tGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTac                                                                                                 >  1:1970084/1‑100 (MQ=255)

 gCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcg                                                                                                <  2:72907/100‑1 (MQ=255)

  ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcga                                                                                               >  1:2747268/1‑100 (MQ=255)

  ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcga                                                                                               >  2:1300641/1‑100 (MQ=255)

   cAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcgac                                                                                              >  2:1631872/1‑100 (MQ=255)

   cAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcgac                                                                                              >  2:2896096/1‑100 (MQ=255)

     gggCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGc                                                                                            >  2:789990/1‑100 (MQ=255)

     gggCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGc                                                                                            >  2:1294920/1‑100 (MQ=255)

       gCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTg                                                                                          <  1:402734/100‑1 (MQ=255)

       gCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTg                                                                                          >  1:1016814/1‑100 (MQ=255)

        cAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGc                                                                                         >  1:1174065/1‑100 (MQ=255)

        cAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGc                                                                                         >  1:2169617/1‑100 (MQ=255)

        cAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATAGAGTTGCGACGCTGc                                                                                         >  1:2989796/1‑100 (MQ=255)

        cAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGAATTATTGAGTTGCGACGCTGc                                                                                         >  1:998201/1‑100 (MQ=255)

         aaGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGcc                                                                                        <  2:1953716/100‑1 (MQ=255)

            cccGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGCCGGGATTATTGAGTTGCGACGCTGCCAGt                                                                                     >  1:2142497/1‑100 (MQ=255)

              cGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCa                                                                                   <  1:2766122/100‑1 (MQ=255)

               ggCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGCGTTGCGACGCTGCCAGTCAg                                                                                  <  2:1855085/100‑1 (MQ=255)

               ggCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAg                                                                                  >  1:2289403/1‑100 (MQ=255)

               ggCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAg                                                                                  <  1:2565762/100‑1 (MQ=255)

                gCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGa                                                                                 >  1:2567965/1‑100 (MQ=255)

                  aCTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGaaa                                                                               <  1:584206/100‑1 (MQ=255)

                   cTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGCTTATTGAGTTGCGACGCTGCCAGTCAGaaaa                                                                              >  2:2316981/1‑100 (MQ=255)

                   cTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGaaaa                                                                              >  1:2500543/1‑100 (MQ=255)

                   cTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGa                                                                                 <  2:1243571/97‑1 (MQ=255)

                    tCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAAt                                                                             >  2:2265988/1‑100 (MQ=255)

                        ccGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGcc                                                                         <  1:960322/100‑1 (MQ=255)

                        ccGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCCGCCAGTCAGAAAATTGcc                                                                         <  2:91736/100‑1 (MQ=255)

                         cGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCa                                                                        <  1:1539712/100‑1 (MQ=255)

                           aTAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGc                                                                      <  1:1732143/100‑1 (MQ=255)

                           aTAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGc                                                                      <  2:765068/100‑1 (MQ=255)

                            tAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCa                                                                     >  2:1499124/1‑100 (MQ=255)

                            tAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCa                                                                     >  2:2440810/1‑100 (MQ=255)

                                cAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGaa                                                                 <  1:1499136/100‑1 (MQ=255)

                                  aTGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAAtt                                                               <  1:789997/100‑1 (MQ=255)

                                   tGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTg                                                              <  1:2847561/100‑1 (MQ=255)

                                   tGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTg                                                              >  1:2363571/1‑100 (MQ=255)

                                     aGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaa                                                            <  2:2674358/100‑1 (MQ=255)

                                     aGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaa                                                            >  1:2354142/1‑100 (MQ=255)

                                      ggATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaaa                                                           >  2:307229/1‑100 (MQ=255)

                                       gatgatGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAg                                                          <  2:2939849/100‑1 (MQ=255)

                                        aTTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGTGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGc                                                         <  2:2851611/100‑1 (MQ=255)

                                         ttATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCGTGTCGGGATTATTGAGTTGAGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCg                                                        <  2:1660614/100‑1 (MQ=255)

                                           aTGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                                                      <  1:602820/100‑1 (MQ=255)

                                            tGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATa                                                     >  1:2664520/1‑100 (MQ=255)

                                            tGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATa                                                     <  1:2217517/100‑1 (MQ=255)

                                             ntCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGACAATTCCCAGCAAGAATTGAAACCGATaa                                                    >  2:697070/2‑100 (MQ=255)

                                             gTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATaa                                                    <  1:307232/100‑1 (MQ=255)

                                              tCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATaaa                                                   <  1:2611189/100‑1 (MQ=255)

                                               cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                  >  2:910106/1‑100 (MQ=255)

                                                aGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATa                                                 >  2:2239057/1‑100 (MQ=255)

                                                 gCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATac                                                >  1:2648265/1‑100 (MQ=255)

                                                   gAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac                                              >  2:2256817/1‑100 (MQ=255)

                                                    aaaGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca                                             <  1:1447452/100‑1 (MQ=255)

                                                    aaaGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca                                             >  2:1616040/1‑100 (MQ=255)

                                                    aaaGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca                                             <  1:2022260/100‑1 (MQ=255)

                                                      aGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACAtt                                           >  1:1351561/1‑100 (MQ=255)

                                                      aGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACAtt                                           <  2:1636/100‑1 (MQ=255)

                                                        cGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaa                                         <  2:1136973/100‑1 (MQ=255)

                                                         ggATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaaa                                        >  1:1898435/1‑100 (MQ=255)

                                                           aTTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAAt                                      <  1:444092/100‑1 (MQ=255)

                                                             tttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg                                    <  2:1340246/100‑1 (MQ=255)

                                                             tttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg                                    <  2:819188/100‑1 (MQ=255)

                                                             tttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg                                    <  1:290176/100‑1 (MQ=255)

                                                              ttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGt                                   <  2:2500521/100‑1 (MQ=255)

                                                              ttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGt                                   >  1:650324/1‑100 (MQ=255)

                                                               tttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGTAAGCGATAAATACACATTAAAATAGTg                                  <  2:2935428/100‑1 (MQ=255)

                                                               tttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTg                                  <  1:2896122/100‑1 (MQ=255)

                                                                  ggCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGACAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAg                                <  2:1351763/99‑1 (MQ=255)

                                                                   gCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTc                              >  1:1144737/1‑100 (MQ=255)

                                                                    cTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAGTTGCCAGCAAGAATTTAAAGCGATAAATACACCTTAAAATAGTGAGTcc                             >  2:2570676/1‑100 (MQ=255)

                                                                    cTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTcc                             <  1:1547125/100‑1 (MQ=255)

                                                                      cGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCgg                           >  2:2756128/1‑100 (MQ=255)

                                                                       gatgatGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCggg                          <  2:1352117/100‑1 (MQ=255)

                                                                       gatgatGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCggg                          <  1:233263/100‑1 (MQ=255)

                                                                        atgatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGa                         >  1:1935982/1‑100 (MQ=255)

                                                                        atgatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGa                         <  2:1926226/100‑1 (MQ=255)

                                                                          gatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAAt                       <  2:769345/100‑1 (MQ=255)

                                                                          gatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAAt                       <  1:86541/100‑1 (MQ=255)

                                                                          gatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAAt                       >  1:1782426/1‑100 (MQ=255)

                                                                           atgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATg                      <  2:413184/100‑1 (MQ=255)

                                                                           atgCATGTAGAGATTATTGAGTTGCGACGCTGCCAGTAAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATg                      <  2:1452062/100‑1 (MQ=255)

                                                                            tgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGa                     >  1:2473461/1‑100 (MQ=255)

                                                                            tgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGa                     <  1:1215774/100‑1 (MQ=255)

                                                                              cATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGaaa                   >  2:435949/1‑100 (MQ=255)

                                                                               aTGCCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAAc                  <  2:2142478/100‑1 (MQ=255)

                                                                                tGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAAcc                 <  2:755190/100‑1 (MQ=255)

                                                                                 gTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCgggtttcattccc                <  1:2127923/100‑11 (MQ=255)

                                                                                  tCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCg               <  2:1783491/100‑1 (MQ=255)

                                                                                   cGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCg               >  2:1946042/1‑99 (MQ=255)

                                                                                    gggATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGa             <  1:1295720/100‑1 (MQ=255)

                                                                                    gggATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGa             <  2:1174054/100‑1 (MQ=255)

                                                                                       attattGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCt          <  2:817438/100‑1 (MQ=255)

                                                                                        ttattGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTg         <  2:2289383/100‑1 (MQ=255)

                                                                                         tattGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGa        <  1:2266008/100‑1 (MQ=255)

                                                                                          attGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGaa       >  2:2535475/1‑100 (MQ=255)

                                                                                            tGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACa     <  1:415128/100‑1 (MQ=255)

                                                                                               gTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGt  >  2:1714268/1‑100 (MQ=255)

                                                                                               gTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGt  <  2:2810295/100‑1 (MQ=255)

                                                                                               gTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGt  >  1:2647887/1‑100 (MQ=255)

TGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGT  >  NC_000913/715751‑715945

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