BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,441,588	C→A	G197C (GGT→TGT)	msrA ←	methionine sulfoxide reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,441,588	0	C	A	100.0%	1451.8 / NA	335	G197C (GGT→TGT)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (160/175); total (160/175)

Only 100 of 335 total aligned reads displayed.

CATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCATCG  >  NC_000913/4441490‑4441685

cATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAAt                                                                                                 >  1:6930799/1‑101 (MQ=255)

 aTAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATa                                                                                                <  2:6930722/101‑1 (MQ=255)

 aTAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATa                                                                                                >  2:6930802/1‑101 (MQ=255)

  tAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATAc                                                                                               >  1:6930806/1‑101 (MQ=255)

  tAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATAc                                                                                               >  2:6930807/1‑101 (MQ=255)

    gACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACgg                                                                                             <  1:6930444/101‑1 (MQ=255)

    gACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACgg                                                                                             <  1:6930180/101‑1 (MQ=255)

       aaaGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGtt                                                                                          >  1:6930814/1‑101 (MQ=255)

       aaaGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGtt                                                                                          >  2:6930813/1‑101 (MQ=255)

         agagCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGtttt                                                                                        >  2:6930819/1‑101 (MQ=255)

         agagCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGtttt                                                                                        <  2:6930538/101‑1 (MQ=255)

         agagCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGtttt                                                                                        <  1:6930566/101‑1 (MQ=255)

          gagCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGttttt                                                                                       >  1:6930821/1‑101 (MQ=255)

             cACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATg                                                                                    >  2:6930827/1‑101 (MQ=255)

              aCTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGc                                                                                   >  2:6930832/1‑101 (MQ=255)

              aCTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGc                                                                                   <  1:6930219/101‑1 (MQ=255)

              aCTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGc                                                                                   >  1:6930828/1‑101 (MQ=255)

                 ttATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGa                                                                                <  1:6930054/101‑1 (MQ=255)

                  tATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGat                                                                               <  2:6930428/101‑1 (MQ=255)

                  tATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGat                                                                               >  2:6930836/1‑101 (MQ=255)

                   aTCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGata                                                                              <  2:6930335/101‑1 (MQ=255)

                    tCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGatat                                                                             <  1:6930306/101‑1 (MQ=255)

                    tCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGatat                                                                             >  1:6930838/1‑101 (MQ=255)

                    tCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGatat                                                                             <  1:6930398/101‑1 (MQ=255)

                     caacaaACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATAtt                                                                            <  1:6930359/101‑1 (MQ=255)

                     caacaaACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATAtt                                                                            >  2:6930844/1‑101 (MQ=255)

                     caacaaACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATAtt                                                                            <  1:6930407/101‑1 (MQ=255)

                      aacaaACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACggg                                                                                            <  2:6930491/84‑1 (MQ=255)

                      aacaaACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTg                                                                           <  1:6930341/101‑1 (MQ=255)

                      aacaaACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTg                                                                           >  2:6930848/1‑101 (MQ=255)

                      aacaaACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTg                                                                           <  1:6930439/101‑1 (MQ=255)

                       acaaACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCATACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGc                                                                          >  1:6930849/1‑101 (MQ=255)

                        caaACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCt                                                                         >  2:6930852/1‑101 (MQ=255)

                         aaaCTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCtg                                                                        <  2:6930190/101‑1 (MQ=255)

                         aaaCTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCtg                                                                        >  1:6930855/1‑101 (MQ=255)

                           aCTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCtggt                                                                      <  2:6930474/101‑1 (MQ=255)

                           aCTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCtggt                                                                      <  1:6930469/101‑1 (MQ=255)

                            cTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCtggtg                                                                     <  2:6930245/101‑1 (MQ=255)

                             tACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCtggtgg                                                                    <  1:6930249/101‑1 (MQ=255)

                              aCATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTggtggt                                                                   >  1:6930862/1‑101 (MQ=255)

                                aTTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCa                                                                 >  2:6930864/1‑101 (MQ=255)

                                aTTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCa                                                                 >  1:6930867/1‑101 (MQ=255)

                                 tttGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCAt                                                                <  2:6930453/101‑1 (MQ=255)

                                  ttGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATc                                                               <  1:6930293/101‑1 (MQ=255)

                                   tGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCt                                                              >  2:6930873/1‑101 (MQ=255)

                                   tGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCt                                                              <  2:6930412/101‑1 (MQ=255)

                                   tGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCt                                                              >  1:6930872/1‑101 (MQ=255)

                                       cccGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                                          >  2:6930879/1‑101 (MQ=255)

                                          gTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCat                                                       >  1:6930882/1‑101 (MQ=255)

                                           tAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCata                                                      >  2:6930885/1‑101 (MQ=255)

                                             aCGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCataat                                                    <  2:6930524/101‑1 (MQ=255)

                                               gCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCa                                                                                  >  1:6930891/1‑69 (MQ=255)

                                               gCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCAtaataa                                                  >  1:6930888/1‑101 (MQ=255)

                                                cTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaa                                                 >  2:6930893/1‑101 (MQ=255)

                                                cTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaa                                                 <  1:6930678/101‑1 (MQ=255)

                                                    gCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACgg                                             >  2:6930898/1‑101 (MQ=255)

                                                     cTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGgt                                            <  2:6930095/101‑1 (MQ=255)

                                                      ttCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGgtg                                           >  2:6930904/1‑101 (MQ=255)

                                                      ttCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGgtg                                           >  1:6930902/1‑101 (MQ=255)

                                                       tCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGgtgt                                          <  2:6930388/101‑1 (MQ=255)

                                                        ccGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGtgtg                                         <  1:6930537/101‑1 (MQ=255)

                                                         cggcggCAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTgg                                        >  1:6930908/1‑101 (MQ=255)

                                                          ggcggcAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGc                                       >  1:6930911/1‑101 (MQ=255)

                                                            cggcAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGt                                     <  1:6930133/101‑1 (MQ=255)

                                                             ggcAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGtt                                    <  1:6929907/101‑1 (MQ=255)

                                                             ggcAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGtt                                    <  2:6930442/101‑1 (MQ=255)

                                                              gcAGACAGACGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTa                                   >  2:6930915/1‑101 (MQ=255)

                                                                agacagacGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGc                                 >  1:6930920/1‑101 (MQ=255)

                                                                agacagacGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGc                                 <  2:6930355/101‑1 (MQ=255)

                                                                agacagacGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGc                                 >  2:6930919/1‑101 (MQ=255)

                                                                    agacGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGAtt                             >  1:6930923/1‑101 (MQ=255)

                                                                      acGCCAATTCCGCCAATTCCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTc                           <  1:6930499/101‑1 (MQ=255)

                                                                       cgccaattccgccaattcCACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTcc                          <  2:6930408/101‑1 (MQ=255)

                                                                         ccaattccgccaattccACAGTAACAATACGGGTTTTTATGCa                                                                                  <  1:6930516/43‑1 (MQ=255)

                                                                          caattccgccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGt                                     <  1:6930334/87‑1 (MQ=255)

                                                                          caattccgccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtg                       <  2:6930484/101‑1 (MQ=255)

                                                                          caattccgccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtg                       >  2:6930929/1‑101 (MQ=255)

                                                                          caattccgccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtg                       <  2:6930511/101‑1 (MQ=255)

                                                                          caattccgccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtg                       >  1:6930930/1‑101 (MQ=255)

                                                                            attccgccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtggt                     >  2:6930935/1‑101 (MQ=255)

                                                                              tccgccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGTTGa                   <  2:6930567/101‑1 (MQ=255)

                                                                              tccgccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGa                   >  1:6930938/1‑101 (MQ=255)

                                                                                cgccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGAtg                 <  2:6930371/101‑1 (MQ=255)

                                                                                 gccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGAtgt                <  1:6930423/101‑1 (MQ=255)

                                                                                 gccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGAtgt                <  1:6930257/101‑1 (MQ=255)

                                                                                  ccaattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGAtgtg               <  1:6930520/101‑1 (MQ=255)

                                                                                   caattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGa              >  2:6930941/1‑101 (MQ=255)

                                                                                   caattccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGa              >  1:6930943/1‑101 (MQ=255)

                                                                                     attccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACg            <  2:6930221/101‑1 (MQ=255)

                                                                                      ttccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGa           <  1:6930370/101‑1 (MQ=255)

                                                                                      ttccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGa           <  1:6930279/101‑1 (MQ=255)

                                                                                      ttccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGa           >  2:6930948/1‑101 (MQ=255)

                                                                                       tccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGAt          <  2:6930642/101‑1 (MQ=255)

                                                                                       tccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGAt          <  1:6930487/101‑1 (MQ=255)

                                                                                       tccACAGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGAt          <  1:6930381/101‑1 (MQ=255)

                                                                                         cacaGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCg        >  1:6930954/1‑101 (MQ=255)

                                                                                         cacaGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCg        >  2:6930952/1‑101 (MQ=255)

                                                                                            aGTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCa     <  1:6930633/101‑1 (MQ=255)

                                                                                             gTAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCAt    <  2:6930618/101‑1 (MQ=255)

                                                                                              tAACAATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCATc   <  2:6930297/101‑1 (MQ=255)

CATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCATCG  >  NC_000913/4441490‑4441685

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