BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,017,921	G→A	C621Y (TGC→TAC)	ygfK →	putative oxidoreductase, Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,017,921	0	G	A	99.6%	1107.2 / NA	275	C621Y (TGC→TAC)	ygfK	putative oxidoreductase, Fe‑S subunit
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (126/148); minor base C (0/1); total (126/149)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 278 total aligned reads displayed.

GGTGAAAAATATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTGCTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGAAACTGGCG  >  NC_000913/3017823‑3018018

ggtgAAAAATATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGgctact                                                                                                 <  1:3911206/101‑1 (MQ=255)

 gtgAAAAATATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTc                                                                                                >  1:3911637/1‑101 (MQ=255)

   gAAAAATATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCAc                                                                                              >  2:3911638/1‑101 (MQ=255)

   gAAAAATATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCAc                                                                                              <  2:3911077/101‑1 (MQ=255)

    aaaaaTATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCAcc                                                                                             >  1:3911643/1‑101 (MQ=255)

     aaaaTATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCAccc                                                                                            >  2:3911644/1‑101 (MQ=255)

      aaaTATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCg                                                                                           >  1:3911648/1‑101 (MQ=255)

      aaaTATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCg                                                                                           <  1:3911297/101‑1 (MQ=255)

      aaaTATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCg                                                                                           <  1:3911360/101‑1 (MQ=255)

          atCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGAttt                                                                                       <  2:3911334/101‑1 (MQ=255)

          atCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGAttt                                                                                       >  2:3911652/1‑101 (MQ=255)

            cATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTaa                                                                                     >  1:3911654/1‑101 (MQ=255)

              ttCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAAcc                                                                                   >  1:3911658/1‑101 (MQ=255)

              ttCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAAcc                                                                                   >  2:3911657/1‑101 (MQ=255)

               accTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCa                                                                                  >  1:3911662/2‑101 (MQ=255)

                ccTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCAt                                                                                 <  2:3911428/101‑1 (MQ=255)

                 cTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCAtt                                                                                >  2:3911664/1‑101 (MQ=255)

                 cTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCAtt                                                                                <  1:3911196/101‑1 (MQ=255)

                 cTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCAtt                                                                                <  1:3911258/101‑1 (MQ=255)

                 cTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCAtt                                                                                <  1:3911155/101‑1 (MQ=255)

                   cAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGa                                                                              <  2:3911571/101‑1 (MQ=255)

                    aGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAg                                                                             >  2:3911668/1‑101 (MQ=255)

                    aGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAg                                                                             >  1:3911670/1‑101 (MQ=255)

                     gTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGc                                                                            <  1:3911332/101‑1 (MQ=255)

                     gTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGc                                                                            <  2:3910982/101‑1 (MQ=255)

                       tCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAg                                                                          >  2:3911674/1‑101 (MQ=255)

                        ccGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGt                                                                         >  2:3911678/1‑101 (MQ=255)

                         cGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGtt                                                                        <  2:3911314/101‑1 (MQ=255)

                          gTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTa                                                                       <  1:3911391/101‑1 (MQ=255)

                          gTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTa                                                                       >  2:3911682/1‑101 (MQ=255)

                           tATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGAATGGCTACTCACCCGATTTAACCATTGAGCAGTTaa                                                                      >  1:3911683/1‑101 (MQ=255)

                            aTTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTaaa                                                                     >  2:3911685/1‑101 (MQ=255)

                              tCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTaaaaa                                                                   >  1:3911689/1‑101 (MQ=255)

                               ccTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGAtt                                                                                        >  1:3911694/1‑79 (MQ=255)

                               ccTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTaaaaaa                                                                  >  2:3911691/1‑101 (MQ=255)

                                cTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAAt                                                                 <  1:3911556/101‑1 (MQ=255)

                                 tGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATc                                                                <  2:3911463/101‑1 (MQ=255)

                                  gCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCa                                                               <  1:3911371/101‑1 (MQ=255)

                                   cAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAg                                                              <  2:3911290/101‑1 (MQ=255)

                                    agagTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAgg                                                             <  1:3911417/101‑1 (MQ=255)

                                    agagTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAgg                                                             <  2:3911444/101‑1 (MQ=255)

                                     gagTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAggg                                                            <  1:3911321/101‑1 (MQ=255)

                                      agTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGc                                                           <  2:3911308/101‑1 (MQ=255)

                                      agTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGc                                                           >  2:3911698/1‑101 (MQ=255)

                                        ttAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCtt                                                         <  2:3911274/101‑1 (MQ=255)

                                        ttAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCtt                                                         <  1:3911337/101‑1 (MQ=255)

                                         tAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTc                                                        >  1:3911700/1‑101 (MQ=255)

                                           aTTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCa                                                      <  2:3911319/101‑1 (MQ=255)

                                           aTTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCa                                                      <  2:3911389/101‑1 (MQ=255)

                                            ttCAGCACGATATCGATTTAGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCAc                                                     <  2:3911193/101‑1 (MQ=255)

                                             tCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACt                                                    <  2:3911513/101‑1 (MQ=255)

                                              cAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTa                                                   >  1:3911704/1‑101 (MQ=255)

                                               aGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTAt                                                  <  2:3911358/101‑1 (MQ=255)

                                                 cACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGt                                                <  1:3911283/101‑1 (MQ=255)

                                                 cACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGt                                                >  1:3911706/1‑101 (MQ=255)

                                                  aCGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTc                                                                                                >  2:3911711/1‑52 (MQ=255)

                                                  aCGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGtt                                               >  1:3911709/1‑101 (MQ=255)

                                                    gATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGc                                                                            <  1:3911241/70‑1 (MQ=255)

                                                    gATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCt                                             >  1:3911714/1‑101 (MQ=255)

                                                     atatCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTg                                            <  1:3911535/101‑1 (MQ=255)

                                                       atCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGAt                                          <  1:3911074/101‑1 (MQ=255)

                                                        tCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGAtt                                         <  2:3911450/101‑1 (MQ=255)

                                                        tCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGAtt                                         <  1:3911349/101‑1 (MQ=255)

                                                         cGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTg                                        >  2:3911717/1‑101 (MQ=255)

                                                          gATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGc                                       <  1:3911031/101‑1 (MQ=255)

                                                           aTTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGcc                                      <  2:3911262/101‑1 (MQ=255)

                                                           aTTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGcc                                      >  2:3911720/1‑101 (MQ=255)

                                                           aTTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGcc                                      >  1:3911721/1‑101 (MQ=255)

                                                                gttgCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCgg                                 >  1:3911725/1‑101 (MQ=255)

                                                                  tgCCGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCa                               >  2:3911727/1‑101 (MQ=255)

                                                                    ccGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACt                             <  1:3911483/101‑1 (MQ=255)

                                                                    ccGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACt                             <  2:3911378/101‑1 (MQ=255)

                                                                    ccGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACt                             <  1:3911465/101‑1 (MQ=255)

                                                                     cGCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTg                            >  2:3911731/1‑101 (MQ=255)

                                                                      gCTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGa                           >  2:3911734/1‑101 (MQ=255)

                                                                       cTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGAt                          <  1:3911214/101‑1 (MQ=255)

                                                                       cTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGAt                          >  1:3911735/1‑101 (MQ=255)

                                                                       cTCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGAt                          <  1:3911448/101‑1 (MQ=255)

                                                                        tCACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATa                         <  1:3911268/101‑1 (MQ=255)

                                                                         cACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATaa                        >  2:3911737/1‑101 (MQ=255)

                                                                         cACGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATaa                        <  1:3911676/101‑1 (MQ=255)

                                                                           cGGCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATaaaa                      <  2:3911292/101‑1 (MQ=255)

                                                                             gCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAAt                    >  2:3911741/1‑101 (MQ=255)

                                                                             gCGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAAt                    <  1:3911565/101‑1 (MQ=255)

                                                                              cGTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATa                   <  2:3911401/101‑1 (MQ=255)

                                                                               gTGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAg                  <  2:3911200/101‑1 (MQ=255)

                                                                                tGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGc                 >  1:3911746/1‑101 (MQ=255)

                                                                                tGAAATTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGc                 >  2:3911745/1‑101 (MQ=255)

                                                                                  aaaTTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCgg               <  2:3911408/101‑1 (MQ=255)

                                                                                  aaaTTTGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCgg               >  1:3911750/1‑101 (MQ=255)

                                                                                     tttGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGgtgt            <  2:3911243/101‑1 (MQ=255)

                                                                                       tGAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGa          >  1:3911753/1‑101 (MQ=255)

                                                                                        gAGTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGaa         <  1:3911398/101‑1 (MQ=255)

                                                                                          gTATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGAAAc       >  2:3911757/1‑101 (MQ=255)

                                                                                           tATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGAAACt      >  1:3911759/1‑101 (MQ=255)

                                                                                           tATGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGAAACt      <  1:3911356/101‑1 (MQ=255)

                                                                                            aTGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGAAACTg     <  2:3911420/101‑1 (MQ=255)

                                                                                            aTGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGAAACTg     <  2:3911300/101‑1 (MQ=255)

                                                                                            aTGGCTACTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGAAACTg     >  1:3911762/1‑101 (MQ=255)

                                                                                               gctactCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGAAACTggcg  <  1:3911603/101‑1 (MQ=255)

GGTGAAAAATATCATTCCTCAGTTCCGTATTCCTGCAGAGTTAATTCAGCACGATATCGATTTTGTTGCCGCTCACGGCGTGAAATTTGAGTATGGCTGCTCACCCGATTTAACCATTGAGCAGTTAAAAAATCAGGGCTTCCACTATGTTCTGATTGCCACCGGCACTGATAAAAATAGCGGTGTGAAACTGGCG  >  NC_000913/3017823‑3018018

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