BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,306,541	G→A	5.2%	A19A (GCC→GCT)	yjcS ←	alkyl sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,306,541	0	G	A	5.2%	806.2 / 17.8	249	A19A (GCC→GCT)	yjcS	alkyl sulfatase
Reads supporting (aligned to +/- strand): ref base G (126/110); new base A (7/6); total (133/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 271 total aligned reads displayed.

GGGCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCGAC  >  NC_000913/4306393‑4306690

gggCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTgg                                                                                                                                                      >  1:489294/1‑150 (MQ=255)

 ggCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggc                                                                                                                                                     <  1:1092219/150‑1 (MQ=255)

  gCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcg                                                                                                                                                    >  1:1105182/1‑150 (MQ=255)

  gCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcg                                                                                                                                                    <  1:453021/150‑1 (MQ=255)

  gCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcg                                                                                                                                                    <  2:681648/150‑1 (MQ=255)

  gCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcg                                                                                                                                                    <  2:489294/150‑1 (MQ=255)

   cGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcgg                                                                                                                                                   <  2:1193208/150‑1 (MQ=255)

   cGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcgg                                                                                                                                                   >  2:513914/1‑150 (MQ=255)

    gATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcggc                                                                                                                                                  <  1:105367/150‑1 (MQ=255)

    gATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcggc                                                                                                                                                  <  1:1203046/150‑1 (MQ=255)

    gATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcggc                                                                                                                                                  <  1:935747/150‑1 (MQ=255)

      tAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAg                                                                                                                                                <  1:1833364/150‑1 (MQ=255)

      tAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAg                                                                                                                                                <  2:815965/150‑1 (MQ=255)

      tAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAg                                                                                                                                                >  2:1573053/1‑150 (MQ=255)

      tAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAg                                                                                                                                                >  2:1635933/1‑150 (MQ=255)

       aaaTCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTTATTGGCTGCGTTGTCGTCCTCTTTCTCGTTAGCGGTATTTACTATCATGCCGCTGGCGGCAGt                                                                                                                                               <  2:156758/150‑1 (MQ=255)

       aaaTCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGt                                                                                                                                               <  2:173616/150‑1 (MQ=255)

       aaaTCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGt                                                                                                                                               >  2:894745/1‑150 (MQ=255)

       aaaTCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTAGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGt                                                                                                                                               >  2:1841626/1‑150 (MQ=255)

         aTCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatgat                                                                                                                                             <  1:1408193/150‑5 (MQ=255)

          tCCAGGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGTTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTaa                                                                                                                                            <  1:1934204/150‑1 (MQ=255)

           ccACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAAc                                                                                                                                           >  2:240859/1‑150 (MQ=255)

              cGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCa                                                                                                                                        <  1:240859/150‑1 (MQ=255)

                   gcgaTCGTCAGGATTTTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAAc                                                                                                                                   >  1:296649/4‑150 (MQ=255)

                    gCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACa                                                                                                                                  <  2:473881/150‑1 (MQ=255)

                        cGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATcc                                                                                                                              >  1:1594781/1‑150 (MQ=255)

                        cGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATcc                                                                                                                              >  1:760867/1‑150 (MQ=255)

                         gTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCt                                                                                                                             >  1:1732127/1‑150 (MQ=255)

                         gTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCt                                                                                                                             >  1:1674745/1‑150 (MQ=255)

                          tCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTg                                                                                                                            <  1:1703710/150‑1 (MQ=255)

                           cAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGc                                                                                                                           >  1:1325654/1‑150 (MQ=255)

                           cAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGc                                                                                                                           >  1:1470713/1‑150 (MQ=255)

                           cAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGc                                                                                                                           <  1:1279452/150‑1 (MQ=255)

                             aaaTCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTc                                                                                                                         <  2:76294/150‑1 (MQ=255)

                             aaaTCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTc                                                                                                                         >  1:1885024/1‑150 (MQ=255)

                                tCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCaaa                                                                                                                      <  1:1898901/150‑1 (MQ=255)

                                tCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGAGGCAGTTAACGCAATAACAATCCTGCTCaaa                                                                                                                      >  1:688011/1‑150 (MQ=255)

                                 ccTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAAc                                                                                                                     >  1:1571630/1‑150 (MQ=255)

                                  cTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatgatg                                                                                                                                            >  1:1259863/1‑121 (MQ=255)

                                   tGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACgg                                                                                                                   <  1:1379300/150‑1 (MQ=255)

                                    gACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCATACGGa                                                                                                                  <  2:309490/150‑1 (MQ=255)

                                    gACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGa                                                                                                                  >  1:434944/1‑150 (MQ=255)

                                      cGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAAt                                                                                                                >  2:124769/1‑150 (MQ=255)

                                       gATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGc                                                                                                                           <  1:578356/138‑1 (MQ=255)

                                          cGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAAcc                                                                                                            <  2:1344983/150‑1 (MQ=255)

                                           gCTAAACGGTAATGATTTGGTGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGACGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCg                                                                                                           >  2:1382255/1‑150 (MQ=255)

                                           gCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCg                                                                                                           <  2:1407157/150‑1 (MQ=255)

                                             tAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgag                                                                                                         <  2:526841/150‑1 (MQ=255)

                                              aaaCGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgaga                                                                                                        >  1:1493994/1‑150 (MQ=255)

                                               aaCGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGagag                                                                                                       >  1:1534225/1‑150 (MQ=255)

                                                aCGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGt                                                                                                      <  1:1860376/150‑1 (MQ=255)

                                                       tGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCat                                                                                               >  1:1590700/1‑150 (MQ=255)

                                                        gATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCata                                                                                              >  1:1016034/1‑150 (MQ=255)

                                                        gATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCata                                                                                              <  2:1083077/150‑1 (MQ=255)

                                                          tttGGCGTAATTCTGATTAAGCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTAGCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCAtata                                                                                            <  1:117558/150‑1 (MQ=255)

                                                          tttGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCAtata                                                                                            >  1:1012806/1‑150 (MQ=255)

                                                             ggCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAt                                                                                         >  1:1747245/1‑150 (MQ=255)

                                                              gCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAtc                                                                                        >  1:379863/1‑150 (MQ=255)

                                                              gCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAtc                                                                                        >  2:1435211/1‑150 (MQ=255)

                                                               cGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAtct                                                                                       >  2:1110165/1‑150 (MQ=255)

                                                               cGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAtct                                                                                       <  1:781725/150‑1 (MQ=255)

                                                                        gATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatc                                                                              >  1:345552/1‑150 (MQ=255)

                                                                           tAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcatt                                                                           <  2:1549461/150‑1 (MQ=255)

                                                                            aaCCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcattt                                                                          >  1:508177/1‑150 (MQ=255)

                                                                               cTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATg                                                                       >  2:263262/1‑150 (MQ=255)

                                                                                  ttGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCg                                                                    <  1:1196960/150‑1 (MQ=255)

                                                                                    gggTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACCATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTa                                                                  <  2:1289109/150‑1 (MQ=255)

                                                                                    gggTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAACGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTa                                                                  >  1:98761/1‑150 (MQ=255)

                                                                                     ggTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTaa                                                                 >  2:646743/1‑150 (MQ=255)

                                                                                      gTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTaaa                                                                <  2:602579/150‑1 (MQ=255)

                                                                                          atTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATACCCGAGCGTTATTCATATAAATCTCCATTTATCATTTATGTAGTAAATATg                                                            >  1:530814/1‑150 (MQ=255)

                                                                                          atTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATg                                                            <  2:497539/150‑1 (MQ=255)

                                                                                              aGTGGCGGCTTTCGCTTCCTCTTTCGCGTTATCGGTATTTACCATCATGCCGCTGGCGGCAGTTAAAGCAATAAAAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAAt                                                        <  1:1925094/150‑1 (MQ=255)

                                                                                              aGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAAt                                                        <  2:1131787/150‑1 (MQ=255)

                                                                                               gTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATc                                                       >  2:360925/1‑150 (MQ=255)

                                                                                                 ggcggcTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCa                                                     >  2:1912666/1‑150 (MQ=255)

                                                                                                    ggcTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGc                                                  <  1:764503/150‑1 (MQ=255)

                                                                                                         tCGCTTCCTCTTTCGCGTTAGCTGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAAt                                             <  2:924845/150‑1 (MQ=255)

                                                                                                               tctctTTCGCGGTATCGGTATTTTCCATCATTCCGGTGGGGGGATGTTACGGAATAACAATCCTGCTCAAACGGAATAACCGAGATTTATTTATATAAATCTCCATTTATGATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCa                                       <  2:1177568/149‑1 (MQ=255)

                                                                                                                 tctTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCAtt                                     <  2:177498/150‑1 (MQ=255)

                                                                                                                  ctTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAAACCTGCTCAAACGGAATAACCGAGAGTTATTCAAATAAATGTCCATATATCAGTTATGTCGTAAATATGTAATCCATGCAAAATAAGTCAttt                                    >  2:853613/1‑150 (MQ=255)

                                                                                                                   tttCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTa                                   <  1:872002/150‑1 (MQ=255)

                                                                                                                     tCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTaaa                                 <  2:1590700/150‑1 (MQ=255)

                                                                                                                      cgcgTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCACACGGAATAACCGAGAGTTATGCATATAAATCTCCATTTATCATTTATGTCGGAAATATGTAATCCATGCAAAATACGTCATTTAAAg                                >  1:886484/1‑150 (MQ=255)

                                                                                                                      cgcgTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAg                                <  2:225445/150‑1 (MQ=255)

                                                                                                                            aGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCg                          >  1:1799772/1‑150 (MQ=255)

                                                                                                                            aGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCg                          >  1:862358/1‑150 (MQ=255)

                                                                                                                                gTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTAt                      >  1:207835/1‑150 (MQ=255)

                                                                                                                                 tATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATc                     <  1:1841136/150‑1 (MQ=255)

                                                                                                                                    ttACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATa                  <  2:1680797/150‑1 (MQ=255)

                                                                                                                                     tACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCa                                       >  2:1840429/1‑128 (MQ=255)

                                                                                                                                     tACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATaa                 >  2:1302154/1‑150 (MQ=255)

                                                                                                                                     tACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATaa                 >  2:579836/1‑150 (MQ=255)

                                                                                                                                       gttaactGCCGCCAGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAAtt               >  2:1260985/7‑150 (MQ=255)

                                                                                                                                         atcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTgg             <  2:1747245/150‑1 (MQ=255)

                                                                                                                                           catGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGa           <  2:1493994/150‑1 (MQ=255)

                                                                                                                                            atGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGaa          >  1:918370/1‑150 (MQ=255)

                                                                                                                                             tGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGaaa         >  2:1951550/1‑150 (MQ=255)

                                                                                                                                              gCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAAt        >  1:812541/1‑150 (MQ=255)

                                                                                                                                                cGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTg      <  1:1979335/150‑1 (MQ=255)

GGGCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCGAC  >  NC_000913/4306393‑4306690

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