BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,958,258	1	.	C	70.7%	347.2 / 163.2	157	intergenic (‑827/‑34)	ppiC/rep	peptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand): ref base . (27/19); new base C (60/51); total (87/70)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 201 total aligned reads displayed.

CCGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CT‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTACCGGCCCCTGCC  >  NC_000913/3958166‑3958346

ccGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCCCATTCAAATGCGATTCTGCTAC‑AA‑t                                                                                                 >  2:5470311/1‑93 (MQ=255)

ccGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑t                                                                                                 >  1:532459/1‑93 (MQ=255)

ccGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑t                                                                                                 >  1:2289082/1‑93 (MQ=255)

 cGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑t‑‑c                                                                                              >  2:912489/1‑93 (MQ=255)

 cGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑t‑‑c                                                                                              <  1:6206404/93‑1 (MQ=255)

 cGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑t‑‑c                                                                                              >  2:4526575/1‑93 (MQ=255)

  gATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑t‑‑c‑c                                                                                            <  1:3129436/93‑1 (MQ=255)

  gATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑t‑‑c‑c                                                                                            <  1:3419073/93‑1 (MQ=255)

   aTTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑t‑‑c‑cc                                                                                           >  1:1944727/1‑92 (MQ=255)

     ttCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑t‑‑c‑cc‑c‑c                                                                                       <  2:3222436/93‑1 (MQ=255)

      tCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑c‑cc‑c‑cc                                                                                      <  1:1517268/93‑1 (MQ=255)

      tCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑c‑cc‑c‑cc                                                                                      <  2:3397415/93‑1 (MQ=255)

      tCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑c‑cc‑c‑cc                                                                                      <  2:5152281/93‑1 (MQ=255)

       cGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑cc‑c‑ccc                                                                                     >  1:3908619/1‑93 (MQ=255)

         gTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑cc‑c‑ccccc                                                                                   >  2:4730720/1‑93 (MQ=255)

          tcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑cc‑c‑cccccc                                                                                  <  2:5059227/93‑2 (MQ=255)

           cttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑cc‑c‑ccccccg                                                                                 >  1:866272/1‑91 (MQ=255)

            ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑cc‑c‑ccccccgt                                                                                >  2:5494984/1‑90 (MQ=255)

               tctctctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑g                                                                             <  2:5461188/93‑1 (MQ=255)

               tctctctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑Ga                                                                            >  1:4568525/1‑93 (MQ=255)

                 tctctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑Gaa                                                                           <  1:2366751/93‑1 (MQ=255)

                 tctctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAg                                                                          >  2:6033048/1‑93 (MQ=255)

                  ctctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAg                                                                          >  1:780701/1‑93 (MQ=255)

                  ctctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAg                                                                          <  2:5180754/93‑1 (MQ=255)

                   tctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGa                                                                         <  2:2286477/93‑1 (MQ=255)

                   tctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGa                                                                         >  2:1054372/1‑93 (MQ=255)

                   tctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑CACC‑C‑CCCCCGGTC‑GAAGa                                                                         >  2:1028135/1‑93 (MQ=255)

                     tGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGAtt                                                                       >  1:569238/1‑93 (MQ=255)

                     tGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGAtt                                                                       >  1:508031/1‑93 (MQ=255)

                      gATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGa                                                                     <  1:5813985/93‑1 (MQ=255)

                       atatAAACATCCGCAGCCAACCGGTTAGCGACTTACACGCGGTCACATACAAATGTGATTTTACTACAAA‑C‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGa                                                                     <  2:2096987/90‑1 (MQ=255)

                         ttAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGc                                                                   <  1:4460096/93‑1 (MQ=255)

                         ttAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCCTTC‑GAAGATTGAga                                                                   >  2:1924796/1‑92 (MQ=255)

                          tAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAgga                                                                  >  1:4938165/1‑91 (MQ=255)

                            aaCATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAAt                                                                <  2:296505/93‑1 (MQ=255)

                            aaCATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAAt                                                                >  1:1884909/1‑93 (MQ=255)

                            aaCATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCCTTC‑GCAGATTGTGCaca                                                                >  1:2233134/1‑91 (MQ=255)

                             aCATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATa                                                               >  2:1975649/1‑93 (MQ=255)

                              cATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATac                                                              <  2:5343172/93‑1 (MQ=255)

                              cATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATaa                                                              >  2:2008983/1‑92 (MQ=255)

                               aTCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATaca                                                             >  1:1652218/1‑93 (MQ=255)

                               aTCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATacac                                                            >  2:4689968/1‑93 (MQ=255)

                                tCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGCTTGGGAAATaaac                                                            >  2:394526/1‑93 (MQ=255)

                                tCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATacac                                                            <  1:3734966/93‑1 (MQ=255)

                                 ccGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACAcc                                                           <  2:4651119/93‑1 (MQ=255)

                                 ccGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACAcc                                                           <  1:1426890/93‑1 (MQ=255)

                                 ccGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCt                                                          <  1:3480332/93‑1 (MQ=255)

                                  cGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTa                                                         <  2:865186/93‑1 (MQ=255)

                                   gCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCtt                                                         >  2:4942509/1‑92 (MQ=255)

                                   gCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTAt                                                        >  1:4219998/1‑93 (MQ=255)

                                   ccAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTa                                                         <  2:3698042/92‑1 (MQ=255)

                                    cAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTAt                                                        <  2:1093673/93‑1 (MQ=255)

                                    cAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCCTGCAGAAGATTGAGCAAACCACatat                                                        >  2:5730856/1‑93 (MQ=255)

                                     aGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATg                                                       >  1:1873665/1‑93 (MQ=255)

                                     aGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATg                                                       <  1:5967720/93‑1 (MQ=255)

                                      gCCAACCGTTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGGTC‑GAAGATTGCGGAACAAACCTAttgg                                                     >  2:6138860/1‑90 (MQ=255)

                                      gCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCg                                                     >  2:2726824/1‑93 (MQ=255)

                                         aaCCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTTAGCAATACACCCATGCggc                                                   >  1:1293117/1‑91 (MQ=255)

                                         aaCCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTc                                                   <  2:80487/93‑1 (MQ=255)

                                             ggTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑CTCCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTaaa                                               >  1:3523298/1‑93 (MQ=255)

                                              gTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAAc                                              <  2:1672674/93‑1 (MQ=255)

                                              gTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAAc                                              >  1:1000689/1‑93 (MQ=255)

                                               ttAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTGC‑GAATGTTTAGCACTACACCTCTGCGGCTAAAcc                                             >  2:2488272/1‑93 (MQ=255)

                                               ttAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAAccc                                            <  2:1620070/93‑1 (MQ=255)

                                               gtAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TCCC‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAAc                                              <  1:180257/92‑1 (MQ=255)

                                                      cTTACACGTGGTCACATTCTAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGACCAAAAAACCAATGCGTCTAAACCCCCGCca                                      >  2:6111400/1‑93 (MQ=255)

                                                       ttACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCcaac                                    >  1:3046043/1‑93 (MQ=255)

                                                       ttACACGCGCTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCcaac                                    >  2:4134962/1‑93 (MQ=255)

                                                        tACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGCGTGAGCAATACACCTATTCGTCTTAACCCCGGCcaac                                    >  1:5563755/1‑93 (MQ=255)

                                                            cgcgGTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCcaacaac‑a                               >  1:1767802/1‑93 (MQ=255)

                                                               ggTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGc                            >  1:4015284/1‑93 (MQ=255)

                                                               ggTCACATTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATAAAACTATTCGTCTAAAACCCAGCCCACAAA‑AAGc                            >  2:3634082/1‑93 (MQ=255)

                                                                gTCACATTCAAATGCGATTCTGCTAC‑AACC‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCt                           <  1:4439140/93‑1 (MQ=255)

                                                                gTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CCCC‑CCCCCGTTC‑GAAAATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAAAAAC‑AAGCt                           >  2:2678001/1‑93 (MQ=255)

                                                                gTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCccacaaa‑aaaaag                          >  1:4651800/1‑89 (MQ=255)

                                                                 tCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCCGCCAACAAC‑AAACTGt                         >  1:1206886/1‑93 (MQ=255)

                                                                  cacaTTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGt                         <  1:2736851/93‑1 (MQ=255)

                                                                    caTTCAAATGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATTCGTCTAAACCCCGGCCAACAAC‑AAGCTgccc                       >  1:4770735/1‑90 (MQ=255)

                                                                    caTTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGa                      >  1:5114309/1‑93 (MQ=255)

                                                                      ttaaaaTGCAATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGaa                     <  2:3969260/90‑1 (MQ=255)

                                                                          aaTGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCCTTC‑CAAGGCTTCGCAATACCCCTATACCGCTAAACCCCGGCCAACAAC‑AATCTGGCGAATTCg                 >  2:6005001/1‑93 (MQ=255)

                                                                           aTGCGATTGTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGt                <  1:5853371/93‑1 (MQ=255)

                                                                           aTGCGATTCTGCTAC‑AACC‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGt                <  2:4577520/93‑1 (MQ=255)

                                                                           aTGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTCAGAAAAACACCTTGGGGTTTAAACCCCGGCC‑ACAACAAAGCTGCCGAATTCGt                >  2:1857355/1‑93 (MQ=255)

                                                                           aTGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCATTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGt                >  1:2830557/1‑93 (MQ=255)

                                                                            tGCGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGtt               <  2:2859425/93‑1 (MQ=255)

                                                                             gCGATTCTGCTAC‑AACC‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTa              <  2:3531313/93‑1 (MQ=255)

                                                                             gCGATTCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTAc             <  1:6103262/93‑1 (MQ=255)

                                                                              cGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTAc             >  1:930036/1‑93 (MQ=255)

                                                                              cGATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTAc             <  1:3151731/93‑1 (MQ=255)

                                                                               gATTCTGCTAC‑AACC‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTAcc            <  2:1375154/93‑1 (MQ=255)

                                                                               gATTCTGCTAC‑AACC‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTAcc            <  1:5403388/93‑1 (MQ=255)

                                                                               gATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGTTTGAGCAAAACACCAATGCGTCTAAACCCCGGCCAACAAC‑AAGCTTTCGAATTGGTTaac            >  2:4876620/1‑91 (MQ=255)

                                                                               gATTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTAcc            <  1:1059755/93‑1 (MQ=255)

                                                                                aTTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTACCg           <  2:3181492/93‑1 (MQ=255)

                                                                                aTTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTACCg           <  2:1767840/93‑1 (MQ=255)

                                                                                aTTCTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTACCg           >  1:2116777/1‑93 (MQ=255)

                                                                                 ttCTGCTAC‑AA‑T‑‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTACCGGc         >  1:3348598/1‑93 (MQ=255)

                                                                                   cTGCTAC‑AA‑TC‑C‑CC‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTACCGGcc        >  1:1919533/1‑93 (MQ=255)

                                                                                      cTAC‑AA‑T‑‑C‑CCCC‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAAAAAC‑CAGCTGTCGAATTTGTTACCCGccccc     >  1:2560929/1‑92 (MQ=255)

CCGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTAC‑AA‑T‑‑C‑CT‑C‑CCCCCGTTC‑GAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAAC‑AAGCTGTCGAATTCGTTACCGGCCCCTGCC  >  NC_000913/3958166‑3958346

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