BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,781,027	C→G	intergenic (‑151/+58)	ypjA ← / ← b2648	putative ATP‑binding component of a transport system/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,781,027	0	C	G	100.0%	670.0 / NA	218	intergenic (‑151/+58)	ypjA/b2648	putative ATP‑binding component of a transport system/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (95/123); total (95/123)

Only 100 of 219 total aligned reads displayed.

ACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGA‑G‑CCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTT  >  NC_000913/2780939‑2781116

aCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCaa                                                                                        <  2:160803/93‑1 (MQ=255)

aCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCaa                                                                                        >  1:2310004/1‑93 (MQ=255)

aCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAAGCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCaa                                                                                        <  1:260378/93‑1 (MQ=25)

  gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCaaaa                                                                                      <  1:4326093/93‑1 (MQ=255)

   cAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCaaaat                                                                                     >  2:2205892/1‑93 (MQ=255)

    aGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCaaaata                                                                                    <  2:807415/93‑1 (MQ=25)

       gcgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCaaaataaaa                                                                                 <  2:1111354/93‑1 (MQ=25)

       gcgcTGAAGTGTGGTTGGATAAAAAAGGCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCaaaataaaa                                                                                 <  1:1185183/93‑1 (MQ=25)

        cgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCaaaataaaat                                                                                <  2:1385406/93‑1 (MQ=25)

        cgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCaaaataaaat                                                                                <  1:1790589/93‑1 (MQ=25)

           tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGAt                                                                             >  2:2631845/1‑93 (MQ=255)

           tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGAt                                                                             >  1:1372146/1‑93 (MQ=255)

           tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGAt                                                                             <  2:4428750/93‑1 (MQ=25)

              aGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTa                                                                          >  1:3472739/1‑93 (MQ=255)

              aGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTa                                                                          <  1:2778744/93‑1 (MQ=37)

               gtgtGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTAc                                                                         <  1:3804592/93‑1 (MQ=37)

                 gtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACAt                                                                       <  1:1597768/93‑1 (MQ=37)

                  tgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATa                                                                      <  2:3450081/93‑1 (MQ=37)

                   ggTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATaa                                                                     <  1:93939/93‑1 (MQ=37)

                    gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAAt                                                                    >  2:1248130/1‑93 (MQ=255)

                       ggATAAAAAGGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGt                                                                 >  2:3143336/1‑93 (MQ=37)

                       ggATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGt                                                                 <  1:1588184/93‑1 (MQ=37)

                        gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                                <  1:3452068/93‑1 (MQ=37)

                        gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                                >  2:3233325/1‑93 (MQ=255)

                          tAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTc                                                              <  2:1860317/93‑1 (MQ=255)

                          tAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTc                                                              >  1:2868409/1‑93 (MQ=255)

                           aaaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCt                                                             <  1:3106092/93‑1 (MQ=38)

                            aaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtg                                                            <  2:1843220/93‑1 (MQ=255)

                             aaaaGACAATCCAATCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGA‑GTCCGCTGGCAAAATAAAATGATTTACATAATCGTTTCtga                                                           >  2:2717336/1‑93 (MQ=25)

                              aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                          <  2:4340763/93‑1 (MQ=37)

                              aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                          <  2:4480603/93‑1 (MQ=255)

                              aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                          <  1:570188/93‑1 (MQ=37)

                              aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                          <  1:798166/93‑1 (MQ=255)

                                 gTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGaa                                                       >  2:1148275/1‑93 (MQ=255)

                                  tCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAat                                                      <  2:3610890/93‑1 (MQ=25)

                                   cAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAata                                                     <  2:2003596/93‑1 (MQ=255)

                                    aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCATAATAAAATGATTTACATAATTGTTGCTGATGAatat                                                    >  2:4224642/1‑93 (MQ=21)

                                    aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                                    <  1:2013191/93‑1 (MQ=25)

                                     aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                                   >  2:3959300/1‑93 (MQ=255)

                                     aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                                   >  2:2582257/1‑93 (MQ=255)

                                       ccATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtctt                                                 <  2:1652639/93‑1 (MQ=255)

                                        cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                <  1:614793/93‑1 (MQ=255)

                                        cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                >  1:1221947/1‑93 (MQ=255)

                                        cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                <  1:3932286/93‑1 (MQ=25)

                                          ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                              >  1:222485/1‑93 (MQ=255)

                                          ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                              <  1:3989688/93‑1 (MQ=255)

                                          ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                              >  2:1541946/1‑93 (MQ=255)

                                            cAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCt                                            >  2:2321651/1‑93 (MQ=255)

                                               gAAATACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTcac                                         >  2:3281744/1‑93 (MQ=255)

                                                aaaTACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTcaca                                        <  1:2261603/93‑1 (MQ=21)

                                                 aaTACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACAt                                       >  2:4654077/1‑93 (MQ=255)

                                                  aTACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                                      >  1:2694153/1‑93 (MQ=18)

                                                  aTACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                                      >  1:2343350/1‑93 (MQ=18)

                                                  aTACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                                      >  2:1063023/1‑93 (MQ=18)

                                                  aTACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                                      <  1:4484353/93‑1 (MQ=255)

                                                   tACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaa                                     >  2:4526807/1‑93 (MQ=18)

                                                   tACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaa                                     >  1:843905/1‑93 (MQ=18)

                                                   tACGGGCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaa                                     >  1:3877031/1‑93 (MQ=18)

                                                      gggCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                                  <  2:2145000/93‑1 (MQ=18)

                                                      gggCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCAGATGAATATCTTCTGCTCACATaaaaa                                  <  2:205233/93‑1 (MQ=17)

                                                       ggCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                                 <  2:856964/93‑1 (MQ=255)

                                                       ggCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                                 <  1:4618913/93‑1 (MQ=255)

                                                       ggCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                                 <  1:3331867/93‑1 (MQ=255)

                                                       ggCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                                 <  1:1238786/93‑1 (MQ=18)

                                                        gCGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATc                                <  1:1912103/93‑1 (MQ=255)

                                                         cGTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATca                               >  2:3003012/1‑93 (MQ=255)

                                                          gTATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcac                              >  1:2496608/1‑93 (MQ=17)

                                                           tATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                             >  1:1162807/1‑93 (MQ=18)

                                                           tATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                             <  2:3861980/93‑1 (MQ=255)

                                                           tATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                             <  2:3164217/93‑1 (MQ=17)

                                                           tATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                             >  1:2591196/1‑93 (MQ=18)

                                                           tATTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                             >  2:598679/1‑93 (MQ=255)

                                                            aTTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                            <  1:3319376/93‑1 (MQ=255)

                                                            aTTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                            <  2:2841564/93‑1 (MQ=255)

                                                            aTTCTTTTCTTTCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                            <  1:2169056/93‑1 (MQ=255)

                                                              tcttttctttCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACaa                          >  2:4088562/1‑93 (MQ=255)

                                                               cttttctttCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACaaa                         >  2:4220081/1‑92 (MQ=14)

                                                                ttttctttCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa                        >  1:3067313/1‑93 (MQ=17)

                                                                ttttctttCGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa                        <  2:3699169/93‑1 (MQ=255)

                                                                 tttctttcGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATaa                       <  1:1233722/93‑1 (MQ=255)

                                                                 tttctttcGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATaa                       <  2:1936396/93‑1 (MQ=255)

                                                                  ttctttcGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAAc                      <  2:2241352/93‑1 (MQ=255)

                                                                   tctttcGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt                     <  1:4724182/93‑1 (MQ=255)

                                                                   tctttcGACAATGAGG‑CCGGTGGCAAAATAAAATGATTTACATAATCGTTTGTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt                     >  2:787237/1‑93 (MQ=14)

                                                                    ctttcGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt                    <  1:1897116/93‑1 (MQ=255)

                                                                     tttcGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt                   >  1:1660320/1‑93 (MQ=17)

                                                                     tttcGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt                   >  1:4045636/1‑93 (MQ=17)

                                                                      ttcGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTGCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTg                  <  2:3217021/93‑1 (MQ=14)

                                                                       tcGACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga                 <  1:1516781/93‑1 (MQ=255)

                                                                         gACAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgaga               >  1:257878/1‑93 (MQ=18)

                                                                          aCAATGAGG‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt              >  2:3846325/1‑93 (MQ=17)

                                                                               gagg‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAg         <  1:50410/90‑1 (MQ=255)

                                                                               gagg‑CCGTTGGCAAAATAAAATGATTTACACAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAg         >  1:1459082/4‑93 (MQ=21)

                                                                                agg‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGa        >  1:2099323/3‑93 (MQ=255)

                                                                                agg‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGa        <  2:662939/91‑1 (MQ=255)

                                                                                agg‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGa        <  1:4648614/91‑1 (MQ=255)

                                                                                agg‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGa        >  1:2480456/3‑93 (MQ=255)

                                                                                 gg‑CCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGAt       <  2:3949324/92‑1 (MQ=255)

                                                                                      gTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGtt   <  1:3046543/93‑1 (MQ=255)

                                                                                       ttGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGttt  <  1:1007558/93‑1 (MQ=255)

ACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGA‑G‑CCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTT  >  NC_000913/2780939‑2781116

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