BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,430,330	0	G	A	35.2%	397.7 / inf	324	A170V (GCG→GTG)	aroE	shikimate dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (94/116); new base A (63/51); total (157/167)
Fisher's exact test for biased strand distribution p-value = 8.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 329 total aligned reads displayed.

GAGCGATGACGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGG  >  NC_000913/3430230‑3430430

gaGCGATGACGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACg                                                                                                      >  2:6612250/1‑101 (MQ=255)

  gCGATGACGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGcc                                                                                                    >  2:6612253/1‑101 (MQ=255)

  gCGATGACGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACAcc                                                                                                    >  1:6612254/1‑101 (MQ=255)

    gATGACGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTg                                                                                                  <  1:6611962/101‑1 (MQ=255)

     aTGACGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGa                                                                                                 >  2:6612260/1‑101 (MQ=255)

     aTGACGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGa                                                                                                 >  1:6612258/1‑101 (MQ=255)

       gACGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAAt                                                                                               >  2:6612265/1‑101 (MQ=255)

       gACGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAAt                                                                                               <  1:6611783/101‑1 (MQ=255)

        aCGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATa                                                                                              <  1:6611874/101‑1 (MQ=255)

        aCGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATa                                                                                              >  1:6612266/1‑101 (MQ=255)

              tCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCa                                                                                        <  1:6611622/101‑1 (MQ=255)

               cGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAg                                                                                       <  1:6611774/101‑1 (MQ=255)

                gCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAgt                                                                                      <  2:6611913/101‑1 (MQ=255)

                gCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAgt                                                                                      >  2:6612272/1‑101 (MQ=255)

                  cGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAgtgt                                                                                    >  1:6612275/1‑101 (MQ=255)

                  cGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAgtgt                                                                                    >  2:6612277/1‑101 (MQ=255)

                    gAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTgc                                                                                  >  2:6612281/1‑101 (MQ=255)

                     aaTATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTgcg                                                                                 <  2:6611346/101‑1 (MQ=255)

                        atCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCaa                                                                              >  1:6612283/1‑101 (MQ=255)

                         tCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCaaa                                                                             <  2:6611649/101‑1 (MQ=255)

                           accacTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACa                                                                           >  1:6612287/1‑101 (MQ=255)

                            ccacTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGtgtg                                                                                   <  2:6611922/92‑1 (MQ=255)

                             cacTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACaa                                                                          <  1:6611841/100‑1 (MQ=255)

                             cacTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAAt                                                                         >  2:6612291/1‑101 (MQ=255)

                             cacTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAAt                                                                         <  1:6612153/101‑1 (MQ=255)

                             cacTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAAt                                                                         <  2:6611969/101‑1 (MQ=255)

                              acTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAAtt                                                                        >  1:6612293/1‑101 (MQ=255)

                               cTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAAttt                                                                       <  2:6611803/101‑1 (MQ=255)

                               cTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAAttt                                                                       <  1:6611754/101‑1 (MQ=255)

                                tGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTa                                                                      <  1:6611677/101‑1 (MQ=255)

                                   tGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGtgtg                                                                                   >  1:6612297/1‑85 (MQ=255)

                                    gCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCa                                                                  <  1:6611948/101‑1 (MQ=255)

                                    gCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCa                                                                  <  1:6611810/101‑1 (MQ=255)

                                      cACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAAc                                                                >  2:6612302/1‑101 (MQ=255)

                                        cTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAActc                                                              >  1:6612305/1‑101 (MQ=255)

                                            aTGTTGCATTAATAATGAGATCAAACTCCTGACCTTCCAGTTCGTCCATACTCAACGCCTg                                                                                                  >  2:6612308/1‑61 (MQ=255)

                                            aTGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTcc                                                          <  1:6611703/101‑1 (MQ=255)

                                             tGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCcg                                                         <  1:6612043/101‑1 (MQ=255)

                                              gTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCcgc                                                        <  2:6611761/101‑1 (MQ=255)

                                               ttGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCcgcg                                                       >  1:6612311/1‑101 (MQ=255)

                                                tGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCcgcgc                                                      <  1:6611821/101‑1 (MQ=255)

                                                tGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCcgcgc                                                      <  1:6611867/101‑1 (MQ=255)

                                                tGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCcgcgc                                                      >  1:6612314/1‑101 (MQ=255)

                                                   aTTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCggg                                                   >  2:6612316/1‑101 (MQ=255)

                                                   aTTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCggg                                                   <  2:6612026/101‑1 (MQ=255)

                                                    ttAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGa                                                  >  1:6612322/1‑101 (MQ=255)

                                                    ttAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGa                                                  >  2:6612321/1‑101 (MQ=255)

                                                     taataaTGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGAt                                                 >  1:6612327/1‑101 (MQ=255)

                                                     taataaTGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGAt                                                 <  2:6611773/101‑1 (MQ=255)

                                                     taataaTGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGAt                                                 <  2:6611827/101‑1 (MQ=255)

                                                     taataaTGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGAt                                                 <  1:6611381/101‑1 (MQ=255)

                                                      aataatGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATa                                                >  1:6612331/1‑101 (MQ=255)

                                                      aataatGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATa                                                <  1:6612000/101‑1 (MQ=255)

                                                         aatGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCg                                             <  2:6611747/101‑1 (MQ=255)

                                                         aatGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCg                                             >  2:6612334/1‑101 (MQ=255)

                                                          atGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGt                                            <  2:6611569/101‑1 (MQ=255)

                                                          atGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGt                                            >  1:6612335/1‑101 (MQ=255)

                                                           tGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTc                                           <  2:6611901/101‑1 (MQ=255)

                                                             agaTCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCg                                         <  2:6611855/101‑1 (MQ=255)

                                                               aTCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGAt                                       >  1:6612340/1‑101 (MQ=255)

                                                               aTCAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGAt                                       >  2:6612341/1‑101 (MQ=255)

                                                                tCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGAtt                                      <  1:6611861/101‑1 (MQ=255)

                                                                 cAAACTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTa                                     <  1:6611920/101‑1 (MQ=255)

                                                                   aaCTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGt                                   >  2:6612348/1‑101 (MQ=255)

                                                                     cTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTa                                 >  2:6612352/1‑101 (MQ=255)

                                                                     cTCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTa                                 >  1:6612351/1‑101 (MQ=255)

                                                                      tCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTAt                                <  1:6611726/101‑1 (MQ=255)

                                                                      tCATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTAt                                <  2:6612053/101‑1 (MQ=255)

                                                                       cATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTAtt                               >  2:6612356/1‑101 (MQ=255)

                                                                       cATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTAtt                               >  1:6612362/1‑101 (MQ=255)

                                                                       cATGACCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTAtt                               >  2:6612359/1‑101 (MQ=255)

                                                                        aTGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTAt                                >  2:6612366/1‑99 (MQ=255)

                                                                         tGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGt                             <  1:6611881/101‑1 (MQ=255)

                                                                          gACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTc                            <  2:6611787/101‑1 (MQ=255)

                                                                           aCCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCa                           >  2:6612370/1‑101 (MQ=255)

                                                                           aCCTTCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCa                           <  2:6611945/101‑1 (MQ=255)

                                                                            ccTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCAc                          <  2:6611680/101‑1 (MQ=255)

                                                                             cTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCAcc                         >  1:6612372/1‑101 (MQ=255)

                                                                              ttCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCg                        <  2:6612110/101‑1 (MQ=255)

                                                                              ttCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCg                        >  2:6612374/1‑101 (MQ=255)

                                                                               tCCAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCAAACTCTTCCGCGCGGGa                                                  >  1:6612379/1‑74 (MQ=255)

                                                                                 cAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGcac                     <  2:6612063/101‑1 (MQ=255)

                                                                                 cAGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGcac                     <  1:6611899/101‑1 (MQ=255)

                                                                                  aGTTCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGcaca                    >  1:6612383/1‑101 (MQ=255)

                                                                                   gTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAg                   >  2:6612385/1‑101 (MQ=255)

                                                                                    ttCGTCCATACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGt                  >  1:6612388/1‑101 (MQ=255)

                                                                                      cGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTcc                <  2:6611796/101‑1 (MQ=255)

                                                                                      cGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTcc                <  2:6611738/101‑1 (MQ=255)

                                                                                      cGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTcc                >  2:6612391/1‑101 (MQ=255)

                                                                                        tCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAg              >  1:6612395/1‑101 (MQ=255)

                                                                                          cATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAggg            <  1:6611848/101‑1 (MQ=255)

                                                                                          cATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAggg            <  2:6611892/101‑1 (MQ=255)

                                                                                           aTACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGa           >  1:6612401/1‑101 (MQ=255)

                                                                                           aTACTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGa           >  2:6612400/1‑101 (MQ=255)

                                                                                             aCTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGaaa         <  2:6611999/101‑1 (MQ=255)

                                                                                             aCTCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGaaa         >  1:6612406/1‑101 (MQ=255)

                                                                                              cTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAg        <  1:6612013/101‑1 (MQ=255)

                                                                                               tCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAgg       <  1:6611658/101‑1 (MQ=255)

                                                                                               tCAACACCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAgg       <  2:6611981/101‑1 (MQ=255)

                                                                                                   cGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTg   <  1:6611857/101‑1 (MQ=255)

GAGCGATGACGGGATCGCCGGAATATCACCACTGATGCCACTGGATGTTGCATTAATAATGAGATCAAACTCATGACCTTCCAGTTCGTCCATACTCAACGCCTGAATACTGCCAGTGTGCGCAAACAATTTAGCCAACTCTTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGG  >  NC_000913/3430230‑3430430

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: `‑`