BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,241,465	T→C	D397G (GAC→GGC)	malF ←	MalF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,241,465	0	T	C	99.6%	705.4 / ‑5.1	226	D397G (GAC→GGC)	malF	MalF
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (106/119); total (107/119)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

Only 100 of 237 total aligned reads displayed.

GTCAGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGTCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTGCGGGC  >  NC_000913/4241375‑4241553

gTCAGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGccc                                                                                        >  1:2260534/1‑93 (MQ=255)

gTCAGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGccc                                                                                        >  2:471204/1‑93 (MQ=255)

 tCAGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCg                                                                                       >  2:2157227/1‑93 (MQ=255)

   aGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGa                                                                                     >  1:1677475/1‑93 (MQ=255)

   aGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGa                                                                                     >  2:552822/1‑93 (MQ=255)

     cGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAAt                                                                                   <  1:765018/93‑1 (MQ=255)

     cGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAAt                                                                                   <  2:665994/93‑1 (MQ=255)

       gTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCg                                                                                 <  2:799777/93‑1 (MQ=255)

        tttAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGc                                                                                <  1:4192556/93‑1 (MQ=255)

         ttAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCt                                                                               <  1:4123059/93‑1 (MQ=255)

          tAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCtt                                                                              <  2:3200742/93‑1 (MQ=255)

          tAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCtt                                                                              >  1:3341565/1‑93 (MQ=255)

          tAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCtt                                                                              >  2:4098681/1‑93 (MQ=255)

           aaTCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCttt                                                                             >  2:3862175/1‑93 (MQ=255)

            aTCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTc                                                                            <  1:1746375/93‑1 (MQ=255)

             tCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCa                                                                           <  2:2652395/93‑1 (MQ=255)

              cagctgCGGCAGCGTAATCTTAAAGAAGTTCTGGACCGGCCCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAg                                                                          <  1:264341/93‑1 (MQ=255)

               agcagcGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGc                                                                         <  1:2902427/93‑1 (MQ=255)

                gcagcggcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCa                                                                        >  2:348797/1‑93 (MQ=255)

                 cagcggcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCaa                                                                       <  1:2890698/93‑1 (MQ=255)

                   gcggcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGc                                                                     <  1:441657/93‑1 (MQ=255)

                    cggcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAgtc                                                                    >  2:2886255/1‑91 (MQ=255)

                    cggcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGcc                                                                    <  1:3696185/93‑1 (MQ=255)

                       cagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGCGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCAt                                                                 >  1:79066/1‑93 (MQ=255)

                       cagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCAt                                                                 <  1:692365/93‑1 (MQ=255)

                        agcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATg                                                                >  2:1591871/1‑93 (MQ=255)

                           gTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGTCCATGCag                                                             >  2:1607653/1‑93 (MQ=255)

                            tAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCaga                                                            >  2:2974108/1‑93 (MQ=255)

                            tAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCaga                                                            <  1:1628196/93‑1 (MQ=255)

                            tAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCaga                                                            >  2:1245018/1‑93 (MQ=255)

                            tAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCaga                                                            >  1:4520660/1‑93 (MQ=255)

                              aTCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAgg                                                          >  2:3167662/1‑93 (MQ=255)

                               tCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGa                                                         <  1:1810165/93‑1 (MQ=255)

                               tCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGa                                                         >  1:1185761/1‑93 (MQ=255)

                                  tAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGatca                                                      <  1:1013656/93‑1 (MQ=255)

                                    aagaagTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGatcatc                                                    >  1:849683/1‑93 (MQ=255)

                                    aagaagTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGatcatc                                                    <  1:1119732/93‑1 (MQ=255)

                                    aagaagTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGatcatc                                                    <  2:4422957/93‑1 (MQ=255)

                                     agaagTTCTGTAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGTAATCGCTTTCAGCAAGCCTATGCAGAGGAtcatca                                                   <  2:1299660/93‑1 (MQ=255)

                                         gTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTa                                               <  2:757484/93‑1 (MQ=255)

                                         gTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTa                                               <  2:635637/93‑1 (MQ=255)

                                          ttCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTAc                                              >  2:821635/1‑93 (MQ=255)

                                          ttCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTAc                                              <  2:71511/93‑1 (MQ=255)

                                             tGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGGTGATCATGTACGGa                                           >  2:91346/1‑93 (MQ=255)

                                             tGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGa                                           <  2:1098667/93‑1 (MQ=255)

                                             tGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGa                                           <  2:326795/93‑1 (MQ=255)

                                               gAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATa                                         >  2:2025354/1‑93 (MQ=255)

                                                aaCGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATaa                                        <  2:1481611/93‑1 (MQ=255)

                                                 aCGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAAc                                       >  2:2468923/1‑93 (MQ=255)

                                                 aCGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAAc                                       >  2:3613698/1‑93 (MQ=255)

                                                   ggACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAAccc                                     >  1:194321/1‑93 (MQ=255)

                                                    gACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACcca                                    <  1:3236900/93‑1 (MQ=255)

                                                    gACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACcca                                    >  1:1865865/1‑93 (MQ=255)

                                                       cTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAATCCCATGCAGAGGATCATCATGTACGGATAACccagcc                                 <  1:2099224/93‑1 (MQ=255)

                                                       cTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagcc                                 >  1:1385806/1‑93 (MQ=255)

                                                        tGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagcca                                <  2:1716009/93‑1 (MQ=255)

                                                         gcgcCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagccag                               <  2:2803519/93‑1 (MQ=255)

                                                         gcgcCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagccag                               <  1:2321157/93‑1 (MQ=255)

                                                           gcCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGt                             <  2:1349084/93‑1 (MQ=255)

                                                           gcCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGt                             >  2:3993393/1‑93 (MQ=255)

                                                            ccatccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTa                            >  1:4634548/1‑93 (MQ=255)

                                                              atccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTAtt                          <  2:1791303/93‑1 (MQ=255)

                                                              atccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTAtt                          >  1:2924889/1‑93 (MQ=255)

                                                              atccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTAtt                          >  2:4325299/1‑93 (MQ=255)

                                                               tccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTg                         >  1:2554540/1‑93 (MQ=255)

                                                               tccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTg                         >  2:2623456/1‑93 (MQ=255)

                                                                ccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGa                        >  2:4162463/1‑93 (MQ=255)

                                                                  atTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACg                      <  1:2485780/93‑1 (MQ=255)

                                                                  atTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACg                      <  2:3405179/93‑1 (MQ=255)

                                                                  atTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACg                      >  2:2842684/1‑93 (MQ=255)

                                                                   ttGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGAAGa                     >  1:3863069/1‑93 (MQ=255)

                                                                    tGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGAt                    <  1:4350577/93‑1 (MQ=255)

                                                                     gCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATa                   >  1:30907/1‑93 (MQ=255)

                                                                      cTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATaa                  <  1:1535442/93‑1 (MQ=255)

                                                                        gAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAAtt                <  1:3966387/93‑1 (MQ=255)

                                                                        gAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAAtt                >  2:3741758/1‑93 (MQ=255)

                                                                         aGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTa               <  1:283329/93‑1 (MQ=255)

                                                                         aGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTa               <  1:1353834/93‑1 (MQ=255)

                                                                          ggCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAg              >  2:3458741/1‑93 (MQ=255)

                                                                          ggCTTCATACAAATCGCCCGGAATCCCTTTCAGCAAGCCCATTCAGAGGATCATCATGTACCGTTAACCCATCCAGGGATTTACGATAATTAg              >  2:3710853/1‑93 (MQ=255)

                                                                           gCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGc             <  2:2740824/93‑1 (MQ=255)

                                                                           gCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGc             >  2:4602462/1‑93 (MQ=255)

                                                                            cTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCa            <  1:4542341/93‑1 (MQ=255)

                                                                            cTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCa            <  2:499683/93‑1 (MQ=255)

                                                                             ttCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCAt           >  2:1547090/1‑93 (MQ=255)

                                                                             ttCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCAt           <  2:2413822/93‑1 (MQ=255)

                                                                              tCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATc          <  2:1345754/93‑1 (MQ=255)

                                                                              tCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATc          <  2:1599853/93‑1 (MQ=255)

                                                                               cATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCg         >  2:2122099/1‑93 (MQ=255)

                                                                                aTACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGt        >  1:543335/1‑93 (MQ=255)

                                                                                aTACAAATCGCCCGGAATCGCTTTCAGCAAGCCCACGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGt        >  2:4084721/1‑93 (MQ=255)

                                                                                 tACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTg       >  1:1276682/1‑93 (MQ=255)

                                                                                  aCAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTgc      >  1:2133818/1‑93 (MQ=255)

                                                                                  aCAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTgc      <  1:2209588/93‑1 (MQ=255)

                                                                                  aCAAAGCGCCCGGAATCGCGTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCAGCGTgc      <  2:3039473/93‑1 (MQ=255)

                                                                                   cAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTgcg     <  1:303589/93‑1 (MQ=255)

                                                                                    aaaTCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTgcgg    >  2:219476/1‑93 (MQ=255)

                                                                                    aaaTCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTgcgg    <  2:1888553/93‑1 (MQ=255)

                                                                                     aaTCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTgcggg   >  1:3457774/1‑93 (MQ=255)

                                                                                      aTCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGACAACCCAGCCAGGTATTGACGATAATTAGCATCGTgcgggg  >  1:2047344/1‑92 (MQ=255)

GTCAGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGTCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTGCGGGC  >  NC_000913/4241375‑4241553

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