BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,113,473	T→C	19.0%	L61L (TTG→CTG) ‡	uspD →	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,473	0	T	C	19.0%	380.1 / 84.4	185	L61L (TTG→CTG)	uspD	stress‑induced protein
Reads supporting (aligned to +/- strand): ref base T (79/70); new base C (18/17); total (98/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-02

Only 100 of 206 total aligned reads displayed.

AGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAA  >  NC_000913/4113374‑4113572

  aCATAATGACTCTCTCCTGACGTTAATTCATATTGAAGATTGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAtt                                                                                                    <  1:2639520/99‑1 (MQ=255)

  aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                   <  1:352001/100‑1 (MQ=255)

  aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                   <  2:3741756/100‑1 (MQ=255)

  aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                   <  2:1048351/100‑1 (MQ=255)

  aCATAATGACGCTCACCTGACGTTAATTCAGATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                   <  1:3692464/100‑1 (MQ=255)

   cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                  >  1:1129418/1‑100 (MQ=255)

    aTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgtt                                                                                                 <  1:1190922/100‑1 (MQ=255)

     tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                >  1:3483599/1‑100 (MQ=255)

     tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                >  2:4197795/1‑100 (MQ=255)

     tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                <  2:2592211/100‑1 (MQ=255)

        tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                             <  2:83308/100‑1 (MQ=255)

           cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                                                                                          >  2:3495389/1‑100 (MQ=255)

            gCTCACCTGACGTTAAGTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                         <  1:2420147/100‑1 (MQ=255)

              tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatc                                                                                       <  2:3888754/100‑3 (MQ=255)

              tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatc                                                                                       <  1:1032100/100‑3 (MQ=255)

               cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatct                                                                                      <  1:1272991/100‑4 (MQ=38)

                 ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                    <  2:906825/100‑6 (MQ=38)

                 ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                    <  1:3684702/100‑6 (MQ=38)

                  cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                   <  1:1592071/100‑7 (MQ=38)

                  cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                   <  2:1352979/100‑7 (MQ=38)

                   tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctg                                                                                  <  2:4307651/100‑8 (MQ=37)

                   tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctg                                                                                  <  1:3268366/100‑8 (MQ=37)

                    gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctgc                                                                                 <  1:3306280/100‑9 (MQ=37)

                    aactttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGa                                                                                  <  1:4218126/95‑1 (MQ=255)

                       gTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAAc                                                                              >  1:3352020/1‑100 (MQ=255)

                             tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTg                                                                        <  2:2324763/100‑1 (MQ=255)

                               atatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGta                                                                      >  2:1526874/1‑100 (MQ=255)

                                tatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtat                                                                     >  2:3697210/1‑100 (MQ=255)

                                tatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtat                                                                     >  2:1315716/1‑100 (MQ=255)

                                 atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAATATATTCTTGAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                    <  2:3103542/100‑1 (MQ=255)

                                 atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                    <  1:963255/100‑1 (MQ=255)

                                   tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                  >  1:1160554/1‑100 (MQ=255)

                                    gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                  >  2:337114/1‑99 (MQ=255)

                                    gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                 <  2:2243306/100‑1 (MQ=255)

                                        atgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                             >  2:2515550/1‑100 (MQ=255)

                                           gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                          >  2:86194/1‑100 (MQ=255)

                                           gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                          >  2:4086224/1‑100 (MQ=255)

                                            cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaa                                                         >  1:880920/1‑100 (MQ=255)

                                              tAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaaa                                                       >  2:1181201/1‑100 (MQ=255)

                                               aaGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAat                                                      >  1:2108090/1‑100 (MQ=255)

                                               aaGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAat                                                      >  2:1205171/1‑100 (MQ=255)

                                                aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                     >  1:1986848/1‑100 (MQ=255)

                                                aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                     >  2:2706600/1‑100 (MQ=255)

                                                aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                     >  2:3177165/1‑100 (MQ=255)

                                                aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                     >  2:3473514/1‑100 (MQ=255)

                                                aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                     >  1:3160619/1‑100 (MQ=255)

                                                 gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                    >  2:2122224/1‑100 (MQ=255)

                                                  cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATATCAAGCTGTATAAACTGACGAAAAATAtt                                                   >  1:157829/1‑100 (MQ=255)

                                                   gggTTGTATCCGGGTATCTACTTCTCTTTAACAGAATATATTCTTTAATTGTTGAAGAATATGTCGTATAACAAGCTGTATAAACTGACGAAAAATAtt                                                   <  1:2348279/97‑1 (MQ=255)

                                                   gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                  <  2:3501374/100‑1 (MQ=255)

                                                   gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                  >  1:3393236/1‑100 (MQ=255)

                                                   gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                  <  2:101765/100‑1 (MQ=255)

                                                   gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGGTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                  <  1:2692061/100‑1 (MQ=255)

                                                    aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                 <  1:586625/100‑1 (MQ=255)

                                                        gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                             >  2:2678979/1‑100 (MQ=255)

                                                        gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                             >  2:3584745/1‑100 (MQ=255)

                                                          aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATTGc                                            >  1:4103115/1‑99 (MQ=255)

                                                          aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                           >  2:2946651/1‑100 (MQ=255)

                                                           tCCGGGTATCTACTTCCCTGCAACAGCAGATATTCTGCAATTGTTGAAGACTAAGTCGGATACCAAGCTGTATAAACTCACGAAAAATATTCAATGGcc                                           >  1:426622/1‑99 (MQ=255)

                                                           tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGACCg                                          >  2:790205/1‑100 (MQ=255)

                                                            ccGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGa                                         >  2:707888/1‑100 (MQ=255)

                                                               ggTATCTACTTCCCTGCAACAGAAGATATTCGTCAATTGTTGCAGCATAAGTCGGATAAAAAGCTGTATACACTGACGAAAAATATTCAATGGCCGAAGa                                      >  1:806885/1‑100 (MQ=255)

                                                                  aTCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaa                                   <  1:1407648/100‑1 (MQ=255)

                                                                   tCTACTTCCCTGAAACATAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTCTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaa                                  <  1:2959411/100‑1 (MQ=255)

                                                                    cTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATACACTGACGAAAAATATTCAATGGCCGAAGACaaaa                                 >  2:434600/1‑100 (MQ=255)

                                                                     tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                >  1:4159388/1‑100 (MQ=255)

                                                                     tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                <  1:2761935/100‑1 (MQ=255)

                                                                       cTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                              >  2:3762574/1‑100 (MQ=255)

                                                                       cTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                              >  2:1670825/1‑100 (MQ=255)

                                                                        ttCCCTTCAACAGAAGATAGTCTTCAATTGTTGAAGGATAAGTCGGATAATAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaaa                                 <  1:3314750/96‑1 (MQ=255)

                                                                         tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                             <  1:4388731/99‑1 (MQ=255)

                                                                          cccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                           <  1:1786671/100‑1 (MQ=255)

                                                                          cagtGCAACAGAAGATATCCTGCAATGGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                            <  2:4215035/96‑1 (MQ=255)

                                                                           ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGATGACAAAACTGCGt                           <  1:2727959/99‑1 (MQ=255)

                                                                            cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                         >  1:2070502/1‑100 (MQ=255)

                                                                            cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                         >  2:2415709/1‑100 (MQ=255)

                                                                             tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                        >  2:1481987/1‑100 (MQ=255)

                                                                             tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                        >  2:4325597/1‑100 (MQ=255)

                                                                              gCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTg                       <  1:1373000/100‑1 (MQ=255)

                                                                              gCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGAGAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                         <  1:1322781/98‑1 (MQ=255)

                                                                                aaCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGaa                     >  2:2007008/1‑100 (MQ=255)

                                                                                 aCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                    >  2:1822235/1‑100 (MQ=255)

                                                                                  cAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcg                   <  1:1171832/100‑1 (MQ=255)

                                                                                    gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                 >  1:1863105/1‑100 (MQ=255)

                                                                                    cgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                 >  1:3044919/6‑100 (MQ=38)

                                                                                    cgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                 >  2:1067124/6‑100 (MQ=38)

                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                >  2:1981649/5‑100 (MQ=38)

                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                >  2:3212166/5‑100 (MQ=38)

                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                >  2:2138063/5‑100 (MQ=38)

                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                >  1:3976413/5‑100 (MQ=38)

                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATACAATGGCCGAATACAAAACTGCGTATTGaa                     >  1:3995243/5‑95 (MQ=37)

                                                                                      acttatTCTTCAACAATTGAAGGCTAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga               >  1:2567306/4‑100 (MQ=255)

                                                                                       gaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag              <  1:2228589/100‑1 (MQ=255)

                                                                                              ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc       <  2:2699866/100‑1 (MQ=255)

                                                                                                cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg     <  2:3604804/100‑1 (MQ=255)

                                                                                                cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg     <  2:3182665/100‑1 (MQ=255)

                                                                                                cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg     <  2:3042297/100‑1 (MQ=255)

                                                                                                  aTTGTTGAAGAATAAGCCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa   >  2:606986/1‑100 (MQ=255)

                                                                                                   ttgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaaa  <  2:1632617/100‑1 (MQ=255)

                                                                                                   ttgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa   >  2:2665298/1‑99 (MQ=255)

AGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAA  >  NC_000913/4113374‑4113572

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