BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	715,277	A→G	Y86H (TAC→CAC)	ybfG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	715,277	0	A	G	98.6%	449.9 / ‑3.1	145	Y86H (TAC→CAC)	ybfG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (69/74); total (71/74)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 145 total aligned reads displayed.

CTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTATCTAAATCAC  >  NC_000913/715187‑715366

ctttcATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATATCTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGa                                                                                         <  1:3061922/93‑1 (MQ=255)

 tttcATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGaa                                                                                        >  1:1064137/1‑93 (MQ=255)

 tttcATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGaa                                                                                        >  1:3707919/1‑93 (MQ=255)

  ttcATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAAc                                                                                       <  2:895057/93‑1 (MQ=255)

  ttcATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAAc                                                                                       >  1:149393/1‑93 (MQ=255)

  ttcATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAAc                                                                                       <  1:2322411/93‑1 (MQ=255)

   tcATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACg                                                                                      >  2:301353/1‑93 (MQ=255)

    cATGTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGt                                                                                     <  2:4323588/93‑1 (MQ=255)

     aTTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTg                                                                                    >  1:3917938/1‑93 (MQ=255)

        tGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAg                                                                                 <  2:3770270/93‑1 (MQ=255)

        tGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAg                                                                                 <  1:4347169/93‑1 (MQ=255)

          aTACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGc                                                                               <  2:1634230/93‑1 (MQ=255)

           tACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGCGAAGGCt                                                                              >  1:1745141/1‑93 (MQ=255)

            aCTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCtt                                                                             <  1:224518/93‑1 (MQ=255)

             cTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTc                                                                            <  2:1066913/93‑1 (MQ=255)

             cTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTc                                                                            <  1:2470756/93‑1 (MQ=255)

             cTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTc                                                                            <  2:2701019/93‑1 (MQ=255)

             cTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCAATAAGAAAGGGTGGAACGTGAAGGCTTc                                                                            <  1:1212610/93‑1 (MQ=255)

              tGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTcc                                                                           >  2:4151760/1‑93 (MQ=255)

               gCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCg                                                                          >  1:269112/1‑93 (MQ=255)

                cATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGt                                                                         <  1:2711062/93‑1 (MQ=255)

                 aTTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGtt                                                                        <  1:1995437/93‑1 (MQ=255)

                  tttCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTg                                                                       >  1:2240318/1‑93 (MQ=255)

                   ttCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGc                                                                      >  2:661525/1‑93 (MQ=255)

                       cAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGAt                                                                  <  1:691557/93‑1 (MQ=255)

                        agagTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATa                                                                 >  2:1321594/1‑93 (MQ=255)

                         gagTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATaa                                                                <  2:768/93‑1 (MQ=255)

                          agTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAAt                                                               <  2:2382751/93‑1 (MQ=255)

                          agTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAAt                                                               <  2:3681347/93‑1 (MQ=255)

                          agTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAAt                                                               <  1:683448/93‑1 (MQ=255)

                           gTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAAtt                                                              <  2:1521341/93‑1 (MQ=255)

                           gTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAAtt                                                              >  2:3847166/1‑93 (MQ=255)

                              cAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTtta                                                           <  2:1833459/93‑1 (MQ=255)

                               aGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTttat                                                          <  1:3706749/93‑1 (MQ=255)

                               aGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTttat                                                          <  2:431047/93‑1 (MQ=255)

                                gTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTttatt                                                         >  2:685672/1‑93 (MQ=255)

                                gTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTttatt                                                         >  1:3152888/1‑93 (MQ=255)

                                gTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTttata                                                         >  1:4686638/1‑92 (MQ=255)

                                 tAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGGCGATAATTttatta                                                        >  1:3323790/1‑93 (MQ=255)

                                 tAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTttatta                                                        <  2:1223614/93‑1 (MQ=255)

                                 tAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTttatta                                                        >  1:2977856/1‑93 (MQ=255)

                                  aaGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTAc                                                       >  2:2121225/1‑93 (MQ=255)

                                    gCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGt                                                     <  1:3213691/93‑1 (MQ=255)

                                    gCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAACGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGt                                                     >  2:2200435/1‑93 (MQ=255)

                                       aaaTAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGGAGGCTTCCGTTGCCGATAATTTTATTACGtctt                                                  >  2:937282/1‑93 (MQ=255)

                                       aaaTAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGtctt                                                  >  2:3750668/1‑93 (MQ=255)

                                       aaaTAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGtctt                                                  >  1:2193554/1‑93 (MQ=255)

                                       aaaTAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGtctt                                                  >  1:512181/1‑93 (MQ=255)

                                       aaaTAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGtctt                                                  >  1:1235329/1‑93 (MQ=255)

                                        aaTAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGtcttc                                                 <  2:515904/93‑1 (MQ=255)

                                         aTAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGtcttct                                                >  1:2340064/1‑93 (MQ=255)

                                           aaCTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCg                                              <  1:4581238/93‑1 (MQ=255)

                                             cTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGcc                                            >  1:3455290/1‑93 (MQ=255)

                                              ttttCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGcct                                           >  2:4333801/1‑92 (MQ=255)

                                                ttcttgttcttgCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCAtt                                         <  2:3529300/93‑1 (MQ=255)

                                                 tcttgttcttgCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATtg                                        >  1:336470/1‑93 (MQ=255)

                                                 tcttgttcttgCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATtg                                        >  1:1794347/1‑93 (MQ=255)

                                                     gttcttgCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATtgctgc                                    >  2:318429/1‑93 (MQ=255)

                                                      ttcttgCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGcc                                   >  1:1684499/1‑93 (MQ=255)

                                                       tcttgCTTCTTCACTTTGTACGCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCa                                  >  1:4648189/1‑93 (MQ=255)

                                                         ttgCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATc                                >  2:2060340/1‑93 (MQ=255)

                                                         ttgCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATc                                >  1:602424/1‑93 (MQ=255)

                                                         ttgCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATc                                >  1:2486277/1‑93 (MQ=255)

                                                           gCTTCTTCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGt                              >  2:3046987/1‑93 (MQ=255)

                                                            cttcttCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTa                             >  2:2113132/1‑93 (MQ=255)

                                                            cttcttCACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTa                             <  1:2052612/93‑1 (MQ=255)

                                                             ttcttcACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTAc                            <  2:177565/93‑1 (MQ=255)

                                                             ttcttcACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTAc                            <  1:1091062/93‑1 (MQ=255)

                                                             ttcttcACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTAc                            <  2:2012399/93‑1 (MQ=255)

                                                             ttcttcACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTAc                            <  1:3898976/93‑1 (MQ=255)

                                                             ttcttcACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTAc                            <  1:2668867/93‑1 (MQ=255)

                                                                ttcACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACAtt                         >  1:2812382/1‑93 (MQ=255)

                                                                ttcACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACAtt                         >  1:4433057/1‑93 (MQ=255)

                                                                 tcACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACAttt                        >  1:2991792/1‑93 (MQ=255)

                                                                  cACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACAtttt                       <  1:629024/93‑1 (MQ=255)

                                                                  cACTTTGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACAtttt                       >  1:3996901/1‑93 (MQ=255)

                                                                  cACTTTGTACACATAAGAAAGGGTGGAACGTGAAGTCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACAtttt                       <  1:2815806/93‑1 (MQ=255)

                                                                     tttGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCAt                    >  2:2356145/1‑93 (MQ=255)

                                                                      ttGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATa                   >  2:108505/1‑93 (MQ=255)

                                                                       tGTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAg                  <  2:1224874/93‑1 (MQ=255)

                                                                        gTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGa                 <  1:1008566/93‑1 (MQ=255)

                                                                        gTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGa                 <  2:1843259/93‑1 (MQ=255)

                                                                        gTACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGa                 >  1:3547380/1‑93 (MQ=255)

                                                                         tACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGaa                <  2:3516788/93‑1 (MQ=255)

                                                                         tACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGaa                <  2:674074/93‑1 (MQ=255)

                                                                         tACCCATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGaa                <  1:590115/93‑1 (MQ=255)

                                                                           cccATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAgg              <  2:3779343/93‑1 (MQ=255)

                                                                           cccATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAgg              <  1:4033220/93‑1 (MQ=255)

                                                                           cccATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAgg              >  2:218739/1‑93 (MQ=255)

                                                                            ccATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGt             <  2:3909686/93‑1 (MQ=255)

                                                                            ccATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGt             >  1:3592561/1‑93 (MQ=255)

                                                                             cATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTa            <  2:3716854/93‑1 (MQ=255)

                                                                             cATAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTa            <  2:615050/93‑1 (MQ=255)

                                                                               tAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTAtt          >  2:324292/1‑92 (MQ=255)

                                                                               tAAGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTATc          <  2:2241812/93‑1 (MQ=255)

                                                                                aaGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTATCt         <  1:3048150/93‑1 (MQ=255)

                                                                                 aGAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTATCTa        <  2:2812854/93‑1 (MQ=255)

                                                                                  gAAAGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTATCTaa       >  2:4087483/1‑93 (MQ=255)

                                                                                   aaaGGGTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTATCTaaa      >  1:4401831/1‑93 (MQ=255)

                                                                                      gggTGGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTATCTAAATCa   >  2:1332721/1‑93 (MQ=255)

CTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTATCTAAATCAC  >  NC_000913/715187‑715366

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