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	97.7%	456.7 / ‑0.9	129	L146L (TTA→TTG)	ybfP	lipoprotein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (68/58); total (69/60)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 134 total aligned reads displayed.

CCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGT  >  NC_000913/715753‑715945

 ttccgatctcccGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcgac                                                                                              <  2:2722442/91‑1 (MQ=255)

  aGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATTCATGTCGGGATTATTGAGTTgcgacg                                                                                             <  1:2080881/100‑1 (MQ=255)

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

      cAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGc                                                                                         <  2:1056125/100‑1 (MQ=255)

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

        aGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCa                                                                                       >  1:1647986/1‑100 (MQ=255)

         gCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAg                                                                                      <  2:1509399/100‑1 (MQ=255)

         gCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAg                                                                                      <  1:84635/100‑1 (MQ=255)

          cccGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcgacagatcggaa                                                                                     >  1:2722451/1‑91 (MQ=255)

           ccGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTc                                                                                    <  1:2597376/100‑1 (MQ=255)

           ccGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTc                                                                                    <  2:2473406/100‑1 (MQ=255)

           ccGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTc                                                                                    <  1:3049082/100‑1 (MQ=255)

            cGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCa                                                                                   >  1:2620966/1‑100 (MQ=255)

            cGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCa                                                                                   <  2:1159538/100‑1 (MQ=255)

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

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

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

               cACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGaa                                                                                >  1:1180054/1‑100 (MQ=255)

               cACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGaa                                                                                >  1:3331382/1‑100 (MQ=255)

                 cTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGaaaa                                                                              <  1:2801358/100‑1 (MQ=255)

                   cGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAAtt                                                                            >  2:2693067/1‑100 (MQ=255)

                   cGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAAtt                                                                            <  1:1360096/100‑1 (MQ=255)

                    gACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTg                                                                           <  1:2146140/100‑1 (MQ=255)

                     aCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGc                                                                          >  2:815074/1‑100 (MQ=255)

                     aCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGc                                                                          >  2:1616990/1‑100 (MQ=255)

                     aCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGACTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAAtaga                                                                          >  2:1820726/1‑97 (MQ=255)

                      atcttAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGcc                                                                         <  1:2802375/96‑1 (MQ=255)

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

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

                        gATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAg                                                                       >  2:2027467/1‑100 (MQ=255)

                            gTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAg                                                                   <  1:1042702/100‑1 (MQ=255)

                            gTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAg                                                                   >  2:2673148/1‑100 (MQ=255)

                             tCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGa                                                                  <  1:2455354/100‑1 (MQ=255)

                              cAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGaa                                                                 >  1:2377/1‑100 (MQ=255)

                                aTGAGGATTATGTCAGCGAAAGCGGA‑TTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTg                                                              <  1:1886267/100‑1 (MQ=255)

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

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

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

                                     gATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAg                                                          <  2:1005869/100‑1 (MQ=255)

                                       ttATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCg                                                        >  2:2543078/1‑100 (MQ=255)

                                        tATGTCAGCGAAAGCGGATTTTTTTGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa                                                       <  1:1671749/100‑1 (MQ=255)

                                         aTGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                                                      >  2:3015754/1‑100 (MQ=255)

                                         aTGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                                                      >  1:1196650/1‑100 (MQ=255)

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

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

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

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

                                                cGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATaca                                               >  2:526810/1‑100 (MQ=255)

                                                 gAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac                                              <  2:2122018/100‑1 (MQ=255)

                                                 cgatctggATTTTTTGGCTCGATGATGCATGTCGTGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGTGATCAATacac                                              <  2:1791538/94‑1 (MQ=255)

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

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

                                                    tgCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACAtt                                           <  1:346054/99‑1 (MQ=255)

                                                     gCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTa                                          <  2:373116/100‑1 (MQ=255)

                                                     gCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTa                                          >  2:346053/1‑100 (MQ=255)

                                                      cGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaa                                         >  1:3320676/1‑100 (MQ=255)

                                                      cGGATTTTTTGGCTCGATGATGCATGTCGGAATAATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAGTACACATTaa                                         >  2:2240632/1‑100 (MQ=255)

                                                       ggATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTan                                         <  2:2348767/99‑2 (MQ=255)

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

                                                       ggATTTTTTGGCTCGATCATGCATGTCGGGATTATTGCGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaa                                         <  1:2716772/99‑1 (MQ=255)

                                                       ggATTTTATGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTAAAAGCGATAAATAcacagatcg                                        >  1:1791542/1‑95 (MQ=255)

                                                        gATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaaaa                                       <  2:1196647/100‑1 (MQ=255)

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

                                                          ttttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATa                                     >  1:1457124/1‑100 (MQ=255)

                                                          ttttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATa                                     <  2:1206546/100‑1 (MQ=255)

                                                           tttttGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg                                    >  1:1026895/1‑100 (MQ=255)

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

                                                              atctcTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGa                                 <  2:15625/96‑1 (MQ=255)

                                                                 gCTCGATGATGCATGTCGGGCTTATTGAGTTGCGACGCTGCCCGTCAGAAAATTGCCAGCACGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTc                              >  2:676717/1‑100 (MQ=255)

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

                                                                  cTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGaagat                             >  1:15625/1‑96 (MQ=255)

                                                                  cTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTcc                             >  2:2209567/1‑100 (MQ=255)

                                                                  cTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTcc                             >  2:2144795/1‑100 (MQ=255)

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

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

                                                                         atgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATg                      >  1:1602973/1‑100 (MQ=255)

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

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

                                                                            cATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGaga                   >  1:1561692/1‑98 (MQ=255)

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

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

                                                                             aTGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAAc                  >  1:2808013/1‑100 (MQ=255)

                                                                                 cGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCgg              >  2:404003/1‑100 (MQ=255)

                                                                                  gggATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGa             >  1:992676/1‑100 (MQ=255)

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

                                                                                   ggATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGAc            <  1:1693968/100‑1 (MQ=255)

                                                                                      ttattGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTg         >  1:2190179/1‑100 (MQ=255)

                                                                                      ttattGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTg         >  1:1181642/1‑100 (MQ=255)

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

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

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

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

                                                                                        attGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCCGaa       >  1:727181/1‑100 (MQ=255)

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

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

                                                                                           gAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACa     >  1:2205327/1‑99 (MQ=255)

                                                                                           gAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACAt    >  2:252092/1‑100 (MQ=255)

                                                                                           gAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACAt    >  2:1450555/1‑100 (MQ=255)

                                                                                           gAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACAt    >  2:1450254/1‑100 (MQ=255)

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

CCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGT  >  NC_000913/715753‑715945

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