BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,197,340	T→C	T65A (ACT→GCT)	ymfD ←	e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,197,340	0	T	C	97.9%	499.9 / ‑4.9	144	T65A (ACT→GCT)	ymfD	e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (68/73); total (70/74)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

GCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATAT  >  NC_000913/1197243‑1197435

 cATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAAc                                                                                              <  1:1566367/100‑1 (MQ=255)

  aTTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACt                                                                                             <  2:1730920/100‑1 (MQ=255)

  aTTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACt                                                                                             <  2:729326/100‑1 (MQ=255)

  aTTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACt                                                                                             <  2:729777/100‑1 (MQ=255)

    ttttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTtc                                                                                           <  1:2612762/100‑1 (MQ=255)

     tttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtca                                                                                          >  2:2295531/1‑100 (MQ=255)

     tttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtca                                                                                          >  1:488568/1‑100 (MQ=255)

     tttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtca                                                                                          <  1:2169037/100‑1 (MQ=255)

       tataATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtcatc                                                                                        <  2:2765910/100‑1 (MQ=255)

            tatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAtt                                                                                   >  2:1332401/1‑100 (MQ=255)

            tatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAtt                                                                                   <  2:2237817/100‑1 (MQ=255)

             atCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAttt                                                                                  >  2:284017/1‑100 (MQ=255)

             atCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAttt                                                                                  >  2:931115/1‑100 (MQ=255)

                gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                               <  2:2649060/100‑1 (MQ=255)

                gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                               >  1:3036521/1‑100 (MQ=255)

                gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                               >  1:295276/1‑100 (MQ=255)

                gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                               <  2:1442369/100‑1 (MQ=255)

                gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                               >  2:2007633/1‑100 (MQ=255)

                 tGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTg                                                                              >  2:2727789/1‑100 (MQ=255)

                 tGGGACATAGTCAATAATTTTAGTTTTAACTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCCAGTCTAGAAAAGCAACTTCATCAAATTTACTg                                                                              <  2:3299385/100‑1 (MQ=255)

                    gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGAtt                                                                           <  1:1161393/100‑1 (MQ=255)

                    gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGAtt                                                                           >  2:913718/1‑100 (MQ=255)

                      cATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCACTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTaa                                                                         >  2:563610/1‑100 (MQ=255)

                       aTAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAt                                                                        <  1:3021874/100‑1 (MQ=255)

                        tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                       <  1:2486706/100‑1 (MQ=255)

                        tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                       <  1:562002/100‑1 (MQ=255)

                         ggtatataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGTCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtc                                                                      <  2:3247907/95‑1 (MQ=255)

                          gTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtca                                                                     >  1:1385417/1‑100 (MQ=255)

                            cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                   >  2:410593/1‑100 (MQ=255)

                            cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                   >  2:31427/1‑100 (MQ=255)

                            cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                   <  2:64462/100‑1 (MQ=255)

                            cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                   >  1:1677733/1‑100 (MQ=255)

                             aataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                   <  1:1061008/99‑1 (MQ=255)

                              ataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAg                                                                 <  2:2410039/100‑1 (MQ=255)

                              ataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAg                                                                 <  2:1398955/100‑1 (MQ=255)

                              ataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAg                                                                 <  1:1089779/100‑1 (MQ=255)

                              ataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAg                                                                 <  2:1962953/100‑1 (MQ=255)

                                aatTTTAGTTTTAATTCCTCTAATAATTTGCGCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGaa                                                               >  2:1560811/1‑100 (MQ=255)

                                 atTTTAGTTTTAATTCCTCTAATATTTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAat                                                              <  2:179804/100‑1 (MQ=255)

                                 atTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAat                                                              <  1:337714/100‑1 (MQ=255)

                                 atTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAat                                                              <  2:3393609/100‑1 (MQ=255)

                                  ttttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAata                                                             <  2:295280/100‑1 (MQ=255)

                                  ttttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAata                                                             >  2:678358/1‑100 (MQ=255)

                                   tttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAatat                                                            >  2:506009/1‑100 (MQ=255)

                                    ttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATc                                                           <  2:1179971/100‑1 (MQ=255)

                                      aGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTa                                                         <  2:2742256/100‑1 (MQ=255)

                                       gTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTaa                                                        >  2:2415957/1‑100 (MQ=255)

                                        ttttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAAGTTCATCAACTTTTCTGATTAATTCATCAGAATATCTAAg                                                       >  1:2438381/1‑100 (MQ=255)

                                        ttttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAg                                                       <  2:1043607/100‑1 (MQ=255)

                                        ttttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAg                                                       <  1:2215257/100‑1 (MQ=255)

                                         tttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTATAAAAGCAATTTCATCAAATTTACTGATTAATTCATCAGAATATCTaa                                                        >  1:2001130/1‑98 (MQ=255)

                                           tAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGttt                                                    <  2:2093081/100‑1 (MQ=255)

                                              ttCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTcc                                                 >  2:324829/1‑100 (MQ=255)

                                              ttCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTcc                                                 <  2:1334936/100‑1 (MQ=255)

                                                ccTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTATAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGc                                               >  1:1316961/1‑100 (MQ=255)

                                                ccTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGc                                               >  2:585564/1‑100 (MQ=255)

                                                ccTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAA‑CAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCa                                              <  2:1105623/100‑1 (MQ=255)

                                                  tctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAg                                             >  2:1462251/1‑100 (MQ=255)

                                                   ctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAGATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGc                                            >  1:2436776/1‑100 (MQ=255)

                                                    taataaTTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGcc                                           <  1:290752/100‑1 (MQ=255)

                                                    taataaTTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGcc                                           >  2:3079860/1‑100 (MQ=255)

                                                    taataaTTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGcc                                           >  2:2470170/1‑100 (MQ=255)

                                                       taatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCaaa                                        <  2:1624286/100‑1 (MQ=255)

                                                        aatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCaaaa                                       >  1:1206171/1‑100 (MQ=255)

                                                           ttGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCa                                    <  2:1120520/100‑1 (MQ=255)

                                                           ttGCCCGCTTTCACCTTGCCTTTTCGAATCTAGAAAAGCAACTTCATCACATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAACCa                                    <  1:2021117/100‑1 (MQ=255)

                                                            tGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCaa                                   <  2:1385437/100‑1 (MQ=255)

                                                              ctctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAGGTTTTCCGCAGCCAAAATCAAgg                                 <  1:3332647/100‑1 (MQ=255)

                                                                  ctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGt                             >  1:2621931/1‑100 (MQ=255)

                                                                  ctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGt                             >  1:662004/1‑100 (MQ=255)

                                                                    ttCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTgg                           <  1:1371462/100‑1 (MQ=255)

                                                                     tCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGc                          <  2:1419819/100‑1 (MQ=255)

                                                                      cAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGcc                         <  1:3181591/100‑1 (MQ=255)

                                                                         gTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtt                      >  2:3402814/1‑100 (MQ=255)

                                                                         gTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtt                      >  2:2077296/1‑100 (MQ=255)

                                                                             ccTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAttttt                   >  1:569310/1‑99 (MQ=255)

                                                                             ccTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtc                  <  2:1574194/100‑1 (MQ=255)

                                                                             ccTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtc                  <  2:2184446/100‑1 (MQ=255)

                                                                             ccTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtc                  >  2:1618420/1‑100 (MQ=255)

                                                                              cTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtct                 <  2:1677756/100‑1 (MQ=255)

                                                                               ttttCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtctc                <  1:2295499/100‑1 (MQ=255)

                                                                               ttttCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtctc                <  2:1903252/100‑1 (MQ=255)

                                                                                tttCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCa               <  1:2238306/100‑1 (MQ=255)

                                                                                   cGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATg            >  2:1864151/1‑100 (MQ=255)

                                                                                    gAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTGAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGc           >  1:3216179/1‑100 (MQ=255)

                                                                                    gAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGc           <  1:43087/100‑1 (MQ=255)

                                                                                      gTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCtt         >  2:415924/1‑100 (MQ=255)

                                                                                      gTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCtt         >  2:1995420/1‑100 (MQ=255)

                                                                                      gTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCtt         >  2:957427/1‑100 (MQ=255)

                                                                                       tCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTc        <  1:913706/100‑1 (MQ=255)

                                                                                        cTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt       <  1:2893428/100‑1 (MQ=255)

                                                                                        cTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt       >  1:2319125/1‑100 (MQ=255)

                                                                                        cTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt       >  1:1517546/1‑100 (MQ=255)

                                                                                        cTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt       >  2:919254/1‑100 (MQ=255)

                                                                                         tAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTa      >  1:2481167/1‑100 (MQ=255)

                                                                                          aGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTaa     <  2:2934627/100‑1 (MQ=255)

                                                                                           gAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat    <  2:2766001/100‑1 (MQ=255)

                                                                                           gAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat    <  2:1016778/100‑1 (MQ=255)

                                                                                            aaaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAata   <  1:513760/100‑1 (MQ=255)

                                                                                            aaaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAata   <  2:1931777/100‑1 (MQ=255)

GCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATAT  >  NC_000913/1197243‑1197435

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