BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,068,621	T→A	T407S (ACG→TCG)	rutG ←	pyrimidine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,068,621	0	T	A	93.6%	447.4 / 8.2	141	T407S (ACG→TCG)	rutG	pyrimidine permease
Reads supporting (aligned to +/- strand): ref base T (3/6); new base A (70/62); total (73/68)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

TGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTT  >  NC_000913/1068524‑1068718

 gaACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGaaaa                                                                                                >  1:3473723/1‑100 (MQ=255)

 gaACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGaaaa                                                                                                <  2:231796/100‑1 (MQ=255)

   acaacTTCTGGCGGTGGAGCGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAAc                                                                                              <  1:3429202/100‑1 (MQ=255)

   acaacTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACAAATCCCTCCCAACGAAAAAc                                                                                              >  1:2064587/1‑100 (MQ=255)

    caacTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAAcc                                                                                             <  1:1311479/100‑1 (MQ=255)

      acTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGc                                                                                           <  1:2537673/100‑1 (MQ=255)

      acTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGc                                                                                           >  1:3104153/1‑100 (MQ=255)

      acTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGc                                                                                           >  1:57397/1‑100 (MQ=255)

        ttCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGccc                                                                                         >  2:3378273/1‑100 (MQ=255)

         tCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCa                                                                                        >  1:569990/1‑100 (MQ=255)

          cTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAg                                                                                       >  1:3163512/1‑100 (MQ=255)

           tGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGc                                                                                      >  2:1274753/1‑100 (MQ=255)

            ggCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGc                                                                                      >  1:2049230/1‑99 (MQ=255)

             gCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGt                                                                                    <  2:1533782/100‑1 (MQ=255)

             gCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGt                                                                                    <  1:3414413/100‑1 (MQ=255)

              cGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTg                                                                                   <  1:1433353/100‑1 (MQ=255)

              cGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTg                                                                                   <  1:237058/100‑1 (MQ=255)

              cGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAAGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTg                                                                                   >  1:580346/1‑100 (MQ=255)

               ggTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTACAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTga                                                                                  <  1:1252309/100‑1 (MQ=255)

               ggTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTga                                                                                  >  2:3094622/1‑100 (MQ=255)

               ggTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTg                                                                                   >  1:2987251/1‑99 (MQ=255)

                ttGGAACTTCAACCAACTTGCTGCTTAACAACGCATCGTGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTga                                                                                  <  2:2025718/98‑1 (MQ=255)

                gTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTgag                                                                                 <  1:3092481/100‑1 (MQ=255)

                  ggAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGagag                                                                               <  2:1733614/100‑1 (MQ=255)

                  ggAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGagag                                                                               <  2:1733901/100‑1 (MQ=255)

                  ggAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCACAACGAAAAACCGCCCAGCGTGagag                                                                               <  1:1715890/100‑1 (MQ=255)

                  ggAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATAAAAAAAAAAAAAAAACCGCCCAGCGTGagag                                                                               <  1:994772/100‑1 (MQ=255)

                      cGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGAGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCaaa                                                                           >  2:2027980/1‑100 (MQ=255)

                       gTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCaaaa                                                                          >  2:2774468/1‑100 (MQ=255)

                         caaccaacTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATc                                                                        >  2:2608578/1‑100 (MQ=255)

                           accaacTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCAc                                                                      <  1:1452350/100‑1 (MQ=255)

                            ccaacTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCAcc                                                                     <  1:299483/100‑1 (MQ=255)

                                    cTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCa                                                             <  2:419224/100‑1 (MQ=255)

                                    cTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCa                                                             >  2:1877676/1‑100 (MQ=255)

                                     tGCGTAACAACGCATTGAGTAAAATCGCTCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCa                                                             <  1:879240/99‑1 (MQ=255)

                                      gCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGc                                                           <  1:684774/100‑1 (MQ=255)

                                      gCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGc                                                           >  1:424364/1‑100 (MQ=255)

                                       cTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGc                                                           <  1:1941262/99‑1 (MQ=255)

                                          aacaacGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCTCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTa                                                       <  1:1850656/100‑1 (MQ=255)

                                           acaacGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTaa                                                      >  1:184351/1‑100 (MQ=255)

                                            caacGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAg                                                     <  2:1436366/100‑1 (MQ=255)

                                            caacGCATTGAGTAAAATCGCGCCAAAGCTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAg                                                     <  2:2625848/100‑1 (MQ=255)

                                             accGCATTGACTAAAATCGCGCCAAAGGTTGCTGTACCAAAAAAAAAAAAAAAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAgg                                                    <  1:832880/98‑1 (MQ=255)

                                             aacGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAgg                                                    >  1:1518992/1‑100 (MQ=255)

                                             aacGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAgg                                                    >  1:394018/1‑100 (MQ=255)

                                             aacGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAgg                                                    <  1:2736926/100‑1 (MQ=255)

                                             aacGCATCGAGTAAAATCGCGCCAAAGGTTGCTGTACCAAAACCACACAAAGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAgg                                                    <  1:375398/100‑1 (MQ=255)

                                               cGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGtt                                                  <  2:476884/100‑1 (MQ=255)

                                               cGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGtt                                                  <  2:685217/100‑1 (MQ=255)

                                                gCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTa                                                 >  2:3320747/1‑100 (MQ=255)

                                                 cATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTAc                                                >  2:105530/1‑100 (MQ=255)

                                                  aTTGAGTAACATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACt                                               <  1:1613047/100‑1 (MQ=255)

                                                  aTTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACt                                               >  2:498514/1‑100 (MQ=255)

                                                   ttGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTg                                              >  1:1978078/1‑100 (MQ=255)

                                                       ggaaaaTCGCGCCAAAGGTTGCTGTACCAATCCATCACAAGGAAAAACCGCCCAGCGTGAACACAAAATCACCCACGACCAGCACTAAGGGTACTg                                              >  1:3184527/3‑96 (MQ=255)

                                                       gTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGAc                                          <  2:316813/100‑1 (MQ=255)

                                                       gTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGAc                                          <  2:964077/100‑1 (MQ=255)

                                                        tAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGAcc                                         <  2:3104203/100‑1 (MQ=255)

                                                         aaaaTCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCa                                        <  2:2487498/100‑1 (MQ=255)

                                                          aaaTCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCAt                                       >  1:2129859/1‑100 (MQ=255)

                                                           aaTCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATa                                      >  2:2639065/1‑100 (MQ=255)

                                                              cgcgCCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATc                                   <  1:1111148/100‑1 (MQ=255)

                                                               gcgcCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCa                                  <  2:363870/100‑1 (MQ=255)

                                                               gcgcCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCa                                  <  1:1047747/100‑1 (MQ=255)

                                                                cgcCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCaa                                 >  2:2640074/1‑100 (MQ=255)

                                                                 gcCAAAGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCaaa                                >  1:1053257/1‑100 (MQ=255)

                                                                   cAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCa                                  >  1:1524206/1‑96 (MQ=255)

                                                                     aaGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTAc                            <  1:771570/100‑1 (MQ=255)

                                                                     aaGGTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTAc                            >  1:1648886/1‑100 (MQ=255)

                                                                       ggTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTAcc                           <  1:1623828/99‑1 (MQ=255)

                                                                       ggTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCg                          >  1:721829/1‑100 (MQ=255)

                                                                       ggTTGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCg                          >  1:481701/1‑100 (MQ=255)

                                                                          tGCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTc                       >  2:294376/1‑100 (MQ=255)

                                                                           gCTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCt                      >  1:2883516/1‑100 (MQ=255)

                                                                            cTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTg                     >  2:2317679/1‑100 (MQ=255)

                                                                            cTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTg                     <  2:1338330/100‑1 (MQ=255)

                                                                            cTGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTg                     <  1:498509/100‑1 (MQ=255)

                                                                             tGTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTgg                    >  2:2432019/1‑100 (MQ=255)

                                                                              gTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGc                   >  2:386001/1‑100 (MQ=255)

                                                                              gTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGc                   >  2:102793/1‑100 (MQ=255)

                                                                              gTACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGc                   >  2:1972388/1‑100 (MQ=255)

                                                                               tACCAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAGGCTTACTGCGACCATAATCAAATTACCGTTCTGGCt                  >  2:2443196/1‑100 (MQ=255)

                                                                                 ccAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTa                <  2:721838/100‑1 (MQ=255)

                                                                                 ccAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTa                >  2:1936162/1‑100 (MQ=255)

                                                                                  cAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTaa               >  1:690060/1‑100 (MQ=255)

                                                                                  cAATCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTaa               >  2:1463463/1‑100 (MQ=255)

                                                                                   aaTCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTaaa              <  1:526694/100‑1 (MQ=255)

                                                                                   aaTCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTaaa              <  2:1680076/100‑1 (MQ=255)

                                                                                    aTCCCTCCCAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAAt             >  2:1101434/1‑100 (MQ=255)

                                                                                     tccctcccAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATc            <  2:1363733/100‑1 (MQ=255)

                                                                                     tccctcccAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATc            <  2:215007/100‑1 (MQ=255)

                                                                                      ccctcccAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCg           >  2:802130/1‑100 (MQ=255)

                                                                                      ccctcccAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCg           <  2:1713257/100‑1 (MQ=255)

                                                                                       cctcccAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGa          <  1:971515/100‑1 (MQ=255)

                                                                                        ctcccAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGac         >  2:1239742/1‑100 (MQ=255)

                                                                                        ctcccAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGac         >  1:2408760/1‑100 (MQ=255)

                                                                                        ctcccAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAACCAAATTACCGTTCTGGCTTAAATCGac         >  2:2040789/1‑100 (MQ=255)

                                                                                           ccAACGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACg      >  2:285237/1‑100 (MQ=255)

                                                                                             aaCGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACgg     >  1:2293559/1‑99 (MQ=255)

                                                                                             aaCGAAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGt    <  2:287831/100‑1 (MQ=255)

TGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTT  >  NC_000913/1068524‑1068718

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