BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,987,532	A→C	F76C (TTT→TGT)	yecJ ←	DUF2766 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,987,532	0	A	C	98.6%	473.6 / ‑2.5	139	F76C (TTT→TGT)	yecJ	DUF2766 family protein
Reads supporting (aligned to +/- strand): ref base A (0/2); new base C (70/66); total (70/69)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 141 total aligned reads displayed.

TCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGC  >  NC_000913/1987435‑1987627

tcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTacac                                                                                               <  1:3673560/100‑1 (MQ=255)

 cGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACt                                                                                              >  1:3073268/1‑100 (MQ=255)

  ggCGCTAATGCAGTGCGCCGATGAAGGAATTGAGAAACAGCTGCCTCATCTGGTTGAAGCAGCGGGGTGTTTAATGGGATTCACCCTGTGGGGTACACtt                                                                                             >  2:583068/1‑100 (MQ=255)

   gcgcTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACttn                                                                                            <  1:310570/100‑2 (MQ=255)

   gcgcTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTg                                                                                            <  1:2088558/100‑1 (MQ=255)

    cgcTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGa                                                                                           >  2:506094/1‑100 (MQ=255)

     gcTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGan                                                                                          <  2:3580323/100‑2 (MQ=255)

     gcTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAg                                                                                          <  2:2738986/100‑1 (MQ=255)

     gcTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAg                                                                                          <  2:1219055/100‑1 (MQ=255)

      cTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGt                                                                                         >  1:3105899/1‑100 (MQ=255)

       tAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTGGAAGCAGCGGTTTTTTTAATGTGATTCACCCTGTGGGGTACACTTGAGtt                                                                                        >  1:3875383/1‑100 (MQ=255)

        aaTGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTc                                                                                       <  2:887116/100‑1 (MQ=255)

           gCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAAt                                                                                    >  1:3608291/1‑100 (MQ=255)

           gCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAAt                                                                                    <  2:3458264/100‑1 (MQ=255)

             aTTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATaa                                                                                  >  2:317668/1‑100 (MQ=255)

             aTTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGCTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTAACTaa                                                                                  >  2:2974289/1‑100 (MQ=255)

                gcgcCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAgcg                                                                               <  1:1233375/100‑1 (MQ=255)

                 cgcCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAgcgc                                                                              >  1:2275841/1‑100 (MQ=255)

                   ccGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGa                                                                            >  2:1508091/1‑100 (MQ=255)

                     gATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATg                                                                          <  1:287573/100‑1 (MQ=255)

                      aTGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATgg                                                                         >  1:199001/1‑100 (MQ=255)

                       tGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGc                                                                        <  1:331625/100‑1 (MQ=255)

                        gAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCt                                                                       >  1:4066733/1‑100 (MQ=255)

                        gAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCt                                                                       >  1:1157680/1‑100 (MQ=255)

                          aGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCttt                                                                     >  1:3666889/1‑100 (MQ=255)

                            nttttGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCttttt                                                                   >  1:3259899/2‑100 (MQ=255)

                             ttttGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTg                                                                  <  1:902840/100‑1 (MQ=255)

                             ttttGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTg                                                                  <  2:2440746/100‑1 (MQ=255)

                             ttttGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTg                                                                  <  1:2133216/100‑1 (MQ=255)

                             ttttGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTg                                                                  <  2:2667142/100‑1 (MQ=255)

                              tttGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTgg                                                                 >  1:991515/1‑100 (MQ=255)

                              tttGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTgg                                                                 <  2:1028812/100‑1 (MQ=255)

                               ttGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGa                                                                <  1:1996864/100‑1 (MQ=255)

                                tGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGAt                                                               >  1:270140/1‑100 (MQ=255)

                                tGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGAt                                                               >  1:1630672/1‑100 (MQ=255)

                                tGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGAt                                                               >  1:2952509/1‑100 (MQ=255)

                                 gagaTACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGAtt                                                              >  2:3384558/1‑100 (MQ=255)

                                 gagaAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGAtt                                                              <  2:758566/100‑1 (MQ=255)

                                 aaaaaaCCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGAtt                                                              >  1:3280280/4‑100 (MQ=255)

                                  agaAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTg                                                             >  1:2351786/1‑100 (MQ=255)

                                   gaAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACACTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGc                                                            >  1:1456253/1‑100 (MQ=255)

                                     aaCCGCTGCCTCAGCTGTTTGAATGAGCTGTTTTTTTAATGGGATTCACCCTGTTGGGTACACTTGAGTTCAATCAGCGCGATGGCTTTTTGGATTGccc                                                          <  2:1285854/100‑1 (MQ=255)

                                      aCCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCg                                                         >  1:3601972/1‑100 (MQ=255)

                                        cGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCa                                                       >  2:2595018/1‑100 (MQ=255)

                                        cGCTGCCTCATATGATTGAAGCAGCGGTTTTTTTAATGGGATTCGCCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCa                                                       >  1:3580341/1‑100 (MQ=255)

                                         gCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCAt                                                      >  2:1945193/1‑100 (MQ=255)

                                         gCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCAt                                                      <  1:608802/100‑1 (MQ=255)

                                             ccTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTg                                                  <  1:2563768/100‑1 (MQ=255)

                                              cTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGa                                                 >  2:1628977/1‑100 (MQ=255)

                                                cATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACATTTGAGTTCAATAAGCGCGATAGCTTTTTGGATTGCCCGCATGGTGAcc                                               <  2:3983025/100‑1 (MQ=255)

                                                cATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGAcc                                               <  1:3453913/100‑1 (MQ=255)

                                                cATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGAcc                                               >  1:2942005/1‑100 (MQ=255)

                                                 aTCTGTTTGACGGAGCGGTTGTTTTAATGGGATTCACCCTGTGGGGTACACTTGCGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCg                                              <  2:285485/100‑1 (MQ=255)

                                                   cTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCggg                                            <  2:2697341/100‑1 (MQ=255)

                                                     gTTTGTAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGTGTACACTTGAGTTCATTAAGCGCGATGGCTTCTTGGATTACCCGCATGGTGACCGGGTc                                          >  2:3005019/1‑100 (MQ=255)

                                                     gTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTc                                          <  1:2795503/100‑1 (MQ=255)

                                                       ttGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTg                                        <  2:21485/100‑1 (MQ=255)

                                                        tGAAGCAGCGGTTTTTTTACTGGGATTCACCCTGTGGGGTACACTTGAGTGCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgc                                       <  2:288284/100‑1 (MQ=255)

                                                           agcagcGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcggc                                    >  1:967745/1‑100 (MQ=255)

                                                           agcagcGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcggc                                    <  1:2546794/100‑1 (MQ=255)

                                                              agcGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcggcggc                                 >  1:922574/1‑100 (MQ=255)

                                                              agcGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcggcggc                                 <  2:777252/100‑1 (MQ=255)

                                                               gcGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcggcggcg                                <  1:834680/100‑1 (MQ=255)

                                                                 ggTTTTTTTAATGGGATTCACCCTGCGGGGTACACTTGAGTTCAATGAGCGCGATGGCTTTTTGGGTTGCCCGCATGGTGACCGGGTCTGCGGCGGCggg                              <  1:3786157/100‑1 (MQ=255)

                                                                  gTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGa                             >  1:995765/1‑100 (MQ=255)

                                                                     tttttAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATgg                          <  2:1028695/100‑1 (MQ=255)

                                                                      ttttAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGt                         >  2:3960761/1‑100 (MQ=255)

                                                                      ttttAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGt                         <  1:569859/100‑1 (MQ=255)

                                                                       tttAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGtt                        <  1:2614491/100‑1 (MQ=255)

                                                                       tttAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGtt                        <  1:2615270/100‑1 (MQ=255)

                                                                       tttAAGGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGtt                        <  2:1034023/100‑1 (MQ=255)

                                                                         tAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTTTTTTCATTGTGACCGGGTCTGCGGCGGCGGGATGTTTAg                      >  2:538098/1‑100 (MQ=255)

                                                                         tAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAg                      <  1:2299476/100‑1 (MQ=255)

                                                                         tAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAg                      >  1:855838/1‑100 (MQ=255)

                                                                         tAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAg                      <  2:3342797/100‑1 (MQ=255)

                                                                          aaTGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGt                     >  2:963155/1‑100 (MQ=255)

                                                                           attggATTTACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATGGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTa                    <  2:311632/97‑1 (MQ=255)

                                                                           aTGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTa                    >  1:3963716/1‑100 (MQ=255)

                                                                           aTGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTa                    <  1:1976362/100‑1 (MQ=255)

                                                                           aTGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTa                    <  1:2800113/100‑1 (MQ=255)

                                                                            tGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTaa                   <  1:1534379/100‑1 (MQ=255)

                                                                             gggATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTaaa                  >  1:1746121/1‑100 (MQ=255)

                                                                               gATTCACCCTGTGGGATACACTTGAGTTCACTAAGCGCGACGGCTCTTTGGATTGCCCGCATGGTGACCGGGACTGCGGCGGCGGGATGGTTAGTAAAGt                >  1:3541685/1‑100 (MQ=255)

                                                                                aTTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGCTGGTTAGTAAAGTc               >  2:2154391/1‑100 (MQ=255)

                                                                                aTTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTc               <  2:1178258/100‑1 (MQ=255)

                                                                                   cACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGat            <  1:2719254/100‑1 (MQ=255)

                                                                                   cACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGat            >  1:1024610/1‑100 (MQ=255)

                                                                                   cACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGat            <  1:2087515/100‑1 (MQ=255)

                                                                                   cACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGat            >  1:3414025/1‑100 (MQ=255)

                                                                                    aCCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGata           >  1:368656/1‑100 (MQ=255)

                                                                                     cccTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGatat          <  1:2367951/100‑1 (MQ=255)

                                                                                     cccTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGatat          >  1:1633155/1‑100 (MQ=255)

                                                                                     cccTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGatat          <  1:1405314/100‑1 (MQ=255)

                                                                                      ccTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATAtt         <  1:4041165/100‑1 (MQ=255)

                                                                                       cTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATAttc        >  2:179600/1‑100 (MQ=255)

                                                                                         gtgGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATAttctt      >  1:1923047/1‑100 (MQ=255)

                                                                                            gggTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAg   >  1:2619214/1‑100 (MQ=255)

                                                                                            gggTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAg   >  2:211744/1‑100 (MQ=255)

                                                                                             ggTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGc  >  2:3459943/1‑100 (MQ=255)

                                                                                             ggTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTGCAGc  <  2:3383990/100‑1 (MQ=255)

TCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGC  >  NC_000913/1987435‑1987627

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