BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,204,942	T→C	pseudogene (539/1134 nt)	aaaE →	pseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,204,942	0	T	C	100.0%	369.2 / NA	100	pseudogene (539/1134 nt)	aaaE	pseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (57/43); total (57/43)

Only 100 of 101 total aligned reads displayed.

GAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGAAAAATGA  >  NC_000913/1204845‑1205038

gAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACctt                                                                                                <  1:1359481/100‑1 (MQ=255)

  ggATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTtat                                                                                              <  2:762045/100‑1 (MQ=255)

  ggATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTtat                                                                                              >  2:1596676/1‑100 (MQ=255)

   gATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTtata                                                                                             <  2:2118263/100‑1 (MQ=255)

    aTAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTtatgc                                                                                            >  1:1711563/1‑98 (MQ=255)

    aTAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATAc                                                                                            <  2:955817/100‑1 (MQ=255)

     tAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATAcc                                                                                           >  2:419084/1‑100 (MQ=255)

      aaaTCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATAccc                                                                                          >  2:1453672/1‑100 (MQ=255)

      aaaTCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATAcc                                                                                           >  1:37223/1‑99 (MQ=255)

       aaTCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCcg                                                                                         >  1:1847941/1‑100 (MQ=255)

        aTCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCcgc                                                                                        <  1:2145386/100‑1 (MQ=255)

          cACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCga                                                                                      >  1:468226/1‑100 (MQ=255)

            cTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCgaga                                                                                    >  1:2372499/1‑100 (MQ=255)

             tGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGAt                                                                                   <  2:1215952/100‑1 (MQ=255)

               aCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTg                                                                                 >  2:962401/1‑100 (MQ=255)

               aCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTg                                                                                 >  1:951326/1‑100 (MQ=255)

               aCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTg                                                                                 <  2:2717004/100‑1 (MQ=255)

                cccTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGa                                                                                >  1:863126/1‑100 (MQ=255)

                  cTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACg                                                                              <  1:2728525/100‑1 (MQ=255)

                   ttGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACgg                                                                             <  2:1076557/100‑1 (MQ=255)

                    tGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGt                                                                            <  2:1331997/100‑1 (MQ=255)

                      aGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTaa                                                                          >  2:1268603/1‑100 (MQ=255)

                       gCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTaaa                                                                         >  1:1450017/1‑100 (MQ=255)

                       gCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTaaa                                                                         >  1:783568/1‑100 (MQ=255)

                        cAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTaaaa                                                                        <  1:1767564/100‑1 (MQ=255)

                         aGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAAc                                                                       >  1:1389773/1‑100 (MQ=255)

                           ttCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGc                                                                     <  1:1685481/100‑1 (MQ=255)

                             cGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCAt                                                                   >  2:2187446/1‑100 (MQ=255)

                                 tggCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATtact                                                               <  1:1589886/99‑1 (MQ=255)

                                  ggCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTAGACCCGCGAGATTGACGGTAAAACGCATtacta                                                              >  2:1589898/1‑100 (MQ=255)

                                    cAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACa                                                            <  1:1268595/100‑1 (MQ=255)

                                     aGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAg                                                           <  1:1105609/100‑1 (MQ=255)

                                      gCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAgt                                                          >  2:850903/1‑100 (MQ=255)

                                      gCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAgt                                                          >  1:1326655/1‑100 (MQ=255)

                                      gCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAgt                                                          >  2:148242/1‑100 (MQ=255)

                                         gTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGtgtg                                                       >  1:586952/1‑100 (MQ=255)

                                          tGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTgg                                                      >  1:2817764/1‑100 (MQ=255)

                                           gCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGc                                                     <  1:1453660/100‑1 (MQ=255)

                                            cATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGcc                                                    >  2:198048/1‑100 (MQ=255)

                                             aTTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGccc                                                   >  1:1246295/1‑100 (MQ=255)

                                             aTTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGc                                                     <  2:747432/98‑1 (MQ=255)

                                              ttCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGcccc                                                  <  2:1088104/100‑1 (MQ=255)

                                              ttCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGcccc                                                  >  1:179087/1‑100 (MQ=255)

                                              ttCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGcccc                                                  >  2:219295/1‑100 (MQ=255)

                                                cTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCcgc                                                >  1:1173602/1‑100 (MQ=255)

                                                  ggCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGt                                              >  2:958988/1‑100 (MQ=255)

                                                   gCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCGTGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGtt                                             >  2:54820/1‑100 (MQ=255)

                                                   gCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCAAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGtt                                             >  2:1252111/1‑100 (MQ=255)

                                                      tttGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCt                                          <  2:873476/100‑1 (MQ=255)

                                                       ttGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCt                                          >  1:1082061/1‑99 (MQ=255)

                                                        tGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTgg                                        <  1:1794597/100‑1 (MQ=255)

                                                         gACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTggg                                       >  2:1888563/1‑100 (MQ=255)

                                                         gACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTggg                                       <  2:1970074/100‑1 (MQ=255)

                                                         gACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTggg                                       <  2:2303564/100‑1 (MQ=255)

                                                         gACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTggg                                       <  2:1749477/100‑1 (MQ=255)

                                                         gACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTggg                                       >  1:814703/1‑100 (MQ=255)

                                                          aCCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGt                                      <  2:1326663/100‑1 (MQ=255)

                                                           ccTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTa                                     >  2:2802796/1‑100 (MQ=255)

                                                            cTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTAc                                    >  1:2195176/1‑100 (MQ=255)

                                                            cTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTAc                                    >  1:1641414/1‑100 (MQ=255)

                                                             tGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTAcc                                   <  1:2187430/100‑1 (MQ=255)

                                                               gcgcGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCgn                                 <  2:463887/100‑2 (MQ=255)

                                                                cgcgTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTa                                <  1:126622/100‑1 (MQ=255)

                                                                cgcgTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTa                                <  2:1711576/100‑1 (MQ=255)

                                                                 gcgTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTAt                               >  2:1166757/1‑100 (MQ=255)

                                                                 gcgTAAGCTGGAGATTAACAGCATGGCGCGACCTTATCCCCGCGAGATTGACTGTAAAATGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCg                                  <  1:564224/97‑1 (MQ=255)

                                                                  cgTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATg                              <  2:368928/100‑1 (MQ=255)

                                                                   gTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCTGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGa                             >  2:948039/1‑100 (MQ=255)

                                                                       gCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACg                         >  1:700799/1‑100 (MQ=255)

                                                                       gCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACg                         <  2:2683256/100‑1 (MQ=255)

                                                                       gCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACg                         >  2:2575767/1‑100 (MQ=255)

                                                                        cTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACgg                        >  2:1101629/1‑100 (MQ=255)

                                                                        cTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACgg                        >  2:1101401/1‑100 (MQ=255)

                                                                        cTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACgg                        <  1:1596664/100‑1 (MQ=255)

                                                                        cTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACgg                        <  1:845395/100‑1 (MQ=255)

                                                                         tGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGgt                       <  2:783574/100‑1 (MQ=255)

                                                                           gATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGgtgt                     <  2:264354/100‑1 (MQ=255)

                                                                           gATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGgtgt                     <  2:2723250/100‑1 (MQ=255)

                                                                           gATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGgtgt                     <  2:1138604/100‑1 (MQ=255)

                                                                            atatGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTa                    >  2:2652003/1‑100 (MQ=255)

                                                                             tatGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTAc                   <  1:198145/100‑1 (MQ=255)

                                                                               tGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGTATTACTACAGTGTGGCCCCGCGTTTCTAGGTACCGTATGACACGGTGTACAg                 >  1:2287008/1‑100 (MQ=255)

                                                                               tGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAg                 >  1:7620/1‑100 (MQ=255)

                                                                                gAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGc                >  1:1994580/1‑100 (MQ=255)

                                                                                gAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGc                <  2:1921470/100‑1 (MQ=255)

                                                                                 aaCAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCg               >  1:366016/1‑100 (MQ=255)

                                                                                  aCAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCg               >  2:412940/1‑99 (MQ=255)

                                                                                  aCAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGt              <  1:338955/100‑1 (MQ=255)

                                                                                    aGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCg            <  1:958981/100‑1 (MQ=255)

                                                                                     gCATGGCGCGACCTTATACCCGCGAGATTGATGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCga           >  1:688333/1‑100 (MQ=255)

                                                                                     gCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCga           <  2:1104473/100‑1 (MQ=255)

                                                                                     gCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCga           >  2:670398/1‑100 (MQ=255)

                                                                                     gCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCga           >  2:773296/1‑100 (MQ=255)

                                                                                       aTGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCgaga         >  2:2656199/1‑100 (MQ=255)

                                                                                       aTGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCgaga         <  2:1414550/100‑1 (MQ=255)

                                                                                          gtgcgACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGaaaa      >  2:1903813/3‑100 (MQ=255)

                                                                                          gcgcGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGaaaa      >  1:830981/1‑100 (MQ=255)

                                                                                           cgcgACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGaaaaa     <  1:592461/100‑1 (MQ=255)

                                                                                           cgcgACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGaaaaa     >  1:353898/1‑100 (MQ=255)

                                                                                              gACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGAAAAATGa  >  1:810670/1‑100 (MQ=255)

GAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGAAAAATGA  >  NC_000913/1204845‑1205038

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