BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,638,193:1	(C)_6→7	100%	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	96.2%	350.4 / 5.4	104	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (1/3); new base C (53/47); total (54/50)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 121 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                   <  2:2277524/100‑1 (MQ=255)

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                                 >  2:1984924/1‑100 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc                                                                                                >  2:424616/1‑100 (MQ=255)

      ttccgatctCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑cc                                                                                          <  1:1899046/92‑1 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑cc                                                                                          >  1:54289/1‑100 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑cc                                                                                          >  1:209526/1‑100 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑cc                                                                                          >  2:1555552/1‑100 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCGTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑ccc                                                                                         <  1:204176/100‑2 (MQ=255)

         cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCAGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑ccctg                                                                                       >  1:1867273/1‑97 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑ccctgt                                                                                      >  2:284062/1‑96 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑ccctg                                                                                       >  1:2407974/1‑96 (MQ=255)

           aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑ccctgtt                                                                                     >  2:2368498/1‑95 (MQ=255)

           aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑ccctgtt                                                                                     <  1:1566389/100‑6 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTc                                                                                     >  2:2156073/1‑100 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGATGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTg                                                                                   <  1:1617108/100‑1 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTgg                                                                                  >  1:1649895/1‑100 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑c‑ccagatcggaa                                                                                 >  2:1899054/1‑91 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGTaa                                                                                >  2:2317714/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGaa                                                                                <  1:284060/100‑1 (MQ=255)

                   gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGaaaa                                                                              >  2:2769881/1‑100 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATg                                                                            <  1:2889399/100‑1 (MQ=255)

                       cTTGACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCACAGGAGACGACTATGTTGCAATCCCC‑C‑CCTGTTCTGGAAAATGAc                                                                          >  1:1113045/1‑100 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGAc                                                                          >  2:1723737/1‑100 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGAc                                                                          <  1:1988277/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGAc                                                                          <  2:2346854/100‑1 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTTTTCTGGaaaa                                                                              >  1:1784868/1‑93 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAg                                                                       >  2:1805594/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  2:1921894/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  2:1541380/100‑1 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:587790/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  1:198271/100‑1 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             <  1:1892134/100‑1 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:250900/1‑100 (MQ=255)

                                     aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  1:993604/1‑100 (MQ=255)

                                     aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  1:2350640/1‑100 (MQ=255)

                                      aaTTCTTCCCCCATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  2:1192406/100‑1 (MQ=255)

                                       tctcttcCCATATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑C‑CGTGCTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  1:1907876/98‑1 (MQ=255)

                                       aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCACAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          >  1:714604/1‑100 (MQ=255)

                                        ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  1:607382/100‑1 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       <  2:345214/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGTAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTACCCTGTTTGGTGCAGAATTg                                                    >  2:888908/1‑99 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:144983/100‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  1:1225472/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:1154360/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGGGCAGTATTGc                                                   <  1:2394371/99‑1 (MQ=255)

                                                 caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  1:135327/97‑1 (MQ=255)

                                                  aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                                               <  1:1180775/100‑1 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              >  1:85998/1‑100 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  1:2708537/1‑96 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:554749/100‑1 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  2:836393/100‑1 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           >  2:2224023/1‑100 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGAAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAAAGACTAGGCAGTTTGGTGCAGCATTGCTGTTGTTGa                                          >  1:2733308/1‑100 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGATTGGTGCAGTATTGCTGTTGTTGAt                                         >  1:669870/1‑100 (MQ=255)

                                                            tATGCTGCTGCCATTTTAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGATAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:2641095/99‑1 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    <  1:842783/100‑1 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    <  1:128791/100‑1 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CTCCTGTTCTGGATTATGACTAGCGTGTTTAGTGCAGTATTGCTGTTGTTGATTCCGa                                    <  1:872718/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  2:1211394/1‑100 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  <  1:432992/100‑1 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 <  2:206233/100‑1 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               <  2:1425933/100‑1 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  2:1407708/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGCAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:1699616/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             >  1:1090044/1‑100 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  2:3180/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:2691934/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGGTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:590267/100‑1 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            <  2:1126199/100‑1 (MQ=255)

                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           >  1:1297085/1‑100 (MQ=255)

                                                                       cATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                          <  1:2137428/100‑1 (MQ=255)

                                                                        nttGCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCCGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         >  1:1044052/2‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:2144698/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:2156395/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:2355324/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:1743942/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAAT‑CCCAC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:1304641/100‑1 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:1649899/100‑1 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    >  1:1885157/1‑100 (MQ=255)

                                                                               aGGAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  <  2:2007970/100‑1 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 <  1:900624/100‑1 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 <  1:1654967/100‑1 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  1:1663780/1‑100 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 <  1:919887/100‑1 (MQ=255)

                                                                                 gagaAGACTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                >  2:2786793/1‑100 (MQ=255)

                                                                                   gaagaCTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              >  2:1649824/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:749285/1‑100 (MQ=255)

                                                                                      gaCTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           <  2:1758662/100‑1 (MQ=255)

                                                                                      gaCTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  2:2564220/1‑100 (MQ=255)

                                                                                       aCTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  2:54290/100‑1 (MQ=255)

                                                                                        cTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  1:1688240/100‑1 (MQ=255)

                                                                                        cTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  2:176887/1‑100 (MQ=255)

                                                                                        cTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  2:1622140/1‑100 (MQ=255)

                                                                                        cTATGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGAGATTGTc          >  2:1624568/1‑99 (MQ=255)

                                                                                          aTGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       >  1:1050818/1‑100 (MQ=255)

                                                                                          aTGTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       <  1:2708171/100‑1 (MQ=255)

                                                                                            gTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:2126810/1‑100 (MQ=255)

                                                                                            gTTGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGATTGGTGCAGTATTGCTGTTGTAGATTCCGATAATGCACATTCAGCAGGTGATTGTCGaa       >  1:2455655/1‑98 (MQ=255)

                                                                                             ttGAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:419943/1‑100 (MQ=255)

                                                                                               gAAATCCCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  2:937048/100‑1 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑C‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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