BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,379,402:1	+A	58.6%	coding (596/1128 nt)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,402	1	.	A	58.6%	285.4 / 222.9	133	coding (596/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (19/36); new base A (45/33); total (64/69)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

GTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTAATACC  >  NC_000913/3379303‑3379499

 tttCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAAACAACACCAGCGTCCACGTTCATTACATCACAATGCTGtta                                                                                                   >  1:834939/1‑99 (MQ=255)

  ttCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGttat                                                                                                  <  2:2376967/100‑3 (MQ=255)

    cATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGttataa                                                                                                >  1:2302122/1‑96 (MQ=255)

    cATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAAt                                                                                               <  1:1071184/100‑1 (MQ=255)

     aTCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGttataat                                                                                               >  1:1808634/1‑95 (MQ=255)

     aTCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATg                                                                                              <  2:1874169/100‑1 (MQ=255)

       cGTGAAGTGGCGAAAGCCACAGATGCGCCTGGGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAAt                                                                                               <  1:2721388/97‑1 (MQ=255)

       cGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATgg                                                                                             <  2:1469449/100‑1 (MQ=255)

          gAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCAt                                                                                          >  2:3292009/1‑100 (MQ=255)

          gAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATc                                                                                         <  2:2264301/100‑1 (MQ=255)

          gAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATc                                                                                         <  1:3126137/100‑1 (MQ=255)

           aaGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATc                                                                                         <  1:3124809/99‑1 (MQ=255)

           aaGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCg                                                                                        <  2:735357/100‑1 (MQ=255)

            aGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGa                                                                                       >  1:1523736/1‑100 (MQ=255)

               ggCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTg                                                                                    >  1:3420428/1‑100 (MQ=255)

                 cGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgc                                                                                   <  2:1263004/100‑1 (MQ=255)

                 cGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgc                                                                                   <  2:650714/100‑1 (MQ=255)

                 cGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTgca                                                                                  <  2:1536144/100‑1 (MQ=255)

                  gAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgca                                                                                  >  2:1562170/1‑100 (MQ=255)

                  gAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTgcag                                                                                 <  1:1716883/100‑1 (MQ=255)

                   aaaGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgcag                                                                                 <  1:1928692/100‑1 (MQ=255)

                     aGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTgcaggca                                                                              <  2:2888278/100‑1 (MQ=255)

                     aGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTgcaggca                                                                              >  2:2153462/1‑100 (MQ=255)

                      gCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgcaggca                                                                              >  2:1775112/1‑100 (MQ=255)

                      gCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgcaggca                                                                              >  1:367316/1‑100 (MQ=255)

                      gCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTgcaggcag                                                                             >  2:1111682/1‑100 (MQ=255)

                      gCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTgcaggcag                                                                             >  2:3042571/1‑100 (MQ=255)

                       ccACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgcaggcag                                                                             >  2:2540585/1‑100 (MQ=255)

                        cacaGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTgcaggcag                                                                             >  1:1705604/1‑99 (MQ=255)

                           aGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGaaaa                                                                         <  1:3312564/100‑1 (MQ=255)

                            gATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAAc                                                                        >  1:1186869/1‑100 (MQ=255)

                            gATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAAcg                                                                       <  1:1723051/100‑1 (MQ=255)

                             aTGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAAcgc                                                                      <  2:1159334/100‑1 (MQ=255)

                               gcgcCTGAGTGAGTGTACGCAGACGATAATCAAGACCAGCGTCCACGTTCATTACATCACAATGCTGTTAGAATGGCATCGATTGCAGGCAGAAAAcgcg                                                                     >  1:2055504/1‑100 (MQ=255)

                               gcgcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAAcgcg                                                                     >  2:2236518/1‑100 (MQ=255)

                               gcgcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAAcgcg                                                                     <  2:3182017/100‑1 (MQ=255)

                               gcgcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAAcgcg                                                                     <  2:3354994/100‑1 (MQ=255)

                                 gcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCa                                                                   >  1:1786073/1‑100 (MQ=255)

                                      agtgagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGc                                                             >  2:1349816/1‑100 (MQ=255)

                                       gtgagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGAAAGAAAAAAGAAAACGCGCACGTTGc                                                             <  1:1149521/100‑1 (MQ=255)

                                        tgagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCa                                                            <  2:1433491/100‑1 (MQ=255)

                                        tgagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCa                                                            <  2:2738339/100‑1 (MQ=255)

                                         gagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAg                                                           >  2:842840/1‑100 (MQ=255)

                                         gagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAgg                                                          <  2:867237/100‑1 (MQ=255)

                                         gagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAgg                                                          <  1:3368046/100‑1 (MQ=255)

                                          agtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGc                                                         <  2:2024174/100‑1 (MQ=255)

                                          agtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACATTGCAGGc                                                         >  1:1144774/1‑100 (MQ=255)

                                           gtgtACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGc                                                         >  2:3283253/1‑100 (MQ=255)

                                              tACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCAt                                                      >  2:1645531/1‑99 (MQ=255)

                                               aCGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCAttt                                                    >  1:3000993/1‑100 (MQ=255)

                                                 gCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTc                                                   >  2:820468/1‑100 (MQ=255)

                                                    gACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGAt                                                <  2:82511/100‑1 (MQ=255)

                                                     aCGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGTAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATan                                              <  1:3038779/100‑2 (MQ=255)

                                                     aCGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAGAACGCGCACGTTGCAGGCCATTTCGATaa                                              <  1:2656733/100‑1 (MQ=255)

                                                     aCGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATaa                                              <  1:2656794/100‑1 (MQ=255)

                                                      cGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATaa                                              >  1:3419497/1‑100 (MQ=255)

                                                      cGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATaa                                              >  2:1809778/1‑100 (MQ=255)

                                                      cGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATaaa                                             >  1:874832/1‑100 (MQ=255)

                                                      cGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATaaa                                             >  1:876691/1‑100 (MQ=255)

                                                      cGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATaaa                                             >  1:3023410/1‑100 (MQ=255)

                                                       gATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAg                                            <  2:2499319/100‑1 (MQ=255)

                                                         tAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGt                                           <  2:672682/100‑1 (MQ=255)

                                                         tAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGtt                                          <  2:2621704/100‑1 (MQ=255)

                                                           aTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGt                                           <  1:2952079/98‑1 (MQ=255)

                                                             cAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCAt                                       >  1:1040138/1‑100 (MQ=255)

                                                             cAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATc                                      <  1:635682/100‑1 (MQ=255)

                                                               acacCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATcc                                     >  1:2381917/1‑100 (MQ=255)

                                                                cacCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCTTCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCAGTTCGATAAAGTTCATc                                      >  1:2938932/1‑98 (MQ=255)

                                                                 acCAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCcgg                                   >  2:3407472/1‑100 (MQ=255)

                                                                 acCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCcggc                                  <  1:1349810/100‑1 (MQ=255)

                                                                   cAGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCcggcg                                 >  1:1453782/1‑100 (MQ=255)

                                                                    aGCGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAATGTTCATCcggcgg                                <  2:2302134/100‑1 (MQ=255)

                                                                    aGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGa                               <  1:333137/100‑1 (MQ=255)

                                                                      cGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGaa                              >  2:2650371/1‑100 (MQ=255)

                                                                      cGTCCACGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGaa                              >  2:2336525/1‑100 (MQ=255)

                                                                       gTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTTCAGGCCATTTCGATAAAGATCATCCGGCGGAat                             >  2:25419/1‑99 (MQ=255)

                                                                        tCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAatat                           >  1:1655282/1‑100 (MQ=255)

                                                                           aCGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATAttt                         <  1:444712/100‑1 (MQ=255)

                                                                           aCGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTg                        <  1:1120047/100‑1 (MQ=255)

                                                                           aCGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTg                        <  2:475488/100‑1 (MQ=255)

                                                                            cGTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTg                        <  2:2028065/100‑1 (MQ=255)

                                                                            cGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGa                       <  1:1645524/100‑1 (MQ=255)

                                                                             gTTCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGa                       <  2:2731211/100‑1 (MQ=255)

                                                                             gTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGAc                      <  1:1615451/100‑1 (MQ=255)

                                                                               tCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACg                     >  2:2148973/1‑100 (MQ=255)

                                                                               tCATTACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACg                     <  1:393103/100‑1 (MQ=255)

                                                                                 aTTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTACAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCg                  >  2:1257367/1‑100 (MQ=255)

                                                                                  ttACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCCTTTCGCTAAAGTTCATCCCGCGGAATATTTGACGTc                   >  1:3123964/1‑99 (MQ=255)

                                                                                  ttACATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTc                   <  2:3350213/99‑1 (MQ=255)

                                                                                    aCATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCt                >  2:1141858/1‑100 (MQ=255)

                                                                                     cATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTa               >  2:277275/1‑100 (MQ=255)

                                                                                     cATCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTa               >  1:1619815/1‑100 (MQ=255)

                                                                                       tCACAATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTAcc             >  2:2022683/1‑100 (MQ=255)

                                                                                       tCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCa            <  1:1034942/100‑1 (MQ=255)

                                                                                         acaATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAg           >  1:1248816/1‑100 (MQ=255)

                                                                                          caATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAgg          >  2:46561/1‑100 (MQ=255)

                                                                                          caATGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATGTGACGTCGCTACCAgg          >  1:2230115/1‑100 (MQ=255)

                                                                                             tGCTGTTATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTa       <  2:1637190/100‑1 (MQ=255)

                                                                                                 gTT‑TAATGGCATCGATTGCAGGCAGAAAACGGGCACGTTGCAGGCCATTTCGATAAAGGTCATCCGGCGGAATATTTGACGTCGCTACAAGGGTAATAcc  >  1:3043543/1‑100 (MQ=255)

                                                                                                 gTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTAATCCGGCGGAATATTTGACGTCGCTACCAGGGTAATAcc  >  1:3043029/1‑100 (MQ=255)

GTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTT‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTAATACC  >  NC_000913/3379303‑3379499

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