BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,650,355	T→C	pseudogene (781/1158 nt)	intQ →	pseudogene, Qin prophage; predicted defective integrase;Phage or Prophage Related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,650,355	0	T	C	100.0%	1308.4 / NA	294	pseudogene (781/1158 nt)	intQ	pseudogene, Qin prophage; predicted defective integrase;Phage or Prophage Related
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (134/160); total (134/160)

Only 100 of 296 total aligned reads displayed.

CGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAT  >  NC_000913/1650262‑1650444

cGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACt                                                                                           >  1:1153047/1‑94 (MQ=255)

    tcattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCac                                                                                       >  1:6156185/1‑93 (MQ=255)

     aattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCaga                                                                                      <  1:4733941/93‑1 (MQ=255)

      attcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagat                                                                                     >  1:2257203/1‑93 (MQ=255)

       ttcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagagt                                                                                    >  1:3836745/1‑93 (MQ=255)

       ttcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagaga                                                                                    <  1:917479/94‑1 (MQ=255)

         catCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagagaga                                                                                  >  1:516931/1‑93 (MQ=255)

          ttCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTa                                                                                 <  1:841405/93‑1 (MQ=255)

            tttattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATg                                                                               <  1:5990408/93‑1 (MQ=255)

             ttattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATga                                                                              >  1:2321620/1‑93 (MQ=255)

               gttaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGca                                                                            <  1:6551459/93‑1 (MQ=255)

               attaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCa                                                                                        <  2:2865666/82‑1 (MQ=255)

                ttaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcag                                                                           >  1:188032/1‑94 (MQ=255)

                 taAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCaga                                                                                      >  2:7053255/1‑82 (MQ=255)

                  aaGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagag                                                                                     <  2:5663693/82‑1 (MQ=255)

                   aGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagaga                                                                                    <  2:4176218/82‑1 (MQ=255)

                   aGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaat                                                                        >  1:2138175/1‑93 (MQ=255)

                   aGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaaa                                                                        >  1:4589627/1‑93 (MQ=255)

                    gCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagagag                                                                                   <  2:4940419/82‑1 (MQ=255)

                      cAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACAGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacagg                                                                     >  1:3236082/1‑93 (MQ=255)

                       aGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTAt                                                                                <  2:1218474/82‑1 (MQ=255)

                        gCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATg                                                                               <  2:662229/82‑1 (MQ=255)

                        acAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaa                                                                   <  1:3477653/93‑1 (MQ=255)

                         cAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATgg                                                                              <  2:4452611/82‑1 (MQ=255)

                         cAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATgg                                                                              >  2:2214398/1‑82 (MQ=255)

                          aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGc                                                                             >  2:2840978/1‑82 (MQ=255)

                           aTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGca                                                                            >  2:6156984/1‑82 (MQ=255)

                              taCGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaa                                                             <  1:1083983/93‑1 (MQ=255)

                              gACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaat                                                             >  1:1433400/1‑93 (MQ=255)

                               aCGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaat                                                            >  1:3964184/1‑93 (MQ=255)

                               aCGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaac                                                            >  1:5380738/1‑93 (MQ=255)

                                cGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaaat                                                           >  1:4403398/1‑93 (MQ=255)

                                 tctctCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                                          <  1:135232/93‑1 (MQ=255)

                                   gctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGc                                                        <  1:5942203/93‑1 (MQ=255)

                                     tcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACaaaaaa                                                                  >  2:161583/1‑82 (MQ=255)

                                      ccGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaaa                                                                 <  2:1072300/82‑1 (MQ=255)

                                       gggAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCt                                                    <  1:122679/93‑1 (MQ=255)

                                       cGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAAc                                                                <  2:3344847/82‑1 (MQ=255)

                                        agAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCtt                                                   <  1:1916608/93‑1 (MQ=255)

                                         gAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCttt                                                  <  1:3645295/94‑1 (MQ=255)

                                          aGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaa                                                             >  2:3049667/1‑82 (MQ=255)

                                           gTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTgg                                                >  1:4197625/1‑93 (MQ=255)

                                           gTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGt                                                <  1:5705364/94‑1 (MQ=255)

                                            tCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaaaa                                                           <  2:5587497/82‑1 (MQ=255)

                                            tCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtt                                               >  1:48653/1‑94 (MQ=255)

                                             cAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttt                                              <  1:6076818/94‑1 (MQ=255)

                                              aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATg                                                         <  2:3165744/82‑1 (MQ=255)

                                              aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATg                                                         >  2:4589944/1‑82 (MQ=255)

                                              aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtttg                                             >  1:4055375/1‑93 (MQ=255)

                                                aTGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCa                                                       <  2:5229058/82‑1 (MQ=255)

                                                aTGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCa                                                       <  2:6464141/82‑1 (MQ=255)

                                                  taCATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCaa                                         <  1:2538174/93‑1 (MQ=255)

                                                     gttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                      <  1:4644137/93‑1 (MQ=255)

                                                     aTTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCttt                                                  >  2:6515168/1‑82 (MQ=255)

                                                      ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                                                 <  2:3241779/82‑1 (MQ=255)

                                                      ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                                                 >  2:3421034/1‑82 (MQ=255)

                                                      ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                                                 <  2:148119/82‑1 (MQ=255)

                                                       tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAGAACAAAAATGCACCTTTGt                                                >  2:1010817/1‑82 (MQ=255)

                                                       tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGt                                                >  2:70164/1‑82 (MQ=255)

                                                       tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGa                                    >  1:5308526/1‑94 (MQ=255)

                                                         ccGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAg                                  <  1:3350002/93‑1 (MQ=255)

                                                         aCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttt                                              >  2:3001548/1‑82 (MQ=255)

                                                           gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttttt                                            >  2:5538957/1‑82 (MQ=255)

                                                           gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttttt                                            <  2:3540802/82‑1 (MQ=255)

                                                           gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgta                                >  1:4498473/1‑93 (MQ=255)

                                                            ggaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgtgt                               <  1:2868477/93‑1 (MQ=255)

                                                             gaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgtgta                              >  1:7299911/1‑93 (MQ=255)

                                                              aCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCaa                                         >  2:6749300/1‑82 (MQ=255)

                                                                tGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAcc                                       >  2:505536/1‑82 (MQ=255)

                                                                tGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAcc                                       <  2:225552/82‑1 (MQ=255)

                                                                tGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAcc                                       <  2:5971771/82‑1 (MQ=255)

                                                                  aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTg                                     >  2:6547392/1‑82 (MQ=255)

                                                                  aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTg                                     <  2:4365596/82‑1 (MQ=255)

                                                                  aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTg                                     >  2:2614683/1‑82 (MQ=255)

                                                                  aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTg                                     >  2:2031256/1‑82 (MQ=255)

                                                                     caaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGagag                      <  1:3085426/93‑1 (MQ=255)

                                                                       ggagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTa                    <  1:2994079/93‑1 (MQ=255)

                                                                       agagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgtg                                >  2:4026326/1‑82 (MQ=255)

                                                                       agagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGtc                    >  1:1501677/1‑93 (MQ=255)

                                                                       agagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTa                    >  1:6501333/1‑94 (MQ=255)

                                                                         agCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTc                              >  2:2512194/1‑82 (MQ=255)

                                                                           cATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAg                            >  2:972515/1‑82 (MQ=255)

                                                                             tatCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaaaaa              >  1:4954884/1‑93 (MQ=255)

                                                                              atCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCga                         <  2:1709419/82‑1 (MQ=255)

                                                                               acATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTa            <  1:6852122/93‑1 (MQ=255)

                                                                                  ttGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATga         >  1:6375280/1‑93 (MQ=255)

                                                                                    gATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                   <  2:3826778/82‑1 (MQ=255)

                                                                                    aaTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTg       <  1:6949787/93‑1 (MQ=255)

                                                                                     aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                  <  2:2321619/82‑1 (MQ=255)

                                                                                     aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                  >  2:4897302/1‑82 (MQ=255)

                                                                                     aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGAACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGGAAAAAATTATGGTgt      >  1:7174906/1‑93 (MQ=255)

                                                                                      tGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAGGTGTCAGCGAGAGTaaaa                 <  2:7180411/82‑1 (MQ=255)

                                                                                      tGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaa                 <  2:2773948/82‑1 (MQ=255)

                                                                                      tGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaa                 <  2:7443819/82‑1 (MQ=255)

                                                                                      tGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaa                 <  2:5415005/82‑1 (MQ=255)

                                                                                       gTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaaa                >  2:6275793/1‑82 (MQ=255)

                                                                                       gTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaaa                <  2:845250/82‑1 (MQ=255)

                                                                                        ttCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaaaa               >  2:1272890/1‑82 (MQ=255)

                                                                                        ttCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCa   >  1:6665007/1‑94 (MQ=255)

                                                                                          cACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAAtt             <  2:2138174/82‑1 (MQ=255)

CGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAT  >  NC_000913/1650262‑1650444

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